Refer to Sections 2.3-2.5 of the lecture notes on Computational Physics for the theory, and definition of parameters. The integration has to be done using Euler method (see Section 1.1 of the lecture notes).
