The Euler, the Leapfrog and the fourth-order Runge-Kutta method are used to solve ODEs in MATLAB. Stability and accuracy improved for different numerical methods. The false position method and the WKB approximation are used to estimate the eigenvalues for a simple harmonic oscillator.