a_rungekutta function

Runge-Kutta for the calculation of the a vector in a EM step