Recently I made my system identification algorithm SINDY to work. Not it can estimate a nonlinear model from measurement data that comes from a very nonlinear hydraulic system. The input signal increase with 1 step atvyhe time. At the top of the plot, the output signal won't increase so much compared to the beginning.
https://github.com/DanielMartensson/Mataveid
Now I have this model
$$\dot y = f (y, u) = 0.719906 + 2.533281y - 1.175728y^2 - 0.218748u^2 + 0.813850yu - 4.212704\sqrt y + 1.255115\sqrt u $$
I wonder how I can find a control law with lyapunov function for this nonlinear model.