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