### Finding Lyapunov functions in Sage

Is there a way to find Lyapunov functions candidates for a system of ODEs?

see for example: http://www.exampleproblems.com/wiki/index.php/ODELF1

I know I can simply "do the math" given I have a candidate function, but I wonder if there is any way to ~~quiclky ~~quickly explore potential candidates.