This post shows how to utilize sum-of-squares optimization to find Lyapunov function to prove global or local stability.