Machine precision is 1.1102230246251568e-016 Beginning Adams solver root test. 0 0.5 0.5 0.62831853071795862 2 1.5493077741038221 0.4054803768614112 0.30661963271977283 4 1.742294596429921 0.3743267212643196 -0.12890384174769798 6 1.0055382397944168 0.41730609001451435 -0.55067965874121361 7.2883402424829775 -3.1849523018934178e-015 0.47485631008838808 -0.76010665852632853 Should have found a root near 7.288 Beginning Adams solver no root test. 0 0.5 0.5 0.62831853071795862 2 1.5493077741038221 0.4054803768614112 0.30661963271977283 4 1.742294596429921 0.3743267212643196 -0.12890384174769798 6 1.0055382397944168 0.41730609001451435 -0.55067965874121361 8 -0.74588732656237355 0.51298719385415892 -0.85370606303034191 10 -3.6281304248930493 0.63325477046925893 -1.0515213144003499 Should have proceeded to 10.0 Beginning stiff solver non-stiff system test. 0 0.5 0.5 0.62831853071795862 2 1.5491098637141683 0.40553383184262431 0.30677198609920253 4 1.7430215367391404 0.37424350216283175 -0.12832013367503328 6 1.0070521266131176 0.41707652647004545 -0.55073254604270616 7.2895041873210111 3.2959746043559335e-016 0.47475519527746995 -0.76066407404601255 Should have found a root near 7.288 Beginning no root stiff solver non-stiff system test. 0 0.5 0.5 0.62831853071795862 2 1.5491098637141683 0.40553383184262431 0.30677198609920253 4 1.7430215367391404 0.37424350216283175 -0.12832013367503328 6 1.0070521266131176 0.41707652647004545 -0.55073254604270616 8 -0.74530990578378264 0.51288638789853258 -0.85396196163898208 10 -3.630224722965063 0.63328305798451801 -1.0517517951045769 Should have proceeded to 10.0 Beginning stiff solver stiff system test. 0 1 0 0 1 0.96645974426149051 3.0746226088749535e-005 0.033509509512372422 2 0.94160760712206071 2.7017526345273639e-005 0.058365375351448347 3 0.92188160593158297 2.4382910291918935e-005 0.078094011158098137 4 0.90551337261870146 2.2404011190902854e-005 0.094464223370076497 4.376713332739163 0.90000000000000002 2.1777805440668583e-005 0.099978222194517519 Should have found y(1)=0.9 near 4.377 ...and with jacobian... 0 1 0 0 1 0.96645974446465943 3.0746226145534692e-005 0.033509509309195093 2 0.94160760717450132 2.7017526349945626e-005 0.058365375299148606 3 0.92188160604793046 2.4382910317829561e-005 0.078094011041751593 4 0.90551337264506637 2.240401119265076e-005 0.094464223343740827 4.3767133353569809 0.90000000000000002 2.177780544213436e-005 0.099978222194557653 Should have found y(1)=0.9 near 4.377 Beginning RK range test. 0 0 1 0.78539816339744828 0.70708128465267572 0.70707636809260788 1.5707963267948966 0.99991993503028176 -7.6766758265595092e-006 2.3561944901923448 0.70701422156581017 -0.70702994584986811 3.1415926535897931 -1.5285283969407315e-005 -0.99984016872966364 3.9269908169872414 -0.70697894799349914 -0.70695248125718313 4.7123889803846897 -0.99976052309909558 2.2876560670636453e-005 5.497787143782138 -0.70689080486100586 0.70692800364611252 6.2831853071795862 3.0448991355350868e-005 0.99968119828793611 Last should be 6.283 0.0 1.0 Beginning RK one_step test. 3.7058605534259774 -0.53480184315256474 -0.84496482061944178 Should have gone to Pi so far. 4.3639400350516979 -0.93989698042059167 -0.34142430439772253 5.3201474586578206 -0.82089482364598421 0.57104018155859781