powered by
Performce the RK of forth order
runge_kutta(avector, bvector, cmatrix, dt, h, S, t)
the a vector
the b vector
the c matrix
the increment
step-length
sub-intensity
exit rates