Liouville–Green method