runge_kutta {matrixdist}

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.