pub fn implicit_euler(prob: &Problem, opts: &Options) -> Solution
Backward Euler method for stiff ODEs.
Uses fixed-point iteration to solve the implicit equation.