mbd_solve_methode: mbd ODE system integrator for a specific ode method
In Giappo/mbd: Multiple Birth Death Diversification

Integrates "func" in the time interval

mbd_solve_methode(
  vector,
  time_interval,
  func = mbd::mbd_loglik_rhs,
  parms,
  methode = "lsoda"
)

`vector`	a vector
`time_interval`	a time interval
`func`	function for the right hand side of the ODE
`parms`	some parameters to pass on the ode function
`methode`	specifies how the integration must be performed: `lsodes`: use `lsodes` and `deSolve::ode` `ode45`: use `ode45` and `deSolve::ode` `lsoda`: use `lsoda` and `deSolve::ode`