running 1 test step = 0, x = 0.0000, y1and5 = 0.000000000000000E+00 0.000000000000000E+00, diff1and5 = 0.0E+00 0.0E+00 step = 12, x = 0.0010, y1and5 = 1.916311926796952E-01 -2.688742328314606E+00, diff1and5 = 4.1E-08 2.5E-05 step = 16, x = 0.0020, y1and5 = 3.218545459857997E-01 -1.744730225557849E+00, diff1and5 = 8.7E-07 1.5E-04 step = 17, x = 0.0030, y1and5 = 3.335987191469630E-01 -7.596643438825829E-01, diff1and5 = 2.7E-06 1.6E-04 step = 20, x = 0.0040, y1and5 = 2.215636007078548E-01 -1.395458336432019E-01, diff1and5 = 3.1E-08 4.0E-06 step = 23, x = 0.0050, y1and5 = 2.822200240466802E-02 5.035341147888064E-02, diff1and5 = 1.2E-06 9.4E-06 step = 24, x = 0.0060, y1and5 = -1.725991658661986E-01 6.096730363234124E-02, diff1and5 = 4.3E-06 3.0E-04 step = 27, x = 0.0070, y1and5 = -3.052894121289478E-01 -4.074905119116261E-01, diff1and5 = 2.7E-07 2.1E-04 step = 30, x = 0.0080, y1and5 = -3.208689480273689E-01 -2.003438545924757E+00, diff1and5 = 5.2E-08 1.2E-04 step = 33, x = 0.0090, y1and5 = -2.112947302771765E-01 -2.771120812467848E+00, diff1and5 = 7.6E-06 9.8E-05 step = 35, x = 0.0100, y1and5 = -1.964755260832735E-02 -2.909025267207428E+00, diff1and5 = 1.2E-06 9.4E-06 step = 36, x = 0.0110, y1and5 = 1.805509844329725E-01 -2.399268748293620E+00, diff1and5 = 5.9E-06 3.9E-04 step = 38, x = 0.0120, y1and5 = 3.126188413107825E-01 -1.460587294362982E+00, diff1and5 = 4.7E-07 3.3E-04 step = 39, x = 0.0130, y1and5 = 3.258822671433425E-01 -4.922170285774613E-01, diff1and5 = 5.6E-06 1.3E-04 step = 43, x = 0.0140, y1and5 = 2.151031987252550E-01 1.085467149160739E-01, diff1and5 = 3.4E-07 1.6E-05 step = 45, x = 0.0150, y1and5 = 2.281762592295544E-02 2.830230580502269E-01, diff1and5 = 1.7E-06 2.6E-05 step = 46, x = 0.0160, y1and5 = -1.771028609942352E-01 2.901769168707256E-01, diff1and5 = 5.2E-06 3.0E-04 step = 48, x = 0.0170, y1and5 = -3.089261062494018E-01 -1.143019737034166E-01, diff1and5 = 2.9E-07 1.5E-04 step = 53, x = 0.0180, y1and5 = -3.239999209426208E-01 -1.773489735052747E+00, diff1and5 = 1.7E-07 1.5E-05 step = 55, x = 0.0190, y1and5 = -2.139062023861129E-01 -2.549976140009503E+00, diff1and5 = 2.4E-07 4.3E-04 step = 57, x = 0.0200, y1and5 = -2.183476092114947E-02 -2.694693795055462E+00, diff1and5 = 7.8E-07 1.6E-05 step = 58, x = 0.0210, y1and5 = 1.787239963326930E-01 -2.190379451978742E+00, diff1and5 = 5.6E-06 8.5E-05 step = 60, x = 0.0220, y1and5 = 3.110923612041102E-01 -1.257994974815820E+00, diff1and5 = 4.1E-07 3.1E-05 step = 61, x = 0.0230, y1and5 = 3.246041884373977E-01 -2.951538111326443E-01, diff1and5 = 3.2E-06 3.3E-04 step = 64, x = 0.0240, y1and5 = 2.140352114671024E-01 2.997035675328908E-01, diff1and5 = 7.0E-08 5.5E-06 step = 66, x = 0.0250, y1and5 = 2.192275843691317E-02 4.688416135226632E-01, diff1and5 = 5.5E-08 6.9E-07 step = 68, x = 0.0260, y1and5 = -1.778526757699608E-01 4.723850877530753E-01, diff1and5 = 3.3E-08 5.7E-06 step = 71, x = 0.0270, y1and5 = -3.095275377534654E-01 7.592242763952231E-02, diff1and5 = 7.0E-08 3.2E-05 step = 75, x = 0.0280, y1and5 = -3.245175839432364E-01 -1.595988841865317E+00, diff1and5 = 1.5E-07 5.8E-05 step = 77, x = 0.0290, y1and5 = -2.143416391952930E-01 -2.377050934355847E+00, diff1and5 = 2.9E-06 7.8E-05 step = 79, x = 0.0300, y1and5 = -2.219528408240040E-02 -2.524853935988888E+00, diff1and5 = 2.2E-07 8.5E-06 step = 80, x = 0.0310, y1and5 = 1.784315561831284E-01 -2.024353412245143E+00, diff1and5 = 4.0E-06 2.9E-04 step = 82, x = 0.0320, y1and5 = 3.108403881688425E-01 -1.095033084537755E+00, diff1and5 = 7.0E-07 6.6E-05 step = 84, x = 0.0330, y1and5 = 3.243889870027998E-01 -1.351281531066862E-01, diff1and5 = 1.1E-06 5.7E-05 step = 86, x = 0.0340, y1and5 = 2.138589333550945E-01 4.561261291820525E-01, diff1and5 = 3.8E-07 1.3E-05 step = 89, x = 0.0350, y1and5 = 2.177505550254943E-02 6.220899930796396E-01, diff1and5 = 7.3E-08 2.5E-06 step = 90, x = 0.0360, y1and5 = -1.780029402927223E-01 6.219615062635127E-01, diff1and5 = 2.7E-05 8.8E-04 step = 93, x = 0.0370, y1and5 = -3.096268590211072E-01 2.255017983309159E-01, diff1and5 = 1.9E-07 8.0E-05 step = 99, x = 0.0380, y1and5 = -3.246033030857612E-01 -1.450729771171510E+00, diff1and5 = 3.2E-07 6.5E-05 step = 102, x = 0.0390, y1and5 = -2.144090078932315E-01 -2.234602271558864E+00, diff1and5 = 1.3E-06 8.2E-05 step = 103, x = 0.0400, y1and5 = -2.225575794426700E-02 -2.385080348224730E+00, diff1and5 = 4.7E-07 9.6E-06 step = 104, x = 0.0410, y1and5 = 1.783880058813793E-01 -1.886779093725898E+00, diff1and5 = 1.0E-05 1.3E-04 step = 106, x = 0.0420, y1and5 = 3.107983906300928E-01 -9.604321864034761E-01, diff1and5 = 4.7E-07 2.9E-05 step = 108, x = 0.0430, y1and5 = 3.243575044289585E-01 -3.145827679016736E-03, diff1and5 = 2.5E-06 9.6E-05 step = 111, x = 0.0440, y1and5 = 2.138294790615231E-01 5.856716367852630E-01, diff1and5 = 2.5E-07 1.6E-05 step = 114, x = 0.0450, y1and5 = 2.175063895712489E-02 7.492442901224280E-01, diff1and5 = 1.9E-07 1.1E-06 step = 115, x = 0.0460, y1and5 = -1.780245999100694E-01 7.466802560157104E-01, diff1and5 = 2.8E-05 9.8E-04 step = 118, x = 0.0470, y1and5 = -3.096432354425651E-01 3.483520547215561E-01, diff1and5 = 3.3E-07 9.6E-05 step = 123, x = 0.0480, y1and5 = -3.246172295071914E-01 -1.330381955254372E+00, diff1and5 = 3.8E-08 4.6E-05 step = 126, x = 0.0490, y1and5 = -2.144228860345288E-01 -2.116620164596746E+00, diff1and5 = 6.7E-07 2.0E-04 step = 127, x = 0.0500, y1and5 = -2.226517863103583E-02 -2.269179795792819E+00, diff1and5 = 3.9E-08 8.7E-06 Radau5: Radau method (Radau IIA) (implicit, order 5, embedded) Number of function evaluations = 1511 Number of Jacobian evaluations = 126 Number of factorizations = 166 Number of lin sys solutions = 461 Number of performed steps = 166 Number of accepted steps = 127 Number of rejected steps = 6 Number of iterations (maximum) = 5 y1to3 = -2.226517863103583E-02 3.068700099790518E+00 2.898340496735758E+00 y4to5 = 2.033525394061642E+00 -2.269179795792819E+00 h = 7.789307871956906E-04