runge_kutta function

Runge-Kutta for the calculation of the a and b vectors and the c matrix in a EM step