running dopri5.f test step(A) = 1, err = 3.727760587956869E-06, h_new = 5.213103498972100E-04, n_yes = 0, n_no = 1, h*lambda = 9.997470703048382E-02 step(A) = 2, err = 1.240953492473138E-02, h_new = 6.845425066468906E-04, n_yes = 0, n_no = 2, h*lambda = 5.212932725594240E-01 step(A) = 3, err = 2.688194416514796E-02, h_new = 9.558391963414619E-04, n_yes = 0, n_no = 3, h*lambda = 6.844363978402807E-01 step(A) = 4, err = 7.358690218532347E-02, h_new = 1.159977816134844E-03, n_yes = 0, n_no = 4, h*lambda = 9.551832656577894E-01 step(A) = 5, err = 7.734184510533144E-02, h_new = 1.453227956747536E-03, n_yes = 0, n_no = 5, h*lambda = 1.158113848815346E+00 step(A) = 6, err = 8.155247362351542E-02, h_new = 1.807875507561472E-03, n_yes = 0, n_no = 6, h*lambda = 1.449058825392334E+00 step(A) = 7, err = 6.439628755218164E-02, h_new = 2.346184573138853E-03, n_yes = 0, n_no = 7, h*lambda = 1.799776360945848E+00 step(A) = 8, err = 5.074827815250559E-02, h_new = 3.140779839550399E-03, n_yes = 0, n_no = 8, h*lambda = 2.330763902743438E+00 step(A) = 9, err = 5.307555787744915E-02, h_new = 4.132994048947248E-03, n_yes = 0, n_no = 9, h*lambda = 3.111350910805671E+00 step(A) = 10, err = 1.828230653840881E-01, h_new = 4.415275469091488E-03, n_yes = 1, n_no = 0, h*lambda = 4.079076160810598E+00 step(A) = 11, err = 9.983777781302734E-01, h_new = 3.713653838225876E-03, n_yes = 2, n_no = 0, h*lambda = 4.340439785241045E+00 step(R) = 12, err = 1.958559439663303E+00, h_new = 2.981361783324316E-03 step(A) = 13, err = 5.429798041880982E-01, h_new = 2.976567673113339E-03, n_yes = 2, n_no = 1, h*lambda = 2.922935427793524E+00 step(A) = 14, err = 2.604231783693373E-01, h_new = 3.286124999538735E-03, n_yes = 2, n_no = 2, h*lambda = 2.910250138596906E+00 step(A) = 15, err = 2.185708904414213E-01, h_new = 3.629301761792994E-03, n_yes = 2, n_no = 3, h*lambda = 3.203240293388239E+00 step(A) = 16, err = 3.190115285501510E-01, h_new = 3.732520464824205E-03, n_yes = 3, n_no = 0, h*lambda = 3.525998778095849E+00 step(A) = 17, err = 5.522585322536945E-01, h_new = 3.550035927656333E-03, n_yes = 4, n_no = 0, h*lambda = 3.613863991483097E+00 step(A) = 18, err = 7.044683881386191E-01, h_new = 3.311498010544673E-03, n_yes = 5, n_no = 0, h*lambda = 3.425947123552543E+00 step(A) = 19, err = 5.762683559805580E-01, h_new = 3.227566930142037E-03, n_yes = 5, n_no = 1, h*lambda = 3.185913475885057E+00 step(A) = 20, err = 3.923779967127298E-01, h_new = 3.331294038260041E-03, n_yes = 5, n_no = 2, h*lambda = 3.095782642831596E+00 step(A) = 21, err = 3.150183898294197E-01, h_new = 3.514686219593283E-03, n_yes = 5, n_no = 3, h*lambda = 3.185287472114879E+00 step(A) = 22, err = 3.399300730445723E-01, h_new = 3.628492873521982E-03, n_yes = 6, n_no = 0, h*lambda = 3.349542976403127E+00 step(A) = 23, err = 4.377868682190470E-01, h_new = 3.599230854914461E-03, n_yes = 7, n_no = 0, h*lambda = 3.446184457471084E+00 step(A) = 24, err = 5.314409248835892E-01, h_new = 3.489597785948117E-03, n_yes = 8, n_no = 0, h*lambda = 3.406762217146414E+00 step(A) = 25, err = 5.258363701584914E-01, h_new = 3.415792466049110E-03, n_yes = 9, n_no = 0, h*lambda = 3.292033490403043E+00 step(A) = 26, err = 4.464080548569233E-01, h_new = 3.436477991049893E-03, n_yes = 9, n_no = 1, h*lambda = 3.211875898655271E+00 step(A) = 27, err = 3.833585347442178E-01, h_new = 3.524783801346556E-03, n_yes = 9, n_no = 2, h*lambda = 3.220655297222187E+00 step(A) = 28, err = 3.745240530606989E-01, h_new = 3.607676997473170E-03, n_yes = 10, n_no = 0, h*lambda = 3.292190171398073E+00 step(A) = 29, err = 4.125127147932751E-01, h_new = 3.628983534443945E-03, n_yes = 11, n_no = 0, h*lambda = 3.357856986287223E+00 step(A) = 30, err = 4.629698432211843E-01, h_new = 3.593362582513521E-03, n_yes = 12, n_no = 0, h*lambda = 3.365788505507941E+00 step(A) = 31, err = 4.807404728719064E-01, h_new = 3.551737348006160E-03, n_yes = 13, n_no = 0, h*lambda = 3.321068963245464E+00 step(A) = 32, err = 4.552668253450887E-01, h_new = 3.548579653715641E-03, n_yes = 14, n_no = 0, h*lambda = 3.271163656930549E+00 THE PROBLEM SEEMS TO BECOME STIFF AT X = 1.216973774601867E-02, ACCEPTED STEP = 32 step(A) = 33, err = 4.187997643679062E-01, h_new = 3.588282399522332E-03, n_yes = 15, n_no = 0, h*lambda = 3.256824769569957E+00 step(A) = 34, err = 4.023721900725468E-01, h_new = 3.641016347314295E-03, n_yes = 16, n_no = 0, h*lambda = 3.281567253670414E+00 step(A) = 35, err = 4.128776652616866E-01, h_new = 3.672489985558470E-03, n_yes = 17, n_no = 0, h*lambda = 3.317746879037275E+00 step(A) = 36, err = 4.372081522525559E-01, h_new = 3.672137905119492E-03, n_yes = 18, n_no = 0, h*lambda = 3.334163600768078E+00 step(A) = 37, err = 4.534015784165334E-01, h_new = 3.657521664994581E-03, n_yes = 19, n_no = 0, h*lambda = 3.321571853459977E+00 step(A) = 38, err = 4.491089007339857E-01, h_new = 3.654171745492282E-03, n_yes = 20, n_no = 0, h*lambda = 3.296159844234577E+00 step(A) = 39, err = 4.322176200204242E-01, h_new = 3.673297571645196E-03, n_yes = 21, n_no = 0, h*lambda = 3.280955619906580E+00 step(A) = 40, err = 4.191700439023628E-01, h_new = 3.706127699014764E-03, n_yes = 22, n_no = 0, h*lambda = 3.285802075575150E+00 step(A) = 41, err = 4.190025846064984E-01, h_new = 3.734923053989653E-03, n_yes = 23, n_no = 0, h*lambda = 3.302611873943486E+00 step(A) = 42, err = 4.290236789614592E-01, h_new = 3.748789215377386E-03, n_yes = 24, n_no = 0, h*lambda = 3.315509301368293E+00 step(A) = 43, err = 4.392810014279703E-01, h_new = 3.751168478178395E-03, n_yes = 25, n_no = 0, h*lambda = 3.314948466779248E+00 step(A) = 44, err = 4.412782483727336E-01, h_new = 3.754202097903182E-03, n_yes = 26, n_no = 0, h*lambda = 3.304151076523684E+00 step(A) = 45, err = 4.349692823545975E-01, h_new = 3.767130787642624E-03, n_yes = 27, n_no = 0, h*lambda = 3.293884821666195E+00 step(A) = 46, err = 4.271791251749550E-01, h_new = 3.789551797221630E-03, n_yes = 28, n_no = 0, h*lambda = 3.292185526777186E+00 step(A) = 47, err = 4.242905792585751E-01, h_new = 3.813747902134294E-03, n_yes = 29, n_no = 0, h*lambda = 3.298567850297265E+00 step(A) = 48, err = 4.274049868596863E-01, h_new = 3.832289394394829E-03, n_yes = 30, n_no = 0, h*lambda = 3.306234916786778E+00 step(A) = 49, err = 4.327411889339119E-01, h_new = 3.843931062320482E-03, n_yes = 31, n_no = 0, h*lambda = 3.308770400607540E+00 step(A) = 50, err = 4.354784875714566E-01, h_new = 3.853389318092278E-03, n_yes = 32, n_no = 0, h*lambda = 3.305185144258213E+00 step(A) = 51, err = 4.338299364634911E-01, h_new = 3.866337385609927E-03, n_yes = 33, n_no = 0, h*lambda = 3.299597220826963E+00 step(A) = 52, err = 4.299482582564692E-01, h_new = 3.884671380738773E-03, n_yes = 34, n_no = 0, h*lambda = 3.296854674172394E+00 step(A) = 53, err = 4.273038470191889E-01, h_new = 3.905783683190740E-03, n_yes = 35, n_no = 0, h*lambda = 3.298510698645459E+00 step(A) = 54, err = 4.276130562871120E-01, h_new = 3.925558968985641E-03, n_yes = 36, n_no = 0, h*lambda = 3.302291579858245E+00 step(A) = 55, err = 4.299022185899039E-01, h_new = 3.941969018536166E-03, n_yes = 37, n_no = 0, h*lambda = 3.304705624712688E+00 step(A) = 56, err = 4.318278667904522E-01, h_new = 3.956286103097857E-03, n_yes = 38, n_no = 0, h*lambda = 3.304077414865707E+00 step(A) = 57, err = 4.318055393702508E-01, h_new = 3.971399995522599E-03, n_yes = 39, n_no = 0, h*lambda = 3.301511376599565E+00 step(A) = 58, err = 4.301433302689515E-01, h_new = 3.989178098074373E-03, n_yes = 40, n_no = 0, h*lambda = 3.299427218684822E+00 step(A) = 59, err = 4.283980039980445E-01, h_new = 4.009187785386613E-03, n_yes = 41, n_no = 0, h*lambda = 3.299349782509169E+00 step(A) = 60, err = 4.278528298765437E-01, h_new = 4.029514807416680E-03, n_yes = 42, n_no = 0, h*lambda = 3.300883653773780E+00 step(A) = 61, err = 4.285629478102378E-01, h_new = 4.048597068577122E-03, n_yes = 43, n_no = 0, h*lambda = 3.302431999138754E+00 step(A) = 62, err = 4.295742454422359E-01, h_new = 4.066409866812659E-03, n_yes = 44, n_no = 0, h*lambda = 3.302719490357073E+00 step(A) = 63, err = 4.298863519095621E-01, h_new = 4.084181818148554E-03, n_yes = 45, n_no = 0, h*lambda = 3.301742927630345E+00 step(A) = 64, err = 4.292626614581955E-01, h_new = 4.103163227947270E-03, n_yes = 46, n_no = 0, h*lambda = 3.300508019483566E+00 step(A) = 65, err = 4.282627026612512E-01, h_new = 4.123628047632168E-03, n_yes = 47, n_no = 0, h*lambda = 3.300014559491820E+00 step(A) = 66, err = 4.276241233630349E-01, h_new = 4.144859665351713E-03, n_yes = 48, n_no = 0, h*lambda = 3.300460431134471E+00 step(A) = 67, err = 4.276369328994878E-01, h_new = 4.165930719970470E-03, n_yes = 49, n_no = 0, h*lambda = 3.301252668087124E+00 step(A) = 68, err = 4.280092024193676E-01, h_new = 4.186494575649766E-03, n_yes = 50, n_no = 0, h*lambda = 3.301645687404717E+00 step(A) = 69, err = 4.282204831855134E-01, h_new = 4.206953407913545E-03, n_yes = 51, n_no = 0, h*lambda = 3.301367652405347E+00 step(A) = 70, err = 4.279837977259472E-01, h_new = 4.227993036622881E-03, n_yes = 52, n_no = 0, h*lambda = 3.300738226216595E+00 step(A) = 71, err = 4.274266185959482E-01, h_new = 4.249985025858271E-03, n_yes = 53, n_no = 0, h*lambda = 3.300289463936793E+00 step(A) = 72, err = 4.269015487143270E-01, h_new = 4.272761559291803E-03, n_yes = 54, n_no = 0, h*lambda = 3.300296852305110E+00 step(A) = 73, err = 4.266570641942520E-01, h_new = 4.295867290392442E-03, n_yes = 55, n_no = 0, h*lambda = 3.300613643146745E+00 step(A) = 74, err = 4.266619080992640E-01, h_new = 4.318990665845511E-03, n_yes = 56, n_no = 0, h*lambda = 3.300876408473815E+00 step(A) = 75, err = 4.266887927143362E-01, h_new = 4.342193967278273E-03, n_yes = 57, n_no = 0, h*lambda = 3.300840114993769E+00 step(A) = 76, err = 4.265362055841218E-01, h_new = 4.365798379295038E-03, n_yes = 58, n_no = 0, h*lambda = 3.300548680184669E+00 step(A) = 77, err = 4.261836256260697E-01, h_new = 4.390085430569889E-03, n_yes = 59, n_no = 0, h*lambda = 3.300239345765668E+00 step(A) = 78, err = 4.257697261779472E-01, h_new = 4.415090795345218E-03, n_yes = 60, n_no = 0, h*lambda = 3.300112913642361E+00 step(A) = 79, err = 4.254479431183302E-01, h_new = 4.440636731222893E-03, n_yes = 61, n_no = 0, h*lambda = 3.300181318849538E+00 step(A) = 80, err = 4.252609652723980E-01, h_new = 4.466529173762782E-03, n_yes = 62, n_no = 0, h*lambda = 3.300297748188529E+00 step(A) = 81, err = 4.251299620507933E-01, h_new = 4.492728906890641E-03, n_yes = 63, n_no = 0, h*lambda = 3.300306923355136E+00 step(A) = 82, err = 4.249426797540556E-01, h_new = 4.519365146404788E-03, n_yes = 64, n_no = 0, h*lambda = 3.300170619890225E+00 step(A) = 83, err = 4.246496712810864E-01, h_new = 4.546612282543069E-03, n_yes = 65, n_no = 0, h*lambda = 3.299973567617452E+00 step(A) = 84, err = 4.242902785325018E-01, h_new = 4.574555891799050E-03, n_yes = 66, n_no = 0, h*lambda = 3.299830092028819E+00 step(A) = 85, err = 4.239422960826399E-01, h_new = 4.603157382739926E-03, n_yes = 67, n_no = 0, h*lambda = 3.299786765752348E+00 step(A) = 86, err = 4.236516104084949E-01, h_new = 4.632325801040992E-03, n_yes = 68, n_no = 0, h*lambda = 3.299799497146005E+00 step(A) = 87, err = 4.234021294822136E-01, h_new = 4.662017979379640E-03, n_yes = 69, n_no = 0, h*lambda = 3.299787665247680E+00 step(A) = 88, err = 4.231413596539532E-01, h_new = 4.692281342521445E-03, n_yes = 70, n_no = 0, h*lambda = 3.299706001124871E+00 step(A) = 89, err = 4.228295479462742E-01, h_new = 4.723216647528324E-03, n_yes = 71, n_no = 0, h*lambda = 3.299573225257547E+00 step(A) = 90, err = 4.224683252333351E-01, h_new = 4.754906518256444E-03, n_yes = 72, n_no = 0, h*lambda = 3.299443688642741E+00 step(A) = 91, err = 4.220903214023303E-01, h_new = 4.787373833360153E-03, n_yes = 73, n_no = 0, h*lambda = 3.299356575427231E+00 step(A) = 92, err = 4.217262527486801E-01, h_new = 4.820597360350105E-03, n_yes = 74, n_no = 0, h*lambda = 3.299307852672613E+00 step(A) = 93, err = 4.213804668437543E-01, h_new = 4.854560809514677E-03, n_yes = 75, n_no = 0, h*lambda = 3.299262534479612E+00 step(A) = 94, err = 4.210326121934276E-01, h_new = 4.889289530465294E-03, n_yes = 76, n_no = 0, h*lambda = 3.299189029321810E+00 step(A) = 95, err = 4.206585781340944E-01, h_new = 4.924848070603226E-03, n_yes = 77, n_no = 0, h*lambda = 3.299083582428204E+00 step(A) = 96, err = 4.202499803473896E-01, h_new = 4.961308437665509E-03, n_yes = 78, n_no = 0, h*lambda = 3.298967423756141E+00 step(A) = 97, err = 4.198169838146367E-01, h_new = 4.998720385112455E-03, n_yes = 79, n_no = 0, h*lambda = 3.298864619069140E+00 step(A) = 98, err = 4.193757250581133E-01, h_new = 5.037107212991481E-03, n_yes = 80, n_no = 0, h*lambda = 3.298782313284495E+00 step(A) = 99, err = 4.189335980660140E-01, h_new = 5.076485536099562E-03, n_yes = 81, n_no = 0, h*lambda = 3.298708596717074E+00 step(A) = 100, err = 4.184845596148351E-01, h_new = 5.116888642345174E-03, n_yes = 82, n_no = 0, h*lambda = 3.298625996384025E+00 step(A) = 101, err = 4.180160636392798E-01, h_new = 5.158374245907839E-03, n_yes = 83, n_no = 0, h*lambda = 3.298526576377812E+00 step(A) = 102, err = 4.175196485347982E-01, h_new = 5.201013723194422E-03, n_yes = 84, n_no = 0, h*lambda = 3.298416291142185E+00 step(A) = 103, err = 4.169959803376240E-01, h_new = 5.244875313946208E-03, n_yes = 85, n_no = 0, h*lambda = 3.298307381257485E+00 step(A) = 104, err = 4.164515043345261E-01, h_new = 5.290016156001509E-03, n_yes = 86, n_no = 0, h*lambda = 3.298207372746440E+00 step(A) = 105, err = 4.158911347372369E-01, h_new = 5.336488066000710E-03, n_yes = 87, n_no = 0, h*lambda = 3.298114148690524E+00 step(A) = 106, err = 4.153136252589897E-01, h_new = 5.384350066409061E-03, n_yes = 88, n_no = 0, h*lambda = 3.298019744859849E+00 step(A) = 107, err = 4.147127240176202E-01, h_new = 5.433676684227517E-03, n_yes = 89, n_no = 0, h*lambda = 3.297918025572469E+00 step(A) = 108, err = 4.140819615578625E-01, h_new = 5.484556618667384E-03, n_yes = 90, n_no = 0, h*lambda = 3.297809325718488E+00 step(A) = 109, err = 4.134184881002320E-01, h_new = 5.537085173718872E-03, n_yes = 91, n_no = 0, h*lambda = 3.297699007140089E+00 step(A) = 110, err = 4.127230989653087E-01, h_new = 5.591358225306473E-03, n_yes = 92, n_no = 0, h*lambda = 3.297592357070076E+00 step(A) = 111, err = 4.119972187572797E-01, h_new = 5.647472821220519E-03, n_yes = 93, n_no = 0, h*lambda = 3.297490657401036E+00 step(A) = 112, err = 4.112398654530366E-01, h_new = 5.705533306165729E-03, n_yes = 94, n_no = 0, h*lambda = 3.297391313375605E+00 step(A) = 113, err = 4.104469888410856E-01, h_new = 5.765657764230638E-03, n_yes = 95, n_no = 0, h*lambda = 3.297291101569820E+00 step(A) = 114, err = 4.096131826626377E-01, h_new = 5.827980436384252E-03, n_yes = 96, n_no = 0, h*lambda = 3.297189317787390E+00 step(A) = 115, err = 4.087338827100014E-01, h_new = 5.892649949218874E-03, n_yes = 97, n_no = 0, h*lambda = 3.297088360977995E+00 step(A) = 116, err = 4.078061727002866E-01, h_new = 5.959826715836608E-03, n_yes = 98, n_no = 0, h*lambda = 3.296991843789327E+00 step(A) = 117, err = 4.068278160431442E-01, h_new = 6.029683065901680E-03, n_yes = 99, n_no = 0, h*lambda = 3.296902269471037E+00 step(A) = 118, err = 4.057956124292814E-01, h_new = 6.102407023815352E-03, n_yes = 100, n_no = 0, h*lambda = 3.296820238846854E+00 step(A) = 119, err = 4.047044721772358E-01, h_new = 6.178207835895364E-03, n_yes = 101, n_no = 0, h*lambda = 3.296745534749260E+00 step(A) = 120, err = 4.035477273735152E-01, h_new = 6.257320635686081E-03, n_yes = 102, n_no = 0, h*lambda = 3.296678879229823E+00 step(A) = 121, err = 4.023180971205103E-01, h_new = 6.340009206905153E-03, n_yes = 103, n_no = 0, h*lambda = 3.296622888215403E+00 step(A) = 122, err = 4.010083145404920E-01, h_new = 6.426567990759728E-03, n_yes = 104, n_no = 0, h*lambda = 3.296581706300584E+00 step(A) = 123, err = 3.996108814224650E-01, h_new = 6.517325459696554E-03, n_yes = 105, n_no = 0, h*lambda = 3.296559976829165E+00 step(A) = 124, err = 3.981171943075251E-01, h_new = 6.612650233252473E-03, n_yes = 106, n_no = 0, h*lambda = 3.296562218466875E+00 step(A) = 125, err = 3.965166987547910E-01, h_new = 6.712959806131000E-03, n_yes = 107, n_no = 0, h*lambda = 3.296593170406035E+00 step(A) = 126, err = 3.947965257597619E-01, h_new = 6.818730889981018E-03, n_yes = 108, n_no = 0, h*lambda = 3.296658832778334E+00 step(A) = 127, err = 3.929415338950972E-01, h_new = 6.930510789484362E-03, n_yes = 109, n_no = 0, h*lambda = 3.296767499179250E+00 step(A) = 128, err = 3.909343060352350E-01, h_new = 7.048930485777037E-03, n_yes = 110, n_no = 0, h*lambda = 3.296930315357218E+00 step(A) = 129, err = 3.887546851856738E-01, h_new = 7.174721211828735E-03, n_yes = 111, n_no = 0, h*lambda = 3.297161486201307E+00 step(A) = 130, err = 3.863787537776220E-01, h_new = 7.308736671995058E-03, n_yes = 112, n_no = 0, h*lambda = 3.297478671472719E+00 step(A) = 131, err = 3.837774367222587E-01, h_new = 7.451982895310504E-03, n_yes = 113, n_no = 0, h*lambda = 3.297904088915819E+00 step(A) = 132, err = 3.809149091891968E-01, h_new = 7.605657796493046E-03, n_yes = 114, n_no = 0, h*lambda = 3.298466526325542E+00 step(A) = 133, err = 3.777467356030886E-01, h_new = 7.771203595909962E-03, n_yes = 115, n_no = 0, h*lambda = 3.299204218927224E+00 step(A) = 134, err = 3.742173477307187E-01, h_new = 7.950377668671118E-03, n_yes = 116, n_no = 0, h*lambda = 3.300168645487699E+00 step(A) = 135, err = 3.702562431552448E-01, h_new = 8.145351292688528E-03, n_yes = 117, n_no = 0, h*lambda = 3.301429746766151E+00 step(A) = 136, err = 3.657721211404258E-01, h_new = 8.358851738839223E-03, n_yes = 118, n_no = 0, h*lambda = 3.303083723809014E+00 step(A) = 137, err = 3.606438455258906E-01, h_new = 8.594373159045208E-03, n_yes = 119, n_no = 0, h*lambda = 3.305265376009520E+00 step(A) = 138, err = 3.547062731112082E-01, h_new = 8.856500450814979E-03, n_yes = 120, n_no = 0, h*lambda = 3.308168166822950E+00 step(A) = 139, err = 3.477270648763295E-01, h_new = 9.151428041711622E-03, n_yes = 121, n_no = 0, h*lambda = 3.312077654455514E+00 step(A) = 140, err = 3.393664556519749E-01, h_new = 9.487836744273925E-03, n_yes = 122, n_no = 0, h*lambda = 3.317429288356214E+00 step(A) = 141, err = 3.291024947836493E-01, h_new = 9.878481014840141E-03, n_yes = 123, n_no = 0, h*lambda = 3.324913729029330E+00 step(A) = 142, err = 3.160801716838454E-01, h_new = 1.034333028832478E-02, n_yes = 124, n_no = 0, h*lambda = 3.335682132106093E+00 step(A) = 143, err = 2.987727607038674E-01, h_new = 1.091658554969621E-02, n_yes = 125, n_no = 0, h*lambda = 3.351781701060984E+00 step(A) = 144, err = 2.741028248923998E-01, h_new = 1.166534735209562E-02, n_yes = 126, n_no = 0, h*lambda = 3.377188017194614E+00 step(A) = 145, err = 2.346239529623908E-01, h_new = 1.275538071971833E-02, n_yes = 127, n_no = 0, h*lambda = 3.420643301256491E+00 step(A) = 146, err = 1.557890191646284E-01, h_new = 1.486002588345571E-02, n_yes = 128, n_no = 0, h*lambda = 3.504653543083134E+00 step(A) = 147, err = 1.226099301421332E-01, h_new = 1.773838960147024E-02, n_yes = 129, n_no = 0, h*lambda = 3.804979086279828E+00 step(R) = 148, err = 1.865508033799077E+00, h_new = 1.435890055765922E-02 step(A) = 149, err = 5.567258694180440E-01, h_new = 1.312640152132764E-02, n_yes = 130, n_no = 0, h*lambda = 3.336672497878781E+00 step(A) = 150, err = 4.497492134490652E-01, h_new = 1.321931282903480E-02, n_yes = 130, n_no = 1, h*lambda = 2.776090342501430E+00 step(A) = 151, err = 2.958470302177168E-01, h_new = 1.417388091543196E-02, n_yes = 130, n_no = 2, h*lambda = 2.510927312771701E+00 step(A) = 152, err = 2.808093537513418E-01, h_new = 1.507800943483606E-02, n_yes = 130, n_no = 3, h*lambda = 2.346913131343200E+00 step(A) = 153, err = 2.677710048675263E-01, h_new = 1.613627102742321E-02, n_yes = 130, n_no = 4, h*lambda = 2.080888195177723E+00 step(A) = 154, err = 2.265592527176492E-01, h_new = 1.773271868994997E-02, n_yes = 130, n_no = 5, h*lambda = 1.724597526689825E+00 step(A) = 155, err = 1.595875690407286E-01, h_new = 2.054543146096174E-02, n_yes = 0, n_no = 6, h*lambda = 1.234987459122405E+00 step(R) = 156, err = 2.116701634616780E+00, h_new = 1.627779728778869E-02 step(A) = 157, err = 3.609774241177324E-01, h_new = 1.618775044370163E-02, n_yes = 0, n_no = 7, h*lambda = 9.990658038115804E-02 step(R) = 158, err = 9.287356029822718E+00, h_new = 9.974413167062967E-03 step(R) = 159, err = 1.287900027183122E+00, h_new = 8.599038580060907E-03 step(A) = 160, err = 6.026470923287880E-01, h_new = 8.098058083031089E-03, n_yes = 0, n_no = 8, h*lambda = 3.967658571562481E-01 step(A) = 161, err = 8.494145900007243E-01, h_new = 7.343034483198105E-03, n_yes = 0, n_no = 9, h*lambda = 1.828023413537553E+00 step(R) = 162, err = 2.666288518306080E+07, h_new = 1.468606896639621E-03 step(A) = 163, err = 3.230695729117591E-02, h_new = 2.353599329454549E-03, n_yes = 0, n_no = 10, h*lambda = 2.238398995440441E+00 step(A) = 164, err = 3.163001418030857E-01, h_new = 1.401209326574877E-03, n_yes = 0, n_no = 11, h*lambda = 1.436600222723829E+00 step(A) = 165, err = 6.786979900907592E-01, h_new = 1.286365809968859E-03, n_yes = 0, n_no = 12, h*lambda = 8.443327662465148E-01 step(A) = 166, err = 8.209750533747804E-01, h_new = 1.178794137219260E-03, n_yes = 0, n_no = 13, h*lambda = 1.192678783949197E+00 step(R) = 167, err = 1.576228468234167E+00, h_new = 9.819406027265132E-04 step(A) = 168, err = 5.647098916180083E-01, h_new = 9.662524518772293E-04, n_yes = 0, n_no = 14, h*lambda = 9.651650874607449E-01 step(A) = 169, err = 6.412908858800414E-01, h_new = 9.166568172831216E-04, n_yes = 0, n_no = 15, h*lambda = 9.687036860166818E-01 step(A) = 170, err = 5.052997920524998E-01, h_new = 9.101817736226046E-04, n_yes = 0, n_no = 16, h*lambda = 9.239742915135311E-01 step(A) = 171, err = 4.628627891109816E-01, h_new = 9.086271786241421E-04, n_yes = 0, n_no = 17, h*lambda = 9.192354929623804E-01 step(A) = 172, err = 3.911348116724248E-01, h_new = 9.301453586853670E-04, n_yes = 0, n_no = 18, h*lambda = 9.179383411554644E-01 step(A) = 173, err = 1.552925699886899E-01, h_new = 1.106600281751699E-03, n_yes = 0, n_no = 19, h*lambda = 9.393384125637803E-01 step(A) = 174, err = 1.436442279292664E-01, h_new = 1.285702008912774E-03, n_yes = 0, n_no = 20, h*lambda = 1.116711101707389E+00 step(A) = 175, err = 1.031704203432106E-01, h_new = 1.575324310472909E-03, n_yes = 0, n_no = 21, h*lambda = 1.296054764415302E+00 step(A) = 176, err = 8.526097529883833E-02, h_new = 1.967556182218248E-03, n_yes = 0, n_no = 22, h*lambda = 1.585819982675070E+00 step(A) = 177, err = 6.180674762187596E-02, h_new = 2.575867751683374E-03, n_yes = 0, n_no = 23, h*lambda = 1.977216120725463E+00 step(A) = 178, err = 4.869859823172219E-02, h_new = 3.466807351524985E-03, n_yes = 0, n_no = 24, h*lambda = 2.582612950428378E+00 step(A) = 179, err = 7.299927040061541E-02, h_new = 4.314283980065618E-03, n_yes = 1, n_no = 0, h*lambda = 3.465213353735224E+00 step(A) = 180, err = 3.505180381762840E-01, h_new = 4.179104349195821E-03, n_yes = 2, n_no = 0, h*lambda = 4.295805939984726E+00 step(R) = 181, err = 1.485964613744471E+00, h_new = 3.516286758978471E-03 step(A) = 182, err = 5.349973884011944E-01, h_new = 3.375168250427538E-03, n_yes = 3, n_no = 0, h*lambda = 3.490303241682396E+00 step(A) = 183, err = 5.803999025955329E-01, h_new = 3.249663232632390E-03, n_yes = 4, n_no = 0, h*lambda = 3.340113367891196E+00 step(A) = 184, err = 4.892905749726796E-01, h_new = 3.231495793130896E-03, n_yes = 4, n_no = 1, h*lambda = 3.206498018435824E+00 step(A) = 185, err = 3.892551199642380E-01, h_new = 3.318095958868127E-03, n_yes = 4, n_no = 2, h*lambda = 3.179243724110409E+00 step(A) = 186, err = 3.547588633304440E-01, h_new = 3.429668708741459E-03, n_yes = 5, n_no = 0, h*lambda = 3.254612180458219E+00 step(A) = 187, err = 3.866476014858118E-01, h_new = 3.480554241335775E-03, n_yes = 6, n_no = 0, h*lambda = 3.353554756146918E+00 step(A) = 188, err = 4.540719782326855E-01, h_new = 3.448834301218070E-03, n_yes = 7, n_no = 0, h*lambda = 3.392504829840897E+00 step(A) = 189, err = 4.976968146050750E-01, h_new = 3.386224911348940E-03, n_yes = 8, n_no = 0, h*lambda = 3.350968749649525E+00 step(A) = 190, err = 4.794714119559821E-01, h_new = 3.358205296088036E-03, n_yes = 9, n_no = 0, h*lambda = 3.279879410233608E+00 step(A) = 191, err = 4.298737629957955E-01, h_new = 3.387757385223475E-03, n_yes = 9, n_no = 1, h*lambda = 3.242633574384862E+00 step(A) = 192, err = 3.974072197874833E-01, h_new = 3.448405822707422E-03, n_yes = 10, n_no = 0, h*lambda = 3.260875897281569E+00 step(A) = 193, err = 4.004290084585745E-01, h_new = 3.494628133623199E-03, n_yes = 11, n_no = 0, h*lambda = 3.308587187955395E+00 step(A) = 194, err = 4.292829756273832E-01, h_new = 3.500886866145740E-03, n_yes = 12, n_no = 0, h*lambda = 3.341980296964947E+00 step(A) = 195, err = 4.573365638530497E-01, h_new = 3.479286767411650E-03, n_yes = 13, n_no = 0, h*lambda = 3.336964452469972E+00 step(A) = 196, err = 4.607508183960181E-01, h_new = 3.462206219161974E-03, n_yes = 14, n_no = 0, h*lambda = 3.305496641646690E+00 THE PROBLEM SEEMS TO BECOME STIFF AT X = 8.717646581250652E-01, ACCEPTED STEP = 189 step(A) = 197, err = 4.414110580345363E-01, h_new = 3.471448611779810E-03, n_yes = 15, n_no = 0, h*lambda = 3.278481491393737E+00 step(A) = 198, err = 4.205546755873187E-01, h_new = 3.503459967418862E-03, n_yes = 16, n_no = 0, h*lambda = 3.276376434522063E+00 step(A) = 199, err = 4.145522690635501E-01, h_new = 3.537562218866223E-03, n_yes = 17, n_no = 0, h*lambda = 3.295523772606389E+00 step(A) = 200, err = 4.245226729681265E-01, h_new = 3.555548581171666E-03, n_yes = 18, n_no = 0, h*lambda = 3.316314660236634E+00 step(A) = 201, err = 4.396009549220398E-01, h_new = 3.555864434972601E-03, n_yes = 19, n_no = 0, h*lambda = 3.321766052386071E+00 step(A) = 202, err = 4.465523480942285E-01, h_new = 3.551662957641965E-03, n_yes = 20, n_no = 0, h*lambda = 3.310637982341947E+00 step(A) = 203, err = 4.411207989226998E-01, h_new = 3.557086033022827E-03, n_yes = 21, n_no = 0, h*lambda = 3.295316889242648E+00 step(A) = 204, err = 4.303657217348121E-01, h_new = 3.575746941094147E-03, n_yes = 22, n_no = 0, h*lambda = 3.288891810450274E+00 step(A) = 205, err = 4.239645584035698E-01, h_new = 3.600118268548514E-03, n_yes = 23, n_no = 0, h*lambda = 3.294558162339284E+00 step(A) = 206, err = 4.258494625786690E-01, h_new = 3.619752878136625E-03, n_yes = 24, n_no = 0, h*lambda = 3.305257953938904E+00 step(A) = 207, err = 4.327561969283347E-01, h_new = 3.630198018866805E-03, n_yes = 25, n_no = 0, h*lambda = 3.311391416744631E+00 step(A) = 208, err = 4.381118801654675E-01, h_new = 3.635407523031155E-03, n_yes = 26, n_no = 0, h*lambda = 3.308973889359748E+00 step(A) = 209, err = 4.378033206390567E-01, h_new = 3.642852387740564E-03, n_yes = 27, n_no = 0, h*lambda = 3.301702492670826E+00 step(A) = 210, err = 4.331838707049168E-01, h_new = 3.656797891096875E-03, n_yes = 28, n_no = 0, h*lambda = 3.296381732091255E+00 step(A) = 211, err = 4.288082266575673E-01, h_new = 3.675577885545703E-03, n_yes = 29, n_no = 0, h*lambda = 3.296813684138294E+00 step(A) = 212, err = 4.279989888771147E-01, h_new = 3.694140400761753E-03, n_yes = 30, n_no = 0, h*lambda = 3.301420765469028E+00 step(A) = 213, err = 4.305553889676714E-01, h_new = 3.708759583319070E-03, n_yes = 31, n_no = 0, h*lambda = 3.305633300473522E+00 step(A) = 214, err = 4.336721825379390E-01, h_new = 3.719759708803672E-03, n_yes = 32, n_no = 0, h*lambda = 3.306143544613923E+00 step(A) = 215, err = 4.346492908611463E-01, h_new = 3.730441485950846E-03, n_yes = 33, n_no = 0, h*lambda = 3.303290092264038E+00 step(A) = 216, err = 4.330551307192517E-01, h_new = 3.743828608745091E-03, n_yes = 34, n_no = 0, h*lambda = 3.300029807125569E+00 step(A) = 217, err = 4.306225375706228E-01, h_new = 3.760310844181102E-03, n_yes = 35, n_no = 0, h*lambda = 3.299020735246466E+00 step(A) = 218, err = 4.293756907639733E-01, h_new = 3.777876527847335E-03, n_yes = 36, n_no = 0, h*lambda = 3.300566115806182E+00 step(A) = 219, err = 4.299588438174837E-01, h_new = 3.794208535276982E-03, n_yes = 37, n_no = 0, h*lambda = 3.302870523431511E+00 step(A) = 220, err = 4.314321127224678E-01, h_new = 3.808602621895213E-03, n_yes = 38, n_no = 0, h*lambda = 3.303907816511167E+00 step(A) = 221, err = 4.323389026423615E-01, h_new = 3.822209930449244E-03, n_yes = 39, n_no = 0, h*lambda = 3.303085320257956E+00 step(A) = 222, err = 4.319722586880699E-01, h_new = 3.836741351494489E-03, n_yes = 40, n_no = 0, h*lambda = 3.301419147277730E+00 step(A) = 223, err = 4.307927176047864E-01, h_new = 3.852987912337984E-03, n_yes = 41, n_no = 0, h*lambda = 3.300385095950699E+00 step(A) = 224, err = 4.298113342124661E-01, h_new = 3.870380413850536E-03, n_yes = 42, n_no = 0, h*lambda = 3.300644049791831E+00 step(A) = 225, err = 4.296631557704300E-01, h_new = 3.887724647013034E-03, n_yes = 43, n_no = 0, h*lambda = 3.301686993645893E+00 step(A) = 226, err = 4.301851388667979E-01, h_new = 3.904286808013262E-03, n_yes = 44, n_no = 0, h*lambda = 3.302485977631584E+00 step(A) = 227, err = 4.307133054066499E-01, h_new = 3.920292122979136E-03, n_yes = 45, n_no = 0, h*lambda = 3.302422090464993E+00 step(A) = 228, err = 4.307142557317793E-01, h_new = 3.936554777334662E-03, n_yes = 46, n_no = 0, h*lambda = 3.301694001216952E+00 step(A) = 229, err = 4.301851746333406E-01, h_new = 3.953711297252496E-03, n_yes = 47, n_no = 0, h*lambda = 3.300988215441240E+00 step(A) = 230, err = 4.295452210841031E-01, h_new = 3.971752422300344E-03, n_yes = 48, n_no = 0, h*lambda = 3.300828777040461E+00 step(A) = 231, err = 4.292107952935771E-01, h_new = 3.990166572077440E-03, n_yes = 49, n_no = 0, h*lambda = 3.301193521191673E+00 step(A) = 232, err = 4.292646244779778E-01, h_new = 4.008455751590996E-03, n_yes = 50, n_no = 0, h*lambda = 3.301646410466073E+00 step(A) = 233, err = 4.294605789942588E-01, h_new = 4.026536547814937E-03, n_yes = 51, n_no = 0, h*lambda = 3.301771849193995E+00 step(A) = 234, err = 4.294881484821736E-01, h_new = 4.044728598672606E-03, n_yes = 52, n_no = 0, h*lambda = 3.301501875902923E+00 step(A) = 235, err = 4.292334976052844E-01, h_new = 4.063422951925307E-03, n_yes = 53, n_no = 0, h*lambda = 3.301097579500324E+00 step(A) = 236, err = 4.288292849856661E-01, h_new = 4.082760731028622E-03, n_yes = 54, n_no = 0, h*lambda = 3.300871772815385E+00 step(A) = 237, err = 4.284966709701275E-01, h_new = 4.102577075871903E-03, n_yes = 55, n_no = 0, h*lambda = 3.300928638830731E+00 step(A) = 238, err = 4.283501725373395E-01, h_new = 4.122601297893762E-03, n_yes = 56, n_no = 0, h*lambda = 3.301124449621309E+00 step(A) = 239, err = 4.283289637297810E-01, h_new = 4.142701463051705E-03, n_yes = 57, n_no = 0, h*lambda = 3.301234132454486E+00 step(A) = 240, err = 4.282820596349022E-01, h_new = 4.162968884373310E-03, n_yes = 58, n_no = 0, h*lambda = 3.301147770350908E+00 step(A) = 241, err = 4.281074167924426E-01, h_new = 4.183607199895076E-03, n_yes = 59, n_no = 0, h*lambda = 3.300934116181848E+00 step(A) = 242, err = 4.278229899028322E-01, h_new = 4.204754278001979E-03, n_yes = 60, n_no = 0, h*lambda = 3.300748273610640E+00 step(A) = 243, err = 4.275270575766886E-01, h_new = 4.226393038144650E-03, n_yes = 61, n_no = 0, h*lambda = 3.300689871481256E+00 step(A) = 244, err = 4.273027857195540E-01, h_new = 4.248404526255284E-03, n_yes = 62, n_no = 0, h*lambda = 3.300735255137489E+00 step(A) = 245, err = 4.271550363186925E-01, h_new = 4.270692093309106E-03, n_yes = 63, n_no = 0, h*lambda = 3.300782193626054E+00 step(A) = 246, err = 4.270228587459070E-01, h_new = 4.293263069108798E-03, n_yes = 64, n_no = 0, h*lambda = 3.300747843562565E+00 step(A) = 247, err = 4.268417375909908E-01, h_new = 4.316211182245253E-03, n_yes = 65, n_no = 0, h*lambda = 3.300630716059762E+00 step(A) = 248, err = 4.265956213797722E-01, h_new = 4.339633801499685E-03, n_yes = 66, n_no = 0, h*lambda = 3.300494780817451E+00 step(A) = 249, err = 4.263193737176272E-01, h_new = 4.363563362692269E-03, n_yes = 67, n_no = 0, h*lambda = 3.300404709203970E+00 step(A) = 250, err = 4.260600881587809E-01, h_new = 4.387964990585588E-03, n_yes = 68, n_no = 0, h*lambda = 3.300374150266467E+00 step(A) = 251, err = 4.258368123725155E-01, h_new = 4.412788910192447E-03, n_yes = 69, n_no = 0, h*lambda = 3.300365536963663E+00 step(A) = 252, err = 4.256311074016821E-01, h_new = 4.438024743325478E-03, n_yes = 70, n_no = 0, h*lambda = 3.300330755475881E+00 step(A) = 253, err = 4.254094374485623E-01, h_new = 4.463713918291515E-03, n_yes = 71, n_no = 0, h*lambda = 3.300251688463127E+00 step(A) = 254, err = 4.251528637705740E-01, h_new = 4.489918712614180E-03, n_yes = 72, n_no = 0, h*lambda = 3.300148342366199E+00 step(A) = 255, err = 4.248688506907005E-01, h_new = 4.516681435733419E-03, n_yes = 73, n_no = 0, h*lambda = 3.300054879682836E+00 step(A) = 256, err = 4.245795292937601E-01, h_new = 4.544008415283115E-03, n_yes = 74, n_no = 0, h*lambda = 3.299989397775974E+00 step(A) = 257, err = 4.243007359497090E-01, h_new = 4.571886655870232E-03, n_yes = 75, n_no = 0, h*lambda = 3.299942855104354E+00 step(A) = 258, err = 4.240304753682130E-01, h_new = 4.600313341654090E-03, n_yes = 76, n_no = 0, h*lambda = 3.299891949478224E+00 step(A) = 259, err = 4.237540025421365E-01, h_new = 4.629312064912752E-03, n_yes = 77, n_no = 0, h*lambda = 3.299820658193305E+00 step(A) = 260, err = 4.234579704995718E-01, h_new = 4.658925511089058E-03, n_yes = 78, n_no = 0, h*lambda = 3.299731457982905E+00 step(A) = 261, err = 4.231404670605873E-01, h_new = 4.689195217083358E-03, n_yes = 79, n_no = 0, h*lambda = 3.299639449519019E+00 step(A) = 262, err = 4.228098096402270E-01, h_new = 4.720147236696066E-03, n_yes = 80, n_no = 0, h*lambda = 3.299557656336072E+00 step(A) = 263, err = 4.224754266122596E-01, h_new = 4.751794061992471E-03, n_yes = 81, n_no = 0, h*lambda = 3.299487064891211E+00 step(A) = 264, err = 4.221397572250009E-01, h_new = 4.784148090630616E-03, n_yes = 82, n_no = 0, h*lambda = 3.299418345989634E+00 step(A) = 265, err = 4.217974025208784E-01, h_new = 4.817233646248613E-03, n_yes = 83, n_no = 0, h*lambda = 3.299341434462446E+00 step(A) = 266, err = 4.214406700743320E-01, h_new = 4.851088314633154E-03, n_yes = 84, n_no = 0, h*lambda = 3.299253687082273E+00 step(A) = 267, err = 4.210657477187539E-01, h_new = 4.885754748734701E-03, n_yes = 85, n_no = 0, h*lambda = 3.299160419085479E+00 step(A) = 268, err = 4.206744564353846E-01, h_new = 4.921271494404488E-03, n_yes = 86, n_no = 0, h*lambda = 3.299068991663855E+00 step(A) = 269, err = 4.202711170789381E-01, h_new = 4.957670479629963E-03, n_yes = 87, n_no = 0, h*lambda = 3.298982513508011E+00 step(A) = 270, err = 4.198580540936551E-01, h_new = 4.994981975275775E-03, n_yes = 88, n_no = 0, h*lambda = 3.298898254205975E+00 step(A) = 271, err = 4.194335995255202E-01, h_new = 5.033241715662309E-03, n_yes = 89, n_no = 0, h*lambda = 3.298811045041153E+00 step(A) = 272, err = 4.189936537320252E-01, h_new = 5.072494208865829E-03, n_yes = 90, n_no = 0, h*lambda = 3.298717913853457E+00 step(A) = 273, err = 4.185348058702847E-01, h_new = 5.112790510439460E-03, n_yes = 91, n_no = 0, h*lambda = 3.298620005927656E+00 step(A) = 274, err = 4.180561978342951E-01, h_new = 5.154183518165018E-03, n_yes = 92, n_no = 0, h*lambda = 3.298520832098713E+00 step(A) = 275, err = 4.175589093278980E-01, h_new = 5.196725241731624E-03, n_yes = 93, n_no = 0, h*lambda = 3.298423010402977E+00 step(A) = 276, err = 4.170438715306105E-01, h_new = 5.240468066169636E-03, n_yes = 94, n_no = 0, h*lambda = 3.298326412487723E+00 step(A) = 277, err = 4.165102735711282E-01, h_new = 5.285468464512027E-03, n_yes = 95, n_no = 0, h*lambda = 3.298228897593983E+00 step(A) = 278, err = 4.159556178380689E-01, h_new = 5.331789991909682E-03, n_yes = 96, n_no = 0, h*lambda = 3.298128510128849E+00 step(A) = 279, err = 4.153770237895118E-01, h_new = 5.379503623611202E-03, n_yes = 97, n_no = 0, h*lambda = 3.298025088501638E+00 step(A) = 280, err = 4.147724711091176E-01, h_new = 5.428686039233089E-03, n_yes = 98, n_no = 0, h*lambda = 3.297920132091144E+00 step(A) = 281, err = 4.141409969555077E-01, h_new = 5.479418021445208E-03, n_yes = 99, n_no = 0, h*lambda = 3.297815401152697E+00 step(A) = 282, err = 4.134819002790782E-01, h_new = 5.531784671826462E-03, n_yes = 100, n_no = 0, h*lambda = 3.297711633507199E+00 step(A) = 283, err = 4.127937823687496E-01, h_new = 5.585877421317612E-03, n_yes = 101, n_no = 0, h*lambda = 3.297608380073231E+00 step(A) = 284, err = 4.120742054534275E-01, h_new = 5.641796452121301E-03, n_yes = 102, n_no = 0, h*lambda = 3.297504863095371E+00 step(A) = 285, err = 4.113200925774184E-01, h_new = 5.699652167295843E-03, n_yes = 103, n_no = 0, h*lambda = 3.297400969642974E+00 step(A) = 286, err = 4.105283746595030E-01, h_new = 5.759565460900966E-03, n_yes = 104, n_no = 0, h*lambda = 3.297297572992058E+00 step(A) = 287, err = 4.096962958683515E-01, h_new = 5.821667654197287E-03, n_yes = 105, n_no = 0, h*lambda = 3.297196079676924E+00 step(A) = 288, err = 4.088211805508948E-01, h_new = 5.886101189535850E-03, n_yes = 106, n_no = 0, h*lambda = 3.297097739264321E+00 step(A) = 289, err = 4.078999292805139E-01, h_new = 5.953021509572973E-03, n_yes = 107, n_no = 0, h*lambda = 3.297003333553935E+00 step(A) = 290, err = 4.069286643367531E-01, h_new = 6.022599701670972E-03, n_yes = 108, n_no = 0, h*lambda = 3.296913436259845E+00 step(A) = 291, err = 4.059027307932550E-01, h_new = 6.095025172835020E-03, n_yes = 109, n_no = 0, h*lambda = 3.296828947092619E+00 step(A) = 292, err = 4.048169258526370E-01, h_new = 6.170507994936443E-03, n_yes = 110, n_no = 0, h*lambda = 3.296751455234593E+00 step(A) = 293, err = 4.036656653309485E-01, h_new = 6.249281203922641E-03, n_yes = 111, n_no = 0, h*lambda = 3.296683224161232E+00 step(A) = 294, err = 4.024428954343365E-01, h_new = 6.331603700551526E-03, n_yes = 112, n_no = 0, h*lambda = 3.296626944716787E+00 step(A) = 295, err = 4.011417849609941E-01, h_new = 6.417764266935371E-03, n_yes = 113, n_no = 0, h*lambda = 3.296585571191088E+00 step(A) = 296, err = 3.997543801402357E-01, h_new = 6.508086807560819E-03, n_yes = 114, n_no = 0, h*lambda = 3.296562443275698E+00 step(A) = 297, err = 3.982713709433062E-01, h_new = 6.602936651123937E-03, n_yes = 115, n_no = 0, h*lambda = 3.296561653745240E+00 step(A) = 298, err = 3.966819670435605E-01, h_new = 6.702727844177274E-03, n_yes = 116, n_no = 0, h*lambda = 3.296588469298817E+00 step(A) = 299, err = 3.949737553404244E-01, h_new = 6.807931743832539E-03, n_yes = 117, n_no = 0, h*lambda = 3.296649654825813E+00 step(A) = 300, err = 3.931324002247971E-01, h_new = 6.919087605637563E-03, n_yes = 118, n_no = 0, h*lambda = 3.296753719325883E+00 step(A) = 301, err = 3.911411317652012E-01, h_new = 7.036816069002807E-03, n_yes = 119, n_no = 0, h*lambda = 3.296911242802162E+00 step(A) = 302, err = 3.889800524912123E-01, h_new = 7.161836500792817E-03, n_yes = 120, n_no = 0, h*lambda = 3.297135461658360E+00 step(A) = 303, err = 3.866253029065805E-01, h_new = 7.294989285968968E-03, n_yes = 121, n_no = 0, h*lambda = 3.297443211238030E+00 step(A) = 304, err = 3.840480530729704E-01, h_new = 7.437264588288742E-03, n_yes = 122, n_no = 0, h*lambda = 3.297856252398439E+00 step(A) = 305, err = 3.812131839851460E-01, h_new = 7.589839976134515E-03, n_yes = 123, n_no = 0, h*lambda = 3.298403033531530E+00 step(A) = 306, err = 3.780774391405725E-01, h_new = 7.754130675266364E-03, n_yes = 124, n_no = 0, h*lambda = 3.299121075776918E+00 step(A) = 307, err = 3.745867642471649E-01, h_new = 7.931858222387527E-03, n_yes = 125, n_no = 0, h*lambda = 3.300060376571410E+00 step(A) = 308, err = 3.706724549957577E-01, h_new = 8.125146419276591E-03, n_yes = 126, n_no = 0, h*lambda = 3.301288478614701E+00 step(A) = 309, err = 3.662455075628569E-01, h_new = 8.336658774680046E-03, n_yes = 127, n_no = 0, h*lambda = 3.302898201715120E+00 step(A) = 310, err = 3.611880979145950E-01, h_new = 8.569801139885829E-03, n_yes = 128, n_no = 0, h*lambda = 3.305019656714996E+00 step(A) = 311, err = 3.553402126001545E-01, h_new = 8.829031172535121E-03, n_yes = 129, n_no = 0, h*lambda = 3.307839385840251E+00 step(A) = 312, err = 3.484776836108640E-01, h_new = 9.120351759979481E-03, n_yes = 130, n_no = 0, h*lambda = 3.311631938628216E+00 step(A) = 313, err = 3.402741534167765E-01, h_new = 9.452140609354995E-03, n_yes = 131, n_no = 0, h*lambda = 3.316814291706466E+00 step(A) = 314, err = 3.302309183282819E-01, h_new = 9.836641122379494E-03, n_yes = 132, n_no = 0, h*lambda = 3.324044700868102E+00 step(A) = 315, err = 3.175367833680536E-01, h_new = 1.029288353523805E-02, n_yes = 133, n_no = 0, h*lambda = 3.334414182583914E+00 step(A) = 316, err = 3.007574192129444E-01, h_new = 1.085311863759427E-02, n_yes = 134, n_no = 0, h*lambda = 3.349848960293911E+00 step(A) = 317, err = 2.770425480094375E-01, h_new = 1.157958015835199E-02, n_yes = 135, n_no = 0, h*lambda = 3.374052218156454E+00 step(A) = 318, err = 2.396474186032754E-01, h_new = 1.262146633910143E-02, n_yes = 136, n_no = 0, h*lambda = 3.415050308842849E+00 step(A) = 319, err = 1.672169330518306E-01, h_new = 1.454044138087391E-02, n_yes = 137, n_no = 0, h*lambda = 3.493245381393254E+00 step(A) = 320, err = 6.547522351608699E-02, h_new = 1.936503210121610E-02, n_yes = 138, n_no = 0, h*lambda = 3.779558822626685E+00 step(R) = 321, err = 2.758072167696834E+00, h_new = 1.466755665203884E-02 step(A) = 322, err = 5.544936042036211E-01, h_new = 1.308521766712097E-02, n_yes = 139, n_no = 0, h*lambda = 3.439844267455001E+00 step(A) = 323, err = 4.608234807061337E-01, h_new = 1.312134759461917E-02, n_yes = 139, n_no = 1, h*lambda = 2.796465720743475E+00 step(A) = 324, err = 2.969699612327859E-01, h_new = 1.407347043448586E-02, n_yes = 139, n_no = 2, h*lambda = 2.524954751929724E+00 step(A) = 325, err = 2.801906534171543E-01, h_new = 1.497907847147494E-02, n_yes = 139, n_no = 3, h*lambda = 2.369860469525626E+00 step(A) = 326, err = 2.705027714409222E-01, h_new = 1.600134757536991E-02, n_yes = 139, n_no = 4, h*lambda = 2.114647776142334E+00 step(A) = 327, err = 2.322671428927122E-01, h_new = 1.751733417013966E-02, n_yes = 139, n_no = 5, h*lambda = 1.767476442444690E+00 step(A) = 328, err = 1.694881217518184E-01, h_new = 2.010927253694785E-02, n_yes = 0, n_no = 6, h*lambda = 1.306442795376463E+00 step(R) = 329, err = 1.084461510298132E+00, h_new = 1.785058548720665E-02 step(A) = 330, err = 4.164041221219500E-01, h_new = 1.736772732164151E-02, n_yes = 0, n_no = 7, h*lambda = 1.821084854705834E-01 step(R) = 331, err = 5.013796320452609E+00, h_new = 1.188385008552715E-02 step(R) = 332, err = 2.939269588993492E+00, h_new = 8.904283366664307E-03 step(A) = 333, err = 6.344689183230303E-01, h_new = 8.360116372816851E-03, n_yes = 0, n_no = 8, h*lambda = 4.999577635130124E-01 step(R) = 334, err = 1.506978358142714E+00, h_new = 7.017406578821495E-03 step(R) = 335, err = 1.507078084970641E+00, h_new = 5.890281776140481E-03 step(R) = 336, err = 1.056204469623091E+00, h_new = 5.252201992389361E-03 step(A) = 337, err = 7.013961933120292E-01, h_new = 4.930223016397128E-03, n_yes = 0, n_no = 9, h*lambda = 4.458951238252798E-01 step(R) = 338, err = 1.368003160369734E+01, h_new = 2.844292572253510E-03 step(A) = 339, err = 2.154297042023680E-01, h_new = 3.276373188376880E-03, n_yes = 0, n_no = 10, h*lambda = 5.781008123990621E-01 step(A) = 340, err = 6.517791037248070E-01, h_new = 2.589118491383514E-03, n_yes = 1, n_no = 0, h*lambda = 9.662911863007947E+00 step(R) = 341, err = 9.363066039698422E+00, h_new = 1.593137869750087E-03 step(R) = 342, err = 1.103967583380753E+00, h_new = 1.409916185041910E-03 step(A) = 343, err = 5.792252937311234E-01, h_new = 1.368723565504569E-03, n_yes = 1, n_no = 1, h*lambda = 6.116552865099667E-01 step(A) = 344, err = 6.875690543096030E-01, h_new = 1.284483056994523E-03, n_yes = 1, n_no = 2, h*lambda = 1.225537118258073E+00 step(R) = 345, err = 2.533668995469276E+00, h_new = 9.870369488894908E-04 step(A) = 346, err = 5.634250690067383E-01, h_new = 9.647758513983588E-04, n_yes = 1, n_no = 3, h*lambda = 9.696329291321375E-01 step(A) = 347, err = 6.397646579896418E-01, h_new = 9.155434055823949E-04, n_yes = 1, n_no = 4, h*lambda = 9.695144007444304E-01 step(A) = 348, err = 5.120257578337231E-01, h_new = 9.069485495246315E-04, n_yes = 1, n_no = 5, h*lambda = 9.258268721123650E-01 step(A) = 349, err = 4.684576636521415E-01, h_new = 9.040300635595680E-04, n_yes = 0, n_no = 6, h*lambda = 9.192261126382347E-01 step(A) = 350, err = 3.992343733368792E-01, h_new = 9.226637306866363E-04, n_yes = 0, n_no = 7, h*lambda = 9.166675883651417E-01 step(A) = 351, err = 1.739001318661580E-01, h_new = 1.077665775295479E-03, n_yes = 0, n_no = 8, h*lambda = 9.352847703090210E-01 step(A) = 352, err = 1.450216591933589E-01, h_new = 1.255726353495396E-03, n_yes = 0, n_no = 9, h*lambda = 1.091654665463814E+00 step(A) = 353, err = 1.079829569870686E-01, h_new = 1.527300332301485E-03, n_yes = 0, n_no = 10, h*lambda = 1.270717098143420E+00 step(A) = 354, err = 8.757978699270973E-02, h_new = 1.902359116371433E-03, n_yes = 0, n_no = 11, h*lambda = 1.543482436580736E+00 step(A) = 355, err = 6.456233631723808E-02, h_new = 2.474769177299382E-03, n_yes = 0, n_no = 12, h*lambda = 1.919294523613246E+00 step(A) = 356, err = 4.949968416402232E-02, h_new = 3.327315357578232E-03, n_yes = 0, n_no = 13, h*lambda = 2.491356651671514E+00 step(A) = 357, err = 6.438796693672601E-02, h_new = 4.232761223977908E-03, n_yes = 1, n_no = 0, h*lambda = 3.339793883858084E+00 step(A) = 358, err = 2.797730075537238E-01, h_new = 4.238981097626998E-03, n_yes = 2, n_no = 0, h*lambda = 4.232766260040452E+00 step(R) = 359, err = 1.316524838249334E+00, h_new = 3.640835437870546E-03 step(A) = 360, err = 5.341844412925590E-01, h_new = 3.464242470661683E-03, n_yes = 3, n_no = 0, h*lambda = 3.629137340689636E+00 step(A) = 361, err = 7.018131540087504E-01, h_new = 3.229243431640143E-03, n_yes = 4, n_no = 0, h*lambda = 3.442495298916736E+00 step(A) = 362, err = 5.894036583006134E-01, h_new = 3.134887599014504E-03, n_yes = 4, n_no = 1, h*lambda = 3.199693006923365E+00 step(A) = 363, err = 4.026123116021955E-01, h_new = 3.224409542044754E-03, n_yes = 4, n_no = 2, h*lambda = 3.097414849933740E+00 step(A) = 364, err = 3.179969364900750E-01, h_new = 3.399979433851089E-03, n_yes = 4, n_no = 3, h*lambda = 3.176578537857591E+00 step(A) = 365, err = 3.367762714656367E-01, h_new = 3.516961908380592E-03, n_yes = 5, n_no = 0, h*lambda = 3.339234748250505E+00 step(A) = 366, err = 4.309276422965480E-01, h_new = 3.496673589755197E-03, n_yes = 6, n_no = 0, h*lambda = 3.443108463539667E+00 step(A) = 367, err = 5.277919916085688E-01, h_new = 3.391994204439725E-03, n_yes = 7, n_no = 0, h*lambda = 3.412353949788284E+00 step(A) = 368, err = 5.296304036466514E-01, h_new = 3.315283929811912E-03, n_yes = 8, n_no = 0, h*lambda = 3.299927072673363E+00 step(A) = 369, err = 4.526156127222287E-01, h_new = 3.328496725446670E-03, n_yes = 8, n_no = 1, h*lambda = 3.215457491426718E+00 step(A) = 370, err = 3.869624821910816E-01, h_new = 3.410484801370306E-03, n_yes = 8, n_no = 2, h*lambda = 3.218340142079507E+00 step(A) = 371, err = 3.744993934307922E-01, h_new = 3.492035838156108E-03, n_yes = 9, n_no = 0, h*lambda = 3.287189159795818E+00 step(A) = 372, err = 4.100194386256036E-01, h_new = 3.516272226636235E-03, n_yes = 10, n_no = 0, h*lambda = 3.354865484299512E+00 step(A) = 373, err = 4.608766850419500E-01, h_new = 3.483595911389426E-03, n_yes = 11, n_no = 0, h*lambda = 3.367069048243797E+00 step(A) = 374, err = 4.816281309864093E-01, h_new = 3.441538699119594E-03, n_yes = 12, n_no = 0, h*lambda = 3.324891144951819E+00 step(A) = 375, err = 4.584277456682849E-01, h_new = 3.434690339403088E-03, n_yes = 13, n_no = 0, h*lambda = 3.274104202325452E+00 step(A) = 376, err = 4.217010246802865E-01, h_new = 3.470005332839345E-03, n_yes = 14, n_no = 0, h*lambda = 3.256970425629175E+00 THE PROBLEM SEEMS TO BECOME STIFF AT X = 1.744401291692531E+00, ACCEPTED STEP = 357 step(A) = 377, err = 4.035671505294718E-01, h_new = 3.520198475311247E-03, n_yes = 15, n_no = 0, h*lambda = 3.279612978036328E+00 step(A) = 378, err = 4.124998511616665E-01, h_new = 3.551601632330066E-03, n_yes = 16, n_no = 0, h*lambda = 3.315887769755170E+00 step(A) = 379, err = 4.364865333541689E-01, h_new = 3.552128464250999E-03, n_yes = 17, n_no = 0, h*lambda = 3.334099100813723E+00 step(A) = 380, err = 4.537348839191425E-01, h_new = 3.537314207101296E-03, n_yes = 18, n_no = 0, h*lambda = 3.323211423181647E+00 step(A) = 381, err = 4.507776443118724E-01, h_new = 3.531950690399045E-03, n_yes = 19, n_no = 0, h*lambda = 3.298050238123096E+00 step(A) = 382, err = 4.342894091804872E-01, h_new = 3.548078055167689E-03, n_yes = 20, n_no = 0, h*lambda = 3.281767755345274E+00 step(A) = 383, err = 4.206542956351669E-01, h_new = 3.578322989303431E-03, n_yes = 21, n_no = 0, h*lambda = 3.285356720285332E+00 step(A) = 384, err = 4.196258972151440E-01, h_new = 3.605723937723291E-03, n_yes = 22, n_no = 0, h*lambda = 3.301762969345275E+00 step(A) = 385, err = 4.292224611770568E-01, h_new = 3.619040631090450E-03, n_yes = 23, n_no = 0, h*lambda = 3.315260929917143E+00 step(A) = 386, err = 4.397641151446334E-01, h_new = 3.620728012934108E-03, n_yes = 24, n_no = 0, h*lambda = 3.315622571336352E+00 step(A) = 387, err = 4.424148223976933E-01, h_new = 3.622231135717094E-03, n_yes = 25, n_no = 0, h*lambda = 3.305262547419302E+00 step(A) = 388, err = 4.365197194349388E-01, h_new = 3.632881216517289E-03, n_yes = 26, n_no = 0, h*lambda = 3.294705192440800E+00 step(A) = 389, err = 4.286360496829625E-01, h_new = 3.652908418797025E-03, n_yes = 27, n_no = 0, h*lambda = 3.292379763151353E+00 step(A) = 390, err = 4.253576036254677E-01, h_new = 3.675163184206087E-03, n_yes = 28, n_no = 0, h*lambda = 3.298373393203926E+00 step(A) = 391, err = 4.281672664980711E-01, h_new = 3.692283301842687E-03, n_yes = 29, n_no = 0, h*lambda = 3.306152666725394E+00 step(A) = 392, err = 4.335203958450980E-01, h_new = 3.702631088355893E-03, n_yes = 30, n_no = 0, h*lambda = 3.309111987348828E+00 step(A) = 393, err = 4.365379891511153E-01, h_new = 3.710475656904800E-03, n_yes = 31, n_no = 0, h*lambda = 3.305861487498091E+00 step(A) = 394, err = 4.351678259500148E-01, h_new = 3.721356918232027E-03, n_yes = 32, n_no = 0, h*lambda = 3.300274009212814E+00 step(A) = 395, err = 4.313533339230897E-01, h_new = 3.737390428949686E-03, n_yes = 33, n_no = 0, h*lambda = 3.297272979759182E+00 step(A) = 396, err = 4.285761795087774E-01, h_new = 3.756293683258537E-03, n_yes = 34, n_no = 0, h*lambda = 3.298677130890415E+00 step(A) = 397, err = 4.287190833217491E-01, h_new = 3.774103377896347E-03, n_yes = 35, n_no = 0, h*lambda = 3.302416479113768E+00 step(A) = 398, err = 4.309643259122914E-01, h_new = 3.788682301899664E-03, n_yes = 36, n_no = 0, h*lambda = 3.304993020455877E+00 step(A) = 399, err = 4.329973638391673E-01, h_new = 3.801069921999938E-03, n_yes = 37, n_no = 0, h*lambda = 3.304563353713653E+00 step(A) = 400, err = 4.331357163464960E-01, h_new = 3.814008869784482E-03, n_yes = 38, n_no = 0, h*lambda = 3.302070284168640E+00 step(A) = 401, err = 4.315662949195275E-01, h_new = 3.829403146181749E-03, n_yes = 39, n_no = 0, h*lambda = 3.299906818493146E+00 step(A) = 402, err = 4.298065745571851E-01, h_new = 3.846972490640173E-03, n_yes = 40, n_no = 0, h*lambda = 3.299698514896586E+00 step(A) = 403, err = 4.291955448638680E-01, h_new = 3.864925504721692E-03, n_yes = 41, n_no = 0, h*lambda = 3.301170400615810E+00 step(A) = 404, err = 4.298730042492989E-01, h_new = 3.881700430530588E-03, n_yes = 42, n_no = 0, h*lambda = 3.302765846451996E+00 step(A) = 405, err = 4.309260669387576E-01, h_new = 3.897172792658790E-03, n_yes = 43, n_no = 0, h*lambda = 3.303154342426523E+00 step(A) = 406, err = 4.313284316802334E-01, h_new = 3.912468980688625E-03, n_yes = 44, n_no = 0, h*lambda = 3.302246150455781E+00 step(A) = 407, err = 4.307845084319388E-01, h_new = 3.928814528533372E-03, n_yes = 45, n_no = 0, h*lambda = 3.301004969890132E+00 step(A) = 408, err = 4.298165423303336E-01, h_new = 3.946538175855745E-03, n_yes = 46, n_no = 0, h*lambda = 3.300457272966248E+00 step(A) = 409, err = 4.291720917134020E-01, h_new = 3.964996354060733E-03, n_yes = 47, n_no = 0, h*lambda = 3.300860336245485E+00 step(A) = 410, err = 4.291798750075976E-01, h_new = 3.983289498679896E-03, n_yes = 48, n_no = 0, h*lambda = 3.301659001676910E+00 step(A) = 411, err = 4.295785616107266E-01, h_new = 4.001038338066340E-03, n_yes = 49, n_no = 0, h*lambda = 3.302098366435485E+00 step(A) = 412, err = 4.298476998635363E-01, h_new = 4.018587629621689E-03, n_yes = 50, n_no = 0, h*lambda = 3.301867531751385E+00 step(A) = 413, err = 4.296760459012742E-01, h_new = 4.036589093769783E-03, n_yes = 51, n_no = 0, h*lambda = 3.301254848796467E+00 step(A) = 414, err = 4.291671457167899E-01, h_new = 4.055423356777790E-03, n_yes = 52, n_no = 0, h*lambda = 3.300792187767611E+00 step(A) = 415, err = 4.286692624819732E-01, h_new = 4.074956413624082E-03, n_yes = 53, n_no = 0, h*lambda = 3.300779491762202E+00 step(A) = 416, err = 4.284453888845481E-01, h_new = 4.094757061906611E-03, n_yes = 54, n_no = 0, h*lambda = 3.301094355125355E+00 step(A) = 417, err = 4.284820221704877E-01, h_new = 4.114508142343789E-03, n_yes = 55, n_no = 0, h*lambda = 3.301378208179972E+00 step(A) = 418, err = 4.285581027441299E-01, h_new = 4.134243848945743E-03, n_yes = 56, n_no = 0, h*lambda = 3.301371663800388E+00 step(A) = 419, err = 4.284643391207051E-01, h_new = 4.154258248987697E-03, n_yes = 57, n_no = 0, h*lambda = 3.301100608812401E+00 step(A) = 420, err = 4.281676627059799E-01, h_new = 4.174824569557473E-03, n_yes = 58, n_no = 0, h*lambda = 3.300795749714216E+00 step(A) = 421, err = 4.278005791273590E-01, h_new = 4.195988237915218E-03, n_yes = 59, n_no = 0, h*lambda = 3.300665609574369E+00 step(A) = 422, err = 4.275203746049487E-01, h_new = 4.217584256108569E-03, n_yes = 60, n_no = 0, h*lambda = 3.300734901041688E+00 step(A) = 423, err = 4.273784348054304E-01, h_new = 4.239419633725231E-03, n_yes = 61, n_no = 0, h*lambda = 3.300863388402065E+00 step(A) = 424, err = 4.273015686248449E-01, h_new = 4.261441761634310E-03, n_yes = 62, n_no = 0, h*lambda = 3.300892260934269E+00 step(A) = 425, err = 4.271763403969597E-01, h_new = 4.283760915589911E-03, n_yes = 63, n_no = 0, h*lambda = 3.300774630150715E+00 step(A) = 426, err = 6.964951300427906E-03, h_new = 4.190371271724428E-03, n_yes = 63, n_no = 1, h*lambda = 1.599612021680380E+00 DoPri5: Dormand-Prince method (explicit, order 5(4), embedded) Number of function evaluations = 2558 Number of performed steps = 426 Number of accepted steps = 406 Number of rejected steps = 20 y = 1.820788982019278E+00 -7.853646714272298E-01 h = 4.190371271724428E-03