runge_kutta: Runge-Kutta for the calculation of the a and b vectors and the c matrix in a EM step
Description
Performs the Runge-Kutta method of fourth order.
Usage
runge_kutta(avector, bvector, cmatrix, dt, h, S, s)
Arguments
- avector
The a vector.
- bvector
The b vector.
- cmatrix
The c matrix.
- dt
The increment.
- h
Step-length.
- S
Sub-intensity matrix.
- s
Exit rates.