runge_kutta: Runge Kutta for the calculation of the a,b and c vectors in a EM step
Description
Performce the RK of forth order
Usage
runge_kutta(avector, bvector, cmatrix, dt, h, S, t)
Arguments
- avector
the a vector
- bvector
the b vector
- cmatrix
the c matrix
- dt
the increment
- h
step-length
- S
sub-intensity
- t
exit rates