It is often good to determine whether or not an equilibrium point of a differential equation model is globally stable. To do this, a Lyapunov function is often employed. May I know if there are methods of determining these Lyapunov functions or one will have to be doing trial and error until one gets a good function.