An Accurate Block Solver for Stiff Initial Value Problems

<table>Graph of <svg height="17.65" id="M200" style="vertical-align:-3.2955pt" version="1.1" viewbox="0 0 86.875 17.65" width="86.875" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.625,11.925)"><path d="M495 163l30 -6q-19 -99 -36 -157h-453v28q63 5 77.5 20t14.5 78v400q0 63 -13 77t-75 19v28h261v-28q-61 -5 -74.5 -19t-13.5 -77v-395q0 -41 5 -59t18 -25q24 -13 96 -13q85 0 112 31q30 33 51 98z" id="x4C"></path></g><g transform="matrix(.017,-0,0,-.017,9.906,11.925)"><path d="M257 449q92 0 154 -65t62 -158q0 -112 -67 -175t-150 -63q-98 0 -158.5 66.5t-60.5 154.5q0 59 21 106.5t54.5 75.5t71 43t73.5 15zM244 416q-48 0 -81 -47t-33 -128q0 -96 38 -158t99 -62q51 0 82 43.5t31 139.5q0 91 -36 151.5t-100 60.5z" id="x6F"></path></g><g transform="matrix(.017,-0,0,-.017,18.576,11.925)"><path d="M463 437l5 -8q-21 -38 -45 -51l-57 3q30 -35 30 -88q0 -76 -52.5 -118t-122.5 -42q-24 0 -51 5q-30 -21 -30 -42q0 -16 16 -27t47 -11q19 0 54 0.5t52 0.5q36 0 65.5 -8t53.5 -36t24 -74q0 -77 -75 -137.5t-174 -60.5q-79 0 -126 37.5t-48 84.5q0 32 24 56q21 21 83 70
q-49 15 -74 55q-13 20 -11 42q56 21 88 61q-37 16 -62.5 51t-25.5 83q0 78 56.5 122t122.5 44h1q51 0 96 -27q71 4 136 15zM219 418q-34 0 -58 -32t-24 -87q0 -59 26.5 -96.5t67.5 -37.5q35 1 59 32t24 87q0 59 -26 96.5t-69 37.5zM241 -11q-62 0 -87 -21q-46 -39 -46 -78
q0 -41 38 -68.5t94 -27.5q69 0 106.5 31t37.5 80q0 52 -45 71q-33 13 -98 13z" id="x67"></path></g>
<g transform="matrix(.012,-0,0,-.012,26.55,16.888)"><path d="M384 0h-275v27q67 5 81.5 18.5t14.5 68.5v385q0 38 -7.5 47.5t-40.5 10.5l-48 2v24q85 15 178 52v-521q0 -55 14.5 -68.5t82.5 -18.5v-27z" id="x31"></path></g><g transform="matrix(.012,-0,0,-.012,32.261,16.888)"><path d="M241 635q53 0 94 -28.5t63.5 -76t33.5 -102.5t11 -116q0 -58 -11 -112.5t-34 -103.5t-63.5 -78.5t-94.5 -29.5t-95 28t-64.5 75t-34.5 102.5t-11 118.5q0 58 11.5 112.5t34.5 103t64.5 78t95.5 29.5zM238 602q-32 0 -55.5 -25t-35.5 -68t-17.5 -91t-5.5 -105
q0 -76 10 -138.5t37 -107.5t69 -45q32 0 55.5 25t35.5 68.5t17.5 91.5t5.5 105t-5.5 105.5t-18 92t-36 68t-56.5 24.5z" id="x30"></path></g>
<g transform="matrix(.017,-0,0,-.017,38.588,11.925)"><path d="M861 0h-263v28q61 5 73 22t12 97l-4 408h-4l-244 -548h-20l-224 534h-2l-15 -285q-5 -98 -3 -166q1 -34 19.5 -45t72.5 -17v-28h-233v28q45 5 62.5 17t22.5 43q11 65 19 183l14 208q3 36 3.5 55t-2.5 37.5t-6.5 25t-16.5 13t-25.5 8t-39.5 4.5v28h164l228 -501l232 501
h170v-28q-68 -6 -81 -21.5t-11 -91.5l9 -362q2 -80 15 -97t78 -22v-28z" id="x4D"></path></g><g transform="matrix(.017,-0,0,-.017,53.733,11.925)"><path d="M673 0h-245v28q47 6 58 16t1 39l-55 154h-213q-31 -87 -36 -129q-16 -47 -5 -61t68 -19v-28h-226v28q50 7 68.5 23.5t43.5 80.5l214 524l24 9l193 -536q23 -63 41.5 -79.5t68.5 -21.5v-28zM418 280l-81 255h-2q-74 -197 -101 -255h184z" id="x41"></path></g><g transform="matrix(.017,-0,0,-.017,65.48,11.925)"><path d="M655 0h-263v28q46 5 55.5 14.5t-5.5 34.5l-132 214q-69 -104 -122 -202q-19 -32 -10.5 -43.5t58.5 -17.5v-28h-229v28q44 5 70 17t47 40q47 59 164 239l-134 214q-33 50 -54 64t-74 20v28h246v-28q-44 -5 -52 -14t5 -28q32 -55 114 -182q66 102 99 162q19 33 10.5 44.5
t-57.5 17.5v28h228v-28q-56 -6 -80 -19.5t-58 -57.5q-45 -62 -121 -179l157 -246q35 -53 60.5 -70t77.5 -22v-28z" id="x58"></path></g><g transform="matrix(.017,-0,0,-.017,76.597,11.925)"><path d="M517 162l29 -4q-18 -100 -35 -158h-482v28q68 5 83.5 20t15.5 77v398q0 66 -13 80t-73 19v28h446q0 -8 8 -144l-29 -4q-11 51 -32 85q-11 17 -35 23t-78 6h-76q-23 0 -28 -5t-5 -27v-228h94q67 0 83 12t25 63h29v-191h-29q-10 54 -26 65.5t-81 11.5h-95v-189
q0 -65 24 -81q19 -13 104 -13q57 0 83 7t43 27q28 33 50 94z" id="x45"></path></g>
</svg> versus TOL for problem (1).</table>

International Scholarly Research Notices

An Accurate Block Solver for Stiff Initial Value Problems

Figure 4