Abstract Discovering governing equations of complex dynamical systems directly from data is a central problem in scientific machine learning. In recent years, the sparse identification nonlinear dynamics (SINDy) framework, powered by heuristic regression methods, has become dominant tool for learning parsimonious models. We propose an exact formulation SINDy using mixed-integer optimization (MI...