powered by
This function calculates the product integral of a matrix function from s to t. It uses a Runge-Kutta method implemented in C++ for efficiency.
prodint(A, s, t, n)
A matrix (transition probabilities if A is an intensity matrix)
A function returning a matrix (intensity matrix)
Initial timepoint
End timepoint
Number of steps for the Runge-Kutta algorithm