pub fn trbdf2(prob: &Problem, opts: &Options) -> Solution
TR-BDF2 method: two-stage implicit method combining trapezoidal rule with BDF2.