logLikelihoodPH_RKs function

Loglikelihood of PI with phase-type using Runge-Kutta