Function solve_glrk4

pub fn solve_glrk4(
    f: CModule,
    x_span: Tensor,
    y0: Tensor,
    step: Tensor,
    optimizer: &dyn Optimizer,
) -> Result<(Tensor, Tensor)>

Solves ODE using Gauss-Legendre-Runge-Kutta 4th order method