pub fn solve_row1( f: CModule, x_span: Tensor, y0: Tensor, step: Tensor, ) -> Result<(Tensor, Tensor)>
Solves ODE using first-order Rosenbrock method (Row1)