Changed runge-kutta method to make him more compact, i.e. the system
notation has been used (but still problems in the graph)