Function solve_ode 

pub fn solve_ode<const N: usize, const M: usize, Sist: ODESystem<N> + Clone>(
    sist: Sist,
    odeparam: ODEParam,
    odesolver: ODESolver,
) -> (Data<M>, Sist, ODEParam)

Main function for solving ODEs, returns an tuple with (Data, Sist, ODEParam) where Sist and ODEParam have updated values.

inputs

sist : Sist. Any type which implements ODESystem<N> and Clone

odeparam : ODEParam. An value of ODEParam.

odesolver : ODESolver. A choice of an ODE solver.