running radau5.f test step = 0, x = 0.00, y = 1.000000000000000E+00 0.000000000000000E+00 0.000000000000000E+00 step = 1, x = 0.00, y = 9.999999600000008E-01 3.999998320000048E-08 1.599999952000000E-14 step = 2, x = 0.00, y = 9.999996400000648E-01 3.599882716565646E-07 1.166354336752796E-11 step = 3, x = 0.00, y = 9.999970800042659E-01 2.913787382334434E-06 6.208351820680937E-09 step = 4, x = 0.00, y = 9.999863148909263E-01 1.307919521098008E-05 6.059138628105538E-07 step = 5, x = 0.00, y = 9.999755499712033E-01 2.135889538415425E-05 3.091133412644612E-06 step = 6, x = 0.00, y = 9.999647854413735E-01 2.724901208916897E-05 7.965546537469537E-06 step = 7, x = 0.00, y = 9.999478453625211E-01 3.255484201352810E-05 1.959979546541964E-05 step = 8, x = 0.00, y = 9.999253573022412E-01 3.530902486682297E-05 3.933367289202690E-05 step = 9, x = 0.00, y = 9.999028742076045E-01 3.614958793293454E-05 6.097620446265927E-05 step = 10, x = 0.00, y = 9.998490892191750E-01 3.647645985151893E-05 1.144343209735216E-04 step = 11, x = 0.01, y = 9.997567898163745E-01 3.647881802650697E-05 2.067313655991199E-04 step = 12, x = 0.01, y = 9.995023134527544E-01 3.643227719665359E-05 4.612542700490739E-04 step = 13, x = 0.05, y = 9.982091577136255E-01 3.619191783743737E-05 1.754650368537187E-03 step = 14, x = 0.18, y = 9.930472117083733E-01 3.525604364519484E-05 6.917532247981698E-03 step = 15, x = 0.30, y = 9.886740138499884E-01 3.447720471782070E-05 1.129150894529390E-02 Radau5: Radau method (Radau IIA) (implicit, order 5, embedded) Number of function evaluations = 88 Number of Jacobian evaluations = 8 Number of factorizations = 15 Number of lin sys solutions = 24 Number of performed steps = 17 Number of accepted steps = 15 Number of rejected steps = 1 Number of iterations (maximum) = 2 y = 9.886740138499884E-01 3.447720471782070E-05 1.129150894529390E-02 h = 8.160578540333708E-01