running dop853.f test step(A) = 1, err = 1.959638663289740E-12, h_new = 6.000000000000001E-04, n_yes = 0, n_no = 1, h*lambda = 0.000000000000000E+00 step(A) = 2, err = 2.955980245499226E-06, h_new = 2.651893439400644E-03, n_yes = 0, n_no = 2, h*lambda = 5.999638637675994E-01 step(A) = 3, err = 3.289190833225518E-01, h_new = 2.742599827984791E-03, n_yes = 0, n_no = 3, h*lambda = 2.647684181446536E+00 step(A) = 4, err = 2.963282904871433E-02, h_new = 3.832068947656679E-03, n_yes = 0, n_no = 4, h*lambda = 2.731697216875155E+00 step(A) = 5, err = 1.765232930466588E-02, h_new = 5.712491558181956E-03, n_yes = 0, n_no = 5, h*lambda = 3.803790881529834E+00 step(A) = 6, err = 1.357237168717890E-08, h_new = 3.427494934909174E-02, n_yes = 0, n_no = 6, h*lambda = 5.641296505768319E+00 step(R) = 7, err = 3.049983581528165E+02, h_new = 1.508975343207560E-02 step(R) = 8, err = 2.753687427000322E+01, h_new = 8.972951547905639E-03 step(A) = 9, err = 1.190643609945584E-01, h_new = 1.053671198224488E-02, n_yes = 1, n_no = 0, h*lambda = 8.789351829648975E+00 step(R) = 10, err = 3.891356809844428E+00, h_new = 6.814205085743480E-03 step(A) = 11, err = 1.254149001555174E-01, h_new = 7.949943966746855E-03, n_yes = 2, n_no = 0, h*lambda = 6.633281260805703E+00 step(A) = 12, err = 1.807919811765532E-01, h_new = 7.594694594261238E-03, n_yes = 3, n_no = 0, h*lambda = 6.591787703063471E+00 step(R) = 13, err = 1.084024659426730E+00, h_new = 6.766637768784153E-03 step(A) = 14, err = 2.161217856127681E-01, h_new = 7.375273321472055E-03, n_yes = 4, n_no = 0, h*lambda = 6.504692821305968E+00 step(A) = 15, err = 2.488521029659334E-01, h_new = 7.246407845260238E-03, n_yes = 5, n_no = 0, h*lambda = 6.463674121840817E+00 step(A) = 16, err = 7.435797934263277E-01, h_new = 6.767828353418105E-03, n_yes = 6, n_no = 0, h*lambda = 6.874624514790060E+00 step(A) = 17, err = 5.600360888226632E-01, h_new = 6.548846985852363E-03, n_yes = 7, n_no = 0, h*lambda = 6.379539619848891E+00 step(A) = 18, err = 2.689360370853225E-01, h_new = 6.945463407347431E-03, n_yes = 8, n_no = 0, h*lambda = 6.134291480727791E+00 step(A) = 19, err = 4.737739730192550E-01, h_new = 6.862736705124246E-03, n_yes = 9, n_no = 0, h*lambda = 6.462359077921882E+00 step(A) = 20, err = 4.081931733449650E-01, h_new = 6.908466077183188E-03, n_yes = 10, n_no = 0, h*lambda = 6.342614488854702E+00 step(A) = 21, err = 3.936599445388877E-01, h_new = 6.986086931314706E-03, n_yes = 11, n_no = 0, h*lambda = 6.341752994061589E+00 step(A) = 22, err = 4.078712131210029E-01, h_new = 7.033331916449623E-03, n_yes = 12, n_no = 0, h*lambda = 6.368533582524074E+00 step(A) = 23, err = 4.114534445404955E-01, h_new = 7.073160854625421E-03, n_yes = 13, n_no = 0, h*lambda = 6.366738680516228E+00 step(A) = 24, err = 4.033631445442230E-01, h_new = 7.130894601877946E-03, n_yes = 14, n_no = 0, h*lambda = 6.357033171613826E+00 THE PROBLEM SEEMS TO BECOME STIFF AT X = 1.563905377322407E-02, ACCEPTED STEP = 21 step(A) = 25, err = 4.077460795117875E-01, h_new = 7.179394245392615E-03, n_yes = 15, n_no = 0, h*lambda = 6.362627602473077E+00 step(A) = 26, err = 4.037854318248442E-01, h_new = 7.237048469636972E-03, n_yes = 16, n_no = 0, h*lambda = 6.358575193217763E+00 step(A) = 27, err = 4.074225933954899E-01, h_new = 7.286992996411961E-03, n_yes = 17, n_no = 0, h*lambda = 6.361764246185123E+00 step(A) = 28, err = 4.025036243432588E-01, h_new = 7.348431279766412E-03, n_yes = 18, n_no = 0, h*lambda = 6.356714697192896E+00 step(A) = 29, err = 4.066531572312910E-01, h_new = 7.400893045571465E-03, n_yes = 19, n_no = 0, h*lambda = 6.360744201444202E+00 step(A) = 30, err = 4.015562150819599E-01, h_new = 7.465490436421224E-03, n_yes = 20, n_no = 0, h*lambda = 6.355444820557152E+00 step(A) = 31, err = 4.058232560149916E-01, h_new = 7.520708164291287E-03, n_yes = 21, n_no = 0, h*lambda = 6.359535584820661E+00 step(A) = 32, err = 4.004792807135695E-01, h_new = 7.588898415478732E-03, n_yes = 22, n_no = 0, h*lambda = 6.353975475377664E+00 step(A) = 33, err = 4.049042413246206E-01, h_new = 7.647195765549171E-03, n_yes = 23, n_no = 0, h*lambda = 6.358248130028030E+00 step(A) = 34, err = 3.993867275576406E-01, h_new = 7.719168373327975E-03, n_yes = 24, n_no = 0, h*lambda = 6.352456077046539E+00 step(A) = 35, err = 4.038828131676189E-01, h_new = 7.780922719248015E-03, n_yes = 25, n_no = 0, h*lambda = 6.356790279382604E+00 step(A) = 36, err = 3.981946527794150E-01, h_new = 7.857089195073108E-03, n_yes = 26, n_no = 0, h*lambda = 6.350804544551124E+00 step(A) = 37, err = 4.027908750123589E-01, h_new = 7.922627551144631E-03, n_yes = 27, n_no = 0, h*lambda = 6.355226958151572E+00 step(A) = 38, err = 3.968989618768802E-01, h_new = 8.003441114451426E-03, n_yes = 28, n_no = 0, h*lambda = 6.349002719814372E+00 step(A) = 39, err = 4.016091857812493E-01, h_new = 8.073164627669244E-03, n_yes = 29, n_no = 0, h*lambda = 6.353526548383705E+00 step(A) = 40, err = 3.954854478243182E-01, h_new = 8.159151639133593E-03, n_yes = 30, n_no = 0, h*lambda = 6.347032595169746E+00 step(A) = 41, err = 4.003287180194667E-01, h_new = 8.233517650455054E-03, n_yes = 31, n_no = 0, h*lambda = 6.351673639113059E+00 step(A) = 42, err = 3.939366752571402E-01, h_new = 8.325294943726226E-03, n_yes = 32, n_no = 0, h*lambda = 6.344867851517789E+00 step(A) = 43, err = 3.989362262003876E-01, h_new = 8.404835221718525E-03, n_yes = 33, n_no = 0, h*lambda = 6.349646213311302E+00 step(A) = 44, err = 3.922322406642030E-01, h_new = 8.503129665852644E-03, n_yes = 34, n_no = 0, h*lambda = 6.342478281345484E+00 step(A) = 45, err = 3.974165319398817E-01, h_new = 8.588465393941008E-03, n_yes = 35, n_no = 0, h*lambda = 6.347418360429272E+00 step(A) = 46, err = 3.903471745854931E-01, h_new = 8.694141405485170E-03, n_yes = 36, n_no = 0, h*lambda = 6.339826580596097E+00 step(A) = 47, err = 3.957515763723521E-01, h_new = 8.786003607233700E-03, n_yes = 37, n_no = 0, h*lambda = 6.344958758570509E+00 step(A) = 48, err = 3.882509034593817E-01, h_new = 8.900098790260436E-03, n_yes = 38, n_no = 0, h*lambda = 6.336866791986790E+00 step(A) = 49, err = 3.939197892963724E-01, h_new = 8.999354543786694E-03, n_yes = 39, n_no = 0, h*lambda = 6.342229201048759E+00 step(A) = 50, err = 3.859055072050792E-01, h_new = 9.123127602653778E-03, n_yes = 40, n_no = 0, h*lambda = 6.333541511512307E+00 step(A) = 51, err = 3.918952251806554E-01, h_new = 9.230814257720642E-03, n_yes = 41, n_no = 0, h*lambda = 6.339182626112716E+00 step(A) = 52, err = 3.832633846547456E-01, h_new = 9.365810261080880E-03, n_yes = 42, n_no = 0, h*lambda = 6.329778111191833E+00 step(A) = 53, err = 3.896464417389354E-01, h_new = 9.483180686680590E-03, n_yes = 43, n_no = 0, h*lambda = 6.335760400824825E+00 step(A) = 54, err = 3.802639546871396E-01, h_new = 9.631321714049920E-03, n_yes = 44, n_no = 0, h*lambda = 6.325483341323892E+00 step(A) = 55, err = 3.871350479976269E-01, h_new = 9.759904944406365E-03, n_yes = 45, n_no = 0, h*lambda = 6.331888556119084E+00 step(A) = 56, err = 3.768289131875920E-01, h_new = 9.923618730235782E-03, n_yes = 46, n_no = 0, h*lambda = 6.320535467755772E+00 step(A) = 57, err = 3.843138180202694E-01, h_new = 1.006530245373337E-02, n_yes = 47, n_no = 0, h*lambda = 6.327472463992327E+00 step(A) = 58, err = 3.728552563616889E-01, h_new = 1.024770949226123E-02, n_yes = 48, n_no = 0, h*lambda = 6.314772533647831E+00 step(A) = 59, err = 3.811242574519738E-01, h_new = 1.040485400521445E-02, n_yes = 49, n_no = 0, h*lambda = 6.322389158961099E+00 step(A) = 60, err = 3.682047477363228E-01, h_new = 1.061004747695080E-02, n_yes = 50, n_no = 0, h*lambda = 6.307974325744030E+00 step(A) = 61, err = 3.774935233222885E-01, h_new = 1.078564566608427E-02, n_yes = 51, n_no = 0, h*lambda = 6.316476004708351E+00 step(A) = 62, err = 3.626875322468118E-01, h_new = 1.101912423933182E-02, n_yes = 52, n_no = 0, h*lambda = 6.299833751000393E+00 step(A) = 63, err = 3.733307278971418E-01, h_new = 1.121702975772993E-02, n_yes = 53, n_no = 0, h*lambda = 6.309513557016553E+00 step(A) = 64, err = 3.560357713035250E-01, h_new = 1.148639315495337E-02, n_yes = 54, n_no = 0, h*lambda = 6.289909649186166E+00 step(A) = 65, err = 3.685230977022650E-01, h_new = 1.171165033362278E-02, n_yes = 55, n_no = 0, h*lambda = 6.301198914779232E+00 step(A) = 66, err = 3.478595229546738E-01, h_new = 1.202777013980060E-02, n_yes = 56, n_no = 0, h*lambda = 6.277545458237089E+00 step(A) = 67, err = 3.629338985345094E-01, h_new = 1.228709419065767E-02, n_yes = 57, n_no = 0, h*lambda = 6.291103027921865E+00 step(A) = 68, err = 3.375694802238602E-01, h_new = 1.266619879821853E-02, n_yes = 58, n_no = 0, h*lambda = 6.261721378448611E+00 step(A) = 69, err = 3.564085412695703E-01, h_new = 1.296866573696537E-02, n_yes = 59, n_no = 0, h*lambda = 6.278600116502583E+00 step(A) = 70, err = 3.242343501050571E-01, h_new = 1.343632277479681E-02, n_yes = 60, n_no = 0, h*lambda = 6.240767851985754E+00 step(A) = 71, err = 3.488100813437398E-01, h_new = 1.379428863529635E-02, n_yes = 61, n_no = 0, h*lambda = 6.262747484631655E+00 step(A) = 72, err = 3.063010456581432E-01, h_new = 1.439372690752747E-02, n_yes = 62, n_no = 0, h*lambda = 6.211764538341034E+00 step(A) = 73, err = 3.401584201599718E-01, h_new = 1.482366578386025E-02, n_yes = 63, n_no = 0, h*lambda = 6.242077724426380E+00 step(A) = 74, err = 2.810069968550686E-01, h_new = 1.563538150357789E-02, n_yes = 64, n_no = 0, h*lambda = 6.169147305416486E+00 step(A) = 75, err = 3.311667961548223E-01, h_new = 1.615642033930710E-02, n_yes = 65, n_no = 0, h*lambda = 6.214256864081382E+00 step(A) = 76, err = 2.430597376241931E-01, h_new = 1.735295837251806E-02, n_yes = 66, n_no = 0, h*lambda = 6.101034621215406E+00 step(A) = 77, err = 3.256047954119587E-01, h_new = 1.796923901210033E-02, n_yes = 67, n_no = 0, h*lambda = 6.175753177258414E+00 step(A) = 78, err = 1.815900523460503E-01, h_new = 2.001638695004108E-02, n_yes = 67, n_no = 1, h*lambda = 5.977717665587894E+00 step(A) = 79, err = 3.447500895029070E-01, h_new = 2.057975257528611E-02, n_yes = 68, n_no = 0, h*lambda = 6.124456071176638E+00 step(A) = 80, err = 7.494034684253342E-02, h_new = 2.560608578640618E-02, n_yes = 68, n_no = 1, h*lambda = 5.705940812368194E+00 step(A) = 81, err = 6.648330029613858E-01, h_new = 2.425194833240622E-02, n_yes = 69, n_no = 0, h*lambda = 6.133098554590505E+00 step(A) = 82, err = 5.871646231207660E-02, h_new = 3.110958076743863E-02, n_yes = 69, n_no = 1, h*lambda = 4.867384689284566E+00 step(A) = 83, err = 1.224305409033219E-01, h_new = 3.640411243348058E-02, n_yes = 69, n_no = 2, h*lambda = 4.522823329851793E+00 step(A) = 84, err = 7.040557969695867E-02, h_new = 4.565013556890047E-02, n_yes = 70, n_no = 0, h*lambda = 7.528227971169771E+00 step(R) = 85, err = 3.615022310303625E+02, h_new = 1.967527033979558E-02 step(A) = 86, err = 1.763775503304252E-01, h_new = 2.199671381094752E-02, n_yes = 70, n_no = 1, h*lambda = 1.307303765126028E+00 step(R) = 87, err = 1.058639370802532E+03, h_new = 7.414292453421656E-03 step(A) = 88, err = 8.193025328826009E-03, h_new = 1.216554823051712E-02, n_yes = 70, n_no = 2, h*lambda = 1.075904249959075E+00 step(A) = 89, err = 3.247216552609690E-01, h_new = 7.680213897570433E-03, n_yes = 70, n_no = 3, h*lambda = 2.517281269948783E+00 step(R) = 90, err = 4.696476076101794E+02, h_new = 3.203644448936211E-03 step(A) = 91, err = 5.081123195791324E-01, h_new = 3.137919907885860E-03, n_yes = 70, n_no = 4, h*lambda = 9.330798214902337E-01 step(R) = 92, err = 1.109363451506409E+01, h_new = 2.090499808978971E-03 step(A) = 93, err = 2.131027857985464E-02, h_new = 3.043820241309415E-03, n_yes = 70, n_no = 5, h*lambda = 2.978925948900621E+00 step(R) = 94, err = 1.008394292326365E+00, h_new = 1.879484914696249E-03 step(A) = 95, err = 5.003903543931383E-01, h_new = 1.844453613420861E-03, n_yes = 0, n_no = 6, h*lambda = 1.780820386221682E+00 step(A) = 96, err = 3.285995718087074E-03, h_new = 3.392554809898990E-03, n_yes = 0, n_no = 7, h*lambda = 1.841306956507960E+00 step(R) = 97, err = 8.454091748140812E+00, h_new = 2.338223455082580E-03 step(A) = 98, err = 4.898107639146253E-01, h_new = 2.300779406584214E-03, n_yes = 0, n_no = 8, h*lambda = 2.364459947871286E+00 step(A) = 99, err = 1.400696001087453E-01, h_new = 2.647430115769552E-03, n_yes = 0, n_no = 9, h*lambda = 2.323978779541553E+00 step(A) = 100, err = 4.137503866928380E-02, h_new = 3.547927304229437E-03, n_yes = 0, n_no = 10, h*lambda = 2.668153875692219E+00 step(A) = 101, err = 2.165847585282912E-02, h_new = 5.155416245600121E-03, n_yes = 0, n_no = 11, h*lambda = 3.564561766954261E+00 step(A) = 102, err = 2.127984117506283E-03, h_new = 1.001175463256056E-02, n_yes = 0, n_no = 12, h*lambda = 5.156008614706338E+00 step(A) = 103, err = 4.021520223977747E-01, h_new = 1.009726896079593E-02, n_yes = 1, n_no = 0, h*lambda = 9.923844917727425E+00 step(R) = 104, err = 5.101997776428038E+01, h_new = 5.558772508811375E-03 step(A) = 105, err = 1.004589778260953E-02, h_new = 8.891454616236583E-03, n_yes = 1, n_no = 1, h*lambda = 5.482430585146485E+00 step(A) = 106, err = 2.275060301211081E-03, h_new = 1.070526189573793E-02, n_yes = 1, n_no = 2, h*lambda = 5.454889805521367E+00 step(R) = 107, err = 2.885978696859342E+00, h_new = 8.439236657067955E-03 step(A) = 108, err = 2.070056790907456E-01, h_new = 9.248001524194572E-03, n_yes = 2, n_no = 0, h*lambda = 8.217908845733771E+00 step(R) = 109, err = 3.266643901495132E+00, h_new = 6.550628206108270E-03 step(A) = 110, err = 1.083509983474625E-01, h_new = 7.783435765718324E-03, n_yes = 3, n_no = 0, h*lambda = 6.340484772404460E+00 step(A) = 111, err = 9.254395628318393E-02, h_new = 7.938381250603994E-03, n_yes = 4, n_no = 0, h*lambda = 6.302011941423816E+00 step(R) = 112, err = 1.160859196036360E+00, h_new = 7.012566921305652E-03 step(A) = 113, err = 2.256803420494697E-01, h_new = 7.602086585789731E-03, n_yes = 5, n_no = 0, h*lambda = 6.702343634595116E+00 step(A) = 114, err = 3.633185222355420E-01, h_new = 7.162812549441357E-03, n_yes = 6, n_no = 0, h*lambda = 6.658065404312254E+00 step(A) = 115, err = 7.447589964965239E-01, h_new = 6.688429065786689E-03, n_yes = 7, n_no = 0, h*lambda = 6.754678647746046E+00 step(A) = 116, err = 4.270349561007837E-01, h_new = 6.695125258492039E-03, n_yes = 8, n_no = 0, h*lambda = 6.266807641803351E+00 step(A) = 117, err = 3.241872230025924E-01, h_new = 6.936681101854841E-03, n_yes = 9, n_no = 0, h*lambda = 6.232733381749739E+00 step(A) = 118, err = 4.312213324972082E-01, h_new = 6.935163576632131E-03, n_yes = 10, n_no = 0, h*lambda = 6.413926107808549E+00 step(A) = 119, err = 4.187789015969542E-01, h_new = 6.959068626637126E-03, n_yes = 11, n_no = 0, h*lambda = 6.369070235207361E+00 step(A) = 120, err = 3.975480475475240E-01, h_new = 7.028617761148532E-03, n_yes = 12, n_no = 0, h*lambda = 6.346897025088153E+00 step(A) = 121, err = 4.084659132486355E-01, h_new = 7.074861747009606E-03, n_yes = 13, n_no = 0, h*lambda = 6.365529555499854E+00 step(A) = 122, err = 4.060100732388404E-01, h_new = 7.126780218799708E-03, n_yes = 14, n_no = 0, h*lambda = 6.361635455485464E+00 THE PROBLEM SEEMS TO BECOME STIFF AT X = 8.759592223459979E-01, ACCEPTED STEP = 109 step(A) = 123, err = 4.079019965058998E-01, h_new = 7.174908986936162E-03, n_yes = 15, n_no = 0, h*lambda = 6.362085653492173E+00 step(A) = 124, err = 4.032567986990441E-01, h_new = 7.233711661577681E-03, n_yes = 16, n_no = 0, h*lambda = 6.357788804581586E+00 step(A) = 125, err = 4.075247138426195E-01, h_new = 7.283404987172566E-03, n_yes = 17, n_no = 0, h*lambda = 6.362071764095525E+00 step(A) = 126, err = 4.026373250609086E-01, h_new = 7.344508107879011E-03, n_yes = 18, n_no = 0, h*lambda = 6.356878640375506E+00 step(A) = 127, err = 4.066441324506645E-01, h_new = 7.396962385589742E-03, n_yes = 19, n_no = 0, h*lambda = 6.360703437525811E+00 step(A) = 128, err = 4.015773441125634E-01, h_new = 7.461476393626732E-03, n_yes = 20, n_no = 0, h*lambda = 6.355482753882741E+00 step(A) = 129, err = 4.058625590655167E-01, h_new = 7.516573440539908E-03, n_yes = 21, n_no = 0, h*lambda = 6.359594648476476E+00 step(A) = 130, err = 4.005156654303061E-01, h_new = 7.584640069767726E-03, n_yes = 22, n_no = 0, h*lambda = 6.354023251466810E+00 step(A) = 131, err = 4.049351179411302E-01, h_new = 7.642831857938843E-03, n_yes = 23, n_no = 0, h*lambda = 6.358291792682883E+00 step(A) = 132, err = 3.994252920573956E-01, h_new = 7.714670282761915E-03, n_yes = 24, n_no = 0, h*lambda = 6.352510089445839E+00 step(A) = 133, err = 4.039185189030711E-01, h_new = 7.776302712549671E-03, n_yes = 25, n_no = 0, h*lambda = 6.356841223960268E+00 step(A) = 134, err = 3.982362204651791E-01, h_new = 7.852321505048351E-03, n_yes = 26, n_no = 0, h*lambda = 6.350862151023890E+00 step(A) = 135, err = 4.028289394840012E-01, h_new = 7.917726566294398E-03, n_yes = 27, n_no = 0, h*lambda = 6.355281614091100E+00 step(A) = 136, err = 3.969443649199136E-01, h_new = 7.998375772334330E-03, n_yes = 28, n_no = 0, h*lambda = 6.349065935360435E+00 step(A) = 137, err = 4.016504420881893E-01, h_new = 8.067951562448000E-03, n_yes = 29, n_no = 0, h*lambda = 6.353586059912446E+00 step(A) = 138, err = 3.955350426202229E-01, h_new = 8.153755244251135E-03, n_yes = 30, n_no = 0, h*lambda = 6.347101808549010E+00 step(A) = 139, err = 4.003734967315674E-01, h_new = 8.227957034020754E-03, n_yes = 31, n_no = 0, h*lambda = 6.351738620835539E+00 step(A) = 140, err = 3.939911328251167E-01, h_new = 8.319528592135658E-03, n_yes = 32, n_no = 0, h*lambda = 6.344944073888236E+00 step(A) = 141, err = 3.989850010577038E-01, h_new = 8.398885426613493E-03, n_yes = 33, n_no = 0, h*lambda = 6.349717450916532E+00 step(A) = 142, err = 3.922923068837588E-01, h_new = 8.496947646654797E-03, n_yes = 34, n_no = 0, h*lambda = 6.342562622246558E+00 step(A) = 143, err = 3.974698582768102E-01, h_new = 8.582077396113908E-03, n_yes = 35, n_no = 0, h*lambda = 6.347496810177391E+00 step(A) = 144, err = 3.904137727186877E-01, h_new = 8.687489546571957E-03, n_yes = 36, n_no = 0, h*lambda = 6.339920422762179E+00 step(A) = 145, err = 3.958101134434840E-01, h_new = 8.779119156380757E-03, n_yes = 37, n_no = 0, h*lambda = 6.345045574618171E+00 step(A) = 146, err = 3.883251695256497E-01, h_new = 8.892912321973546E-03, n_yes = 38, n_no = 0, h*lambda = 6.336971848328186E+00 step(A) = 147, err = 3.939843275326930E-01, h_new = 8.991903794032542E-03, n_yes = 39, n_no = 0, h*lambda = 6.342325798329268E+00 step(A) = 148, err = 3.859888573585522E-01, h_new = 9.115328304089620E-03, n_yes = 40, n_no = 0, h*lambda = 6.333659933820814E+00 step(A) = 149, err = 3.919667192253866E-01, h_new = 9.222712601171424E-03, n_yes = 41, n_no = 0, h*lambda = 6.339290757261591E+00 step(A) = 150, err = 3.833576079253681E-01, h_new = 9.357302597411634E-03, n_yes = 42, n_no = 0, h*lambda = 6.329912639361193E+00 step(A) = 151, err = 3.897260517275788E-01, h_new = 9.474324461901461E-03, n_yes = 43, n_no = 0, h*lambda = 6.335882260436112E+00 step(A) = 152, err = 3.803713435484557E-01, h_new = 9.621987520663242E-03, n_yes = 44, n_no = 0, h*lambda = 6.325637519545015E+00 step(A) = 153, err = 3.872241957887284E-01, h_new = 9.750165509616092E-03, n_yes = 45, n_no = 0, h*lambda = 6.332026930594409E+00 step(A) = 154, err = 3.769524534465372E-01, h_new = 9.913309733655513E-03, n_yes = 46, n_no = 0, h*lambda = 6.320713961201241E+00 step(A) = 155, err = 3.844142587215670E-01, h_new = 1.005451783931602E-02, n_yes = 47, n_no = 0, h*lambda = 6.327630943162317E+00 step(A) = 156, err = 3.729989106336811E-01, h_new = 1.023623653904493E-02, n_yes = 48, n_no = 0, h*lambda = 6.314981608071641E+00 step(A) = 157, err = 3.812381708884542E-01, h_new = 1.039281688370112E-02, n_yes = 49, n_no = 0, h*lambda = 6.322572443295686E+00 step(A) = 158, err = 3.683738871676462E-01, h_new = 1.059716460195157E-02, n_yes = 50, n_no = 0, h*lambda = 6.308222593956885E+00 step(A) = 159, err = 3.776236289379404E-01, h_new = 1.077208556342378E-02, n_yes = 51, n_no = 0, h*lambda = 6.316690366547824E+00 step(A) = 160, err = 3.628896196535183E-01, h_new = 1.100450432995160E-02, n_yes = 52, n_no = 0, h*lambda = 6.300133377250898E+00 step(A) = 161, err = 3.734804154794846E-01, h_new = 1.120158595867956E-02, n_yes = 53, n_no = 0, h*lambda = 6.309767543248209E+00 step(A) = 162, err = 3.562814773856994E-01, h_new = 1.146958937156355E-02, n_yes = 54, n_no = 0, h*lambda = 6.290278371489547E+00 step(A) = 163, err = 3.686965460528878E-01, h_new = 1.169382918217471E-02, n_yes = 55, n_no = 0, h*lambda = 6.301504456023957E+00 step(A) = 164, err = 3.481646290599441E-01, h_new = 1.200815192900144E-02, n_yes = 56, n_no = 0, h*lambda = 6.278010148914895E+00 step(A) = 165, err = 3.631360706292755E-01, h_new = 1.226619910083735E-02, n_yes = 57, n_no = 0, h*lambda = 6.291477240395897E+00 step(A) = 166, err = 3.379582856956321E-01, h_new = 1.264283971041080E-02, n_yes = 58, n_no = 0, h*lambda = 6.262324685769426E+00 step(A) = 167, err = 3.566446940446351E-01, h_new = 1.294367710420577E-02, n_yes = 59, n_no = 0, h*lambda = 6.279068319108486E+00 step(A) = 168, err = 3.247461179211537E-01, h_new = 1.340778952563008E-02, n_yes = 60, n_no = 0, h*lambda = 6.241581476243811E+00 step(A) = 169, err = 3.490835837448948E-01, h_new = 1.376364666425348E-02, n_yes = 61, n_no = 0, h*lambda = 6.263348396502632E+00 step(A) = 170, err = 3.070030098581453E-01, h_new = 1.435764448171737E-02, n_yes = 62, n_no = 0, h*lambda = 6.212918158244427E+00 step(A) = 171, err = 3.404624903158653E-01, h_new = 1.478485418962688E-02, n_yes = 63, n_no = 0, h*lambda = 6.242872417612695E+00 step(A) = 172, err = 2.820225026082013E-01, h_new = 1.558741452442216E-02, n_yes = 64, n_no = 0, h*lambda = 6.170898881415922E+00 step(A) = 173, err = 3.314517487726826E-01, h_new = 1.610512333618712E-02, n_yes = 65, n_no = 0, h*lambda = 6.215342444019660E+00 step(A) = 174, err = 2.446329688125062E-01, h_new = 1.728391776382303E-02, n_yes = 66, n_no = 0, h*lambda = 6.103967406267166E+00 step(A) = 175, err = 3.256218238131412E-01, h_new = 1.789762946662598E-02, n_yes = 67, n_no = 0, h*lambda = 6.177263852592048E+00 step(A) = 176, err = 1.842322066998068E-01, h_new = 1.990065303253532E-02, n_yes = 67, n_no = 1, h*lambda = 5.983408373653228E+00 step(A) = 177, err = 3.429173867205635E-01, h_new = 2.047439837110997E-02, n_yes = 68, n_no = 0, h*lambda = 6.126272396714477E+00 step(A) = 178, err = 7.954959625265076E-02, h_new = 2.528563770796823E-02, n_yes = 68, n_no = 1, h*lambda = 5.719907969586449E+00 step(A) = 179, err = 6.268881100096639E-01, h_new = 2.412501913194669E-02, n_yes = 69, n_no = 0, h*lambda = 6.125818987017078E+00 step(A) = 180, err = 5.092875291083272E-02, h_new = 3.150212102498735E-02, n_yes = 69, n_no = 1, h*lambda = 4.919484055217311E+00 step(A) = 181, err = 1.868841443849700E-01, h_new = 3.496517836518934E-02, n_yes = 69, n_no = 2, h*lambda = 4.669241473223804E+00 step(A) = 182, err = 3.326367005140728E-02, h_new = 4.815395770099203E-02, n_yes = 69, n_no = 3, h*lambda = 1.170798053973569E+00 step(R) = 183, err = 3.842551388784302E+02, h_new = 2.059667118831029E-02 step(A) = 184, err = 1.290014884507036E-01, h_new = 2.394502005451645E-02, n_yes = 69, n_no = 4, h*lambda = 1.037811969509999E+00 step(R) = 185, err = 1.637439086078181E+03, h_new = 7.349687963259307E-03 step(A) = 186, err = 3.387747535878958E-03, h_new = 1.346705180744629E-02, n_yes = 69, n_no = 5, h*lambda = 8.051118631407250E-01 step(A) = 187, err = 2.392761264869706E-02, h_new = 1.054748677259432E-02, n_yes = 0, n_no = 6, h*lambda = 2.017037316009412E+00 step(R) = 188, err = 5.883494026243961E+02, h_new = 4.277471102581100E-03 step(A) = 189, err = 1.063381786564274E-02, h_new = 6.793497837570425E-03, n_yes = 0, n_no = 7, h*lambda = 1.268282519357536E+00 step(R) = 190, err = 1.442587930668910E+02, h_new = 2.067934634023509E-03 step(A) = 191, err = 5.499728092426808E-02, h_new = 2.674463431230688E-03, n_yes = 0, n_no = 8, h*lambda = 2.355334203844016E-01 step(A) = 192, err = 1.744550738909289E-01, h_new = 2.315095259538681E-03, n_yes = 0, n_no = 9, h*lambda = 1.926758199851226E+00 step(A) = 193, err = 3.666623735391872E-01, h_new = 2.361990148756619E-03, n_yes = 0, n_no = 10, h*lambda = 2.322097313427576E+00 step(A) = 194, err = 4.677791748116633E-01, h_new = 2.337574607402342E-03, n_yes = 0, n_no = 11, h*lambda = 2.388609638868325E+00 step(A) = 195, err = 2.521798830936710E-01, h_new = 2.499160715372906E-03, n_yes = 0, n_no = 12, h*lambda = 2.362450678797520E+00 step(A) = 196, err = 3.993139473097741E-02, h_new = 3.364127018728178E-03, n_yes = 0, n_no = 13, h*lambda = 2.520556544738141E+00 step(A) = 197, err = 2.745227117436805E-02, h_new = 4.745616173047884E-03, n_yes = 0, n_no = 14, h*lambda = 3.382933879894489E+00 step(A) = 198, err = 4.695507521317095E-03, h_new = 8.347853162729110E-03, n_yes = 0, n_no = 15, h*lambda = 4.752191853774447E+00 step(A) = 199, err = 3.874437752120564E-02, h_new = 1.127953574234764E-02, n_yes = 1, n_no = 0, h*lambda = 8.297479848166985E+00 step(R) = 200, err = 1.788249104833765E+01, h_new = 7.079135614604501E-03 step(A) = 201, err = 8.780937804839464E-02, h_new = 8.635353664314799E-03, n_yes = 2, n_no = 0, h*lambda = 6.991810121159718E+00 step(A) = 202, err = 2.296072744807583E-01, h_new = 7.657720945088035E-03, n_yes = 3, n_no = 0, h*lambda = 6.947181426895245E+00 step(R) = 203, err = 1.540362789576154E+00, h_new = 6.529638953185796E-03 step(A) = 204, err = 1.623861307622340E-01, h_new = 7.375866727558742E-03, n_yes = 4, n_no = 0, h*lambda = 6.369816867435181E+00 step(A) = 205, err = 1.471538083519933E-01, h_new = 7.467242226372809E-03, n_yes = 5, n_no = 0, h*lambda = 6.331749647294362E+00 step(A) = 206, err = 9.091148569753943E-01, h_new = 6.801041179803955E-03, n_yes = 6, n_no = 0, h*lambda = 7.190873199670890E+00 step(A) = 207, err = 8.586960162117802E-01, h_new = 6.238612095569134E-03, n_yes = 7, n_no = 0, h*lambda = 6.508063535411213E+00 step(A) = 208, err = 1.522203722643146E-01, h_new = 7.104301674461251E-03, n_yes = 7, n_no = 1, h*lambda = 5.934729718342836E+00 step(A) = 209, err = 8.064588489309779E-01, h_new = 6.568120742556456E-03, n_yes = 8, n_no = 0, h*lambda = 6.713004288885192E+00 step(A) = 210, err = 3.428816463501300E-01, h_new = 6.757570844854359E-03, n_yes = 9, n_no = 0, h*lambda = 6.167297175289476E+00 step(A) = 211, err = 3.716952066558486E-01, h_new = 6.882714369078041E-03, n_yes = 10, n_no = 0, h*lambda = 6.304099123188388E+00 step(A) = 212, err = 4.108452573012789E-01, h_new = 6.922970348586769E-03, n_yes = 11, n_no = 0, h*lambda = 6.377862222260362E+00 step(A) = 213, err = 4.156417345362928E-01, h_new = 6.953365963313394E-03, n_yes = 12, n_no = 0, h*lambda = 6.371888280463802E+00 step(A) = 214, err = 4.032746931646385E-01, h_new = 7.010314075780550E-03, n_yes = 13, n_no = 0, h*lambda = 6.355835742997682E+00 step(A) = 215, err = 4.080227444621372E-01, h_new = 7.057395212359654E-03, n_yes = 14, n_no = 0, h*lambda = 6.363351711621393E+00 THE PROBLEM SEEMS TO BECOME STIFF AT X = 1.749270939102191E+00, ACCEPTED STEP = 196 step(A) = 216, err = 4.051403260117802E-01, h_new = 7.111091449931478E-03, n_yes = 15, n_no = 0, h*lambda = 6.360565720889913E+00 step(A) = 217, err = 4.083105985517257E-01, h_new = 7.158218347879135E-03, n_yes = 16, n_no = 0, h*lambda = 6.362957765869774E+00 step(A) = 218, err = 4.035563235730730E-01, h_new = 7.216214456782076E-03, n_yes = 17, n_no = 0, h*lambda = 6.358114388856899E+00 step(A) = 219, err = 4.075497194216475E-01, h_new = 7.265731855718158E-03, n_yes = 18, n_no = 0, h*lambda = 6.362067208141631E+00 step(A) = 220, err = 4.027779730018146E-01, h_new = 7.326366856221879E-03, n_yes = 19, n_no = 0, h*lambda = 6.357094001500717E+00 step(A) = 221, err = 4.067895963926024E-01, h_new = 7.378361699591563E-03, n_yes = 20, n_no = 0, h*lambda = 6.360918589785125E+00 step(A) = 222, err = 4.017445107270828E-01, h_new = 7.442326291922512E-03, n_yes = 21, n_no = 0, h*lambda = 6.355702205513821E+00 step(A) = 223, err = 4.059972145345996E-01, h_new = 7.496971061328855E-03, n_yes = 22, n_no = 0, h*lambda = 6.359782046983931E+00 step(A) = 224, err = 4.006831844194843E-01, h_new = 7.564464765291072E-03, n_yes = 23, n_no = 0, h*lambda = 6.354257829357048E+00 step(A) = 225, err = 4.050928333183331E-01, h_new = 7.622130738910705E-03, n_yes = 24, n_no = 0, h*lambda = 6.358516102599155E+00 step(A) = 226, err = 3.996066820502495E-01, h_new = 7.693337951817826E-03, n_yes = 25, n_no = 0, h*lambda = 6.352760684694431E+00 step(A) = 227, err = 4.040853787497352E-01, h_new = 7.754399609695395E-03, n_yes = 26, n_no = 0, h*lambda = 6.357079519999597E+00 step(A) = 228, err = 3.984334159525326E-01, h_new = 7.829719756616404E-03, n_yes = 27, n_no = 0, h*lambda = 6.351136028902602E+00 step(A) = 229, err = 4.030094038080632E-01, h_new = 7.894494561424074E-03, n_yes = 28, n_no = 0, h*lambda = 6.355540435749852E+00 step(A) = 230, err = 3.971588512305952E-01, h_new = 7.974368643476979E-03, n_yes = 29, n_no = 0, h*lambda = 6.349364394057486E+00 step(A) = 231, err = 4.018455938617868E-01, h_new = 8.043247204377182E-03, n_yes = 30, n_no = 0, h*lambda = 6.353867465375588E+00 step(A) = 232, err = 3.957694751760499E-01, h_new = 8.128186114056779E-03, n_yes = 31, n_no = 0, h*lambda = 6.347428901129612E+00 step(A) = 233, err = 4.005852761546830E-01, h_new = 8.201613056816588E-03, n_yes = 32, n_no = 0, h*lambda = 6.352045756675778E+00 step(A) = 234, err = 3.942483877368551E-01, h_new = 8.292214821660902E-03, n_yes = 33, n_no = 0, h*lambda = 6.345304039670358E+00 step(A) = 235, err = 6.469446879417725E-05, h_new = 6.908420682852039E-03, n_yes = 33, n_no = 1, h*lambda = 1.773473114548921E+00 DoPri8: Dormand-Prince method (explicit, order 8(5,3), embedded) Number of function evaluations = 2802 Number of performed steps = 235 Number of accepted steps = 215 Number of rejected steps = 20 y = 1.819907445729370E+00 -7.866363461162956E-01 h = 6.908420682852039E-03