Automatic learning of (linear and nonlinear) control functions for nonlinear dynamical systems with stability certificates. The stability of the closed-loop system is certified by means of a Lyapunov Function, whose correctness is guaranteed via SMT solvers.