10 November 2023 1:23:17.220 PM fsolve_be_test(): Fortran77 version Test fsolve_be() which finds a vector ym that satisfies the backward Euler residual: dydt(tm,ym) = (ym-yo)/(tm-to) predator_prey_test(): Initial ||ym-yo-(tm-to)*dydt(tm,ym)||: 9500.00 Final ||ym-yo-(tm-to)*dydt(tm,ym)||: 0.370114E-04 stiff_test(): Initial ||ym-yo-(tm-to)*dydt(tm,ym)||: 12.1114 fsolve_be() returned error flag info = 4 Final ||ym-yo-(tm-to)*dydt(tm,ym)||: 0.666134E-15 fsolve_be_test(): Normal end of execution. 10 November 2023 1:23:17.220 PM