A new method, called the intrinsic method, is proposed to derive suitable Lyapunov functions for a general class of non-linear systems expressed in state variables as n first-order nonlinear differential equations. This method, which applies the integration-by-parts procedure, derives a Lyapunov function directly from the differential equations under study.