Numerical solution by Runge Kutta Method