The problem of constructing a Lyapunov function for continuous-time nonlinear dynamical systems is tackled in this paper via a sampling-based approach. The main idea of the sampling-based method is to verify a Lyapunov-type inequality for a finite number of points (known state vectors) in the state-space and then to extend the validity of the Lyapunov inequality to a neighborhood around these points. In this way, the validity of a Lyapunov function candidate can be certified for a region of interest in the state-space in a systematic way. A candidate Lyapunov function is computed for each sample point using a recent converse Lyapunov theorem for continuous-time nonlinear systems. For certifying the candidate Lyapunov function on a neighborh...