###################################################################### # # A better `efficient' Runge--Kutta (16:8,9) pair # # These are approximate REAL coefficients computed # using MAPLE with 40 digits for # # a slightly more efficient SIXTEEN-stage conventional # pair of methods of orders p=8 and p=9, with # dominant stage order = 5, # # together with approximate 40-digit coefficients # for two interpolants of orders 8 and 9 that # require 5 and 5 extra stages respectively. # # Exact coefficients of the method are possible using # surds in terms of sqrt(6). Exact coefficients for the # interpolants require that c[18] must be a zero of a # cubic polynomial with coefficients that are surds in # terms of sqrt(6). Hence, exact coefficients for the # interpolants cannot be conveniently represented. # # This procedure is a better efficient pair in the # sense that for an unrestricted maximum coefficient # from b and A, it has a propagating formula with # a minimum 2-norm of the local truncation error as # T_10,2 ~ .00000034399 # This 2-norm is slightly smaller than T10,2 for the # (16:8,9) pair tabulated above. This 2-norm has two # local maximum values on [0,1]. # # Additional stages and interpolating weights allow # for the computation of an approximation of order up # to order p, at any point of the domain of solution. # These interpolants have continuous derivatives. # # Nodes c[17]=1, c[18] - c[21] were selected to # minimize the maximum of the 2-norm of the local # truncation error over the interval [0,1] for the # interpolant of order 8, and this value is # # Ti_9,2 ~ .0000009164 # # This 2-norm has three local maximum values on [0,1]. # # The remaining five nodes were selected in an attempt # to minimize the maximum of the 2-norm of the local # truncation error on the interval [0,1] for the # interpolant of order 9, and this value is # # Ti_10,2 ~ .000001857 # # This 2-norm has two local maximum values on [0,1]. # # The formulas scanned for this optimal formula are # those developed in J.H. Verner, SIAM J NA 1978, 772-790, # "Explicit Runge--Kutta methods with estimates of the # Local Truncation Error". It is conceivable that the # pairs in J.H. Verner, Applied Numerical Mathematics 22 # 1996, 345--357, "High-order explicit Runge--Kutta pairs # with low stage order", which also require 16 stages, # but solve the order conditions in a different way # may yield particular pairs of equivalant or more # efficiency. Coefficients for the interpolants are # those developed in J.H. Verner, SIAM J NA 1993, # 1446-1466, "Differentiable Interpolants for high order # Runge--Kutta methods". # # Instructions for using the interpolants are contained # in J.H. Verner, SIAM NA 30, 1993, 1446-1466, "Differentiable # Interpolants for high-order Runge--Kutta methods". # ###################################################################### # NODES # ----- c[1] = 0 c[2] = .3571e-1 c[3] = .9906028091267414072062290548990445721420e-1 c[4] = .1485904213690112110809343582348566858213 c[5] = .6134 c[6] = .2327359473605626756631680289030288192566 c[7] = .5538640526394373243368319710969711807434 c[8] = .6555 c[9] = .491625 c[10] = .6858e-1 c[11] = .253 c[12] = .6620641795412045944786226689660879617864 c[13] = .8309 c[14] = .8998 c[15] = 1 c[16] = 1 # # ******************************************************** # COUPLING COEFFICIENTS # --------------------- # for c[1] = 0 # # for c[2] = .3571e-1 a[2,1] = .3571e-1 # # for c[3] = .9906028091267414072062290548990445721420e-1 a[3,1] = -.3833735636677017025757228807792426014327e-1 a[3,2] = .1373976372794443109781951935678287173575 # # for c[4] = .1485904213690112110809343582348566858213 a[4,1] = .3714760534225280277023358955871417145532e-1 a[4,2] = 0 a[4,3] = .1114428160267584083107007686761425143660 # # for c[5] = .6134 a[5,1] = 2.674764429871505119194043347886848033362 a[5,2] = 0 a[5,3] = -9.982382134885293836441318015951354323506 a[5,4] = 7.921017705013788717247274668064506290145 # # for c[6] = .2327359473605626756631680289030288192566 a[6,1] = .5242104050577351069841401193318449756747e-1 a[6,2] = 0 a[6,3] = 0 a[6,4] = .1796911189175953081279466649977364298948 a[6,5] = .6237879371938568368073519721078917943307e-3 # # for c[7] = .5538640526394373243368319710969711807434 a[7,1] = .1592492223647632083060991150129732726703 a[7,2] = 0 a[7,3] = 0 a[7,4] = -.4298429877241087508189185066909803751595 a[7,5] = .6665266542726088051243445942394404211672e-1 a[7,6] = .7578051525715219863372169033510342411158 # # for c[8] = .6555 a[8,1] = .7283333333333333333333333333333333333333e-1 a[8,2] = 0 a[8,3] = 0 a[8,4] = 0 a[8,5] = 0 a[8,6] = .3359344590665103678713422141936031057620 a[8,7] = .2467322076001562987953244524730635609046 # # for c[9] = .491625 a[9,1] = .729755859375e-1 a[9,2] = 0 a[9,3] = 0 a[9,4] = 0 a[9,5] = 0 a[9,6] = .3348009729699333533129043555243262838993 a[9,7] = .1184158239050666466870956444756737161007 a[9,8] = -.345673828125e-1 # # for c[10] = .6858e-1 a[10,1] = .4911213663452096382929799914888692571185e-1 a[10,2] = 0 a[10,3] = 0 a[10,4] = 0 a[10,5] = 0 a[10,6] = .3983857361308652347066808423646954244674e-1 a[10,7] = .1069675288939354812077139476183615284408 a[10,8] = -.2174259165458647598455650055976622663086e-1 a[10,9] = -.1055956474869564925231235304439517699685 # # for c[11] = .253 a[11,1] = -.2707988818641280489175862864774930636734e-1 a[11,2] = 0 a[11,3] = 0 a[11,4] = 0 a[11,5] = 0 a[11,6] = .333e-1 a[11,7] = -.1645526070036057075627040519129748090051 a[11,8] = .3428266306497389946775103799535438271025e-1 a[11,9] = .1585264064439221046855989221358760453749 a[11,10] = .2185234256811225083011127204294936872873 # # for c[12] = .6620641795412045944786226689660879617864 a[12,1] = .5584657769108862727526281006220448941548e-1 a[12,2] = 0 a[12,3] = 0 a[12,4] = 0 a[12,5] = 0 a[12,6] = .9166533166672539087479545790553151183505e-1 a[12,7] = .2392399655523627049569179599100121305027 a[12,8] = .1023834712248414879740316672186561113729e-1 a[12,9] = -.2679331322859542442991625976085644567159e-2 a[12,10] = .4235624181474284646680744174197106205439e-1 a[12,11] = .2253970470166604185504274586005888014087 # # for c[13] = .8309 a[13,1] = -.4802510512725195893103945790000307739284 a[13,2] = 0 a[13,3] = 0 a[13,4] = 0 a[13,5] = 0 a[13,6] = -6.359610162555930097843606206238737368598 a[13,7] = -.2762313898040841385163666588363199062208 a[13,8] = -6.500796633979846747011407331804389650124 a[13,9] = .5734765877040956875980719881692630407295 a[13,10] = 1.347125994868138849716719646646049357523 a[13,11] = 5.936840409706221306933249144503479736029 a[13,12] = 6.590346245333924728433733996560685564589 # # for c[14] = .8998 a[14,1] = .3307533067671401079975082476263855100019 a[14,2] = 0 a[14,3] = 0 a[14,4] = 0 a[14,5] = 0 a[14,6] = 5.956207776829962111958631943398154563177 a[14,7] = -.4868316400481527952968983900652017328476 a[14,8] = 4.462055288206771191616083778225319083232 a[14,9] = .7410258231442071778525406966672954745201 a[14,10] = -.7118192034575913119318346652024763899638 a[14,11] = -5.454619594516665440649148499095762263275 a[14,12] = -4.140803729244709693305868672684077643030 a[14,13] = .2038319723190386517589855611303633981843 # # for c[15] = 1 a[15,1] = -.5847111122998944389651719860007729221315 a[15,2] = 0 a[15,3] = 0 a[15,4] = 0 a[15,5] = 0 a[15,6] = -12.41268417116267068926060969930886386670 a[15,7] = 1.360245445660928146977003413153919520189 a[15,8] = -22.42610531111868229005965344123311160860 a[15,9] = -.8828857055865458252669099251378431470909 a[15,10] = 1.770155128538230466813900572142546329490 a[15,11] = 12.15809651918533877218739262890752294243 a[15,12] = 22.23037520407760703834156132941550000589 a[15,13] = -.6634483760201249006317401316871863808959 a[15,14] = .4509623787258137198642272397482891274306 # # for c[16] = 1 a[16,1] = 1.940575549810648717395482898478643635949 a[16,2] = 0 a[16,3] = 0 a[16,4] = 0 a[16,5] = 0 a[16,6] = 21.97798408114556310164832153356570429956 a[16,7] = .8230747326984728145128172720328075817426 a[16,8] = 68.16441683626354817206843139989369741474 a[16,9] = -3.117097463620266666801525529574905518624 a[16,10] = -4.568841021822439620303730003755793384193 a[16,11] = -18.74190987126264955904919353670543412590 a[16,12] = -66.57711839637831878350036856645874492472 a[16,13] = 1.098915553165441824029764532524025021446 a[16,14] = 0 a[16,15] = 0 # # ******************************************************** # High order weights c[17] = 1 # i.e. This is the propagating stage, and stage 17 as well. # ------------------------------------------------------- # b[1] = .15006690149797247957662887123770409800223878439125635509305e-1 b[2] = 0 b[3] = 0 b[4] = 0 b[5] = 0 b[6] = 0 b[7] = 0 b[8] = -1.05518099274638127859438168508018447470958212998793507706900 b[9] = .238494726378218311263814055777485195395241336180778947138354 b[10] = .128815177428299135462251588714408117237111752488648622524665 b[11] = .227662311104621561461491789666574650857365714720222465635242 b[12] = 1.22953258743751744316032181517357642189586257623541890214594 b[13] = .462497666281038348739730827525612863770345382666305044524257e-1 b[14] = .138619631936629390362830862120184349246824228027665359337989 b[15] = .308001016831943540520356037516240438999181056294446403250787e-1 b[16] = 0 # # ******************************************************** # Low order weights c[extra] = 1 # ------------------------------------------------------- # bh[1] = .1897210532481101330735875918987801428423e-1 bh[2] = 0 bh[3] = 0 bh[4] = 0 bh[5] = 0 bh[6] = 0 bh[7] = 0 bh[8] = 3.408110314549493848404398964228060776034 bh[9] = .1260323883820920906560270507988839661845 bh[10] = .1188375063451149770930378540709095182997 bh[11] = .2491041997838687569190073177760326138188 bh[12] = -3.269966219928978218713853055139116510335 bh[13] = .3023798100228882907409723963501699715077 bh[14] = 0 bh[15] = 0 bh[16] = .4652989552070924159305071272518165020637e-1 # #*******************************************************` # # Largest coefficient in b or A has magnitude `, 68.16442 # #*******************************************************` # SUMMARY OF NORMS OF ERRORS: A101, A102, A10inf` # ----------------------------------------------------` # A_[10, 1] = `, .6298637441e-5 # A_[10, 2] = `, .3439902103e-6 # A_[10,oo] = `, .5849302357e-7 #**************************************************** # # END OF GENERATION OF A PAIR OF RK METHODS` # #############################################################` # ` # START OF GENERATION OF STABILITY INTERVALS ` # ` #############################################################` # # Stability Boundaries of High Order Method` # -----------------------------------------` # Real Stability Interval is nearly [`, -4.431330766, `,0]` # # Stability Boundaries of Low Order Method` # ----------------------------------------` # Real Stability Interval is nearly [`, -3.828196857, `,0]` # ############################################################ # # START OF GENERATION OF INTERPOLANT # # ******************************************************** # # FIVE ADDITIONAL STAGES FOR INTERPOLANT OF ORDER 8 # # Coupling coefficients for j from 1 to 16 # for node c[17] = 1 # ---------------------------------------------------- a[17,1] = .1500669014979724795766288712377040980022e-1 a[17,2] = 0 a[17,3] = 0 a[17,4] = 0 a[17,5] = 0 a[17,6] = 0 a[17,7] = 0 a[17,8] = -1.055180992746381278594381685080184474710 a[17,9] = .2384947263782183112638140557774851953952 a[17,10] = .1288151774282991354622515887144081172371 a[17,11] = .2276623111046215614614917896665746508574 a[17,12] = 1.229532587437517443160321815173576421896 a[17,13] = .4624976662810383487397308275256128637703e-1 a[17,14] = .1386196319366293903628308621201843492468 a[17,15] = .3080010168319435405203560375162404389992e-1 a[17,16] = 0 # # ******************************************************** # # Coupling coefficients for j from 1 to 17 # for node c[18] = .7375018139988810429214526568476200466933 # ---------------------------------------------------- a[18,1] = .1571801061417788021580499748292804743973e-1 a[18,2] = 0 a[18,3] = 0 a[18,4] = 0 a[18,5] = 0 a[18,6] = 0 a[18,7] = 0 a[18,8] = .4853403452657363080900487157345670578815 a[18,9] = .2107787568904546684574678306538806242458 a[18,10] = .1269802413053354152906551816812865612400 a[18,11] = .2319687014513919157472314236852033704833 a[18,12] = -.3620214714069096774945761015915844860748 a[18,13] = .5366106712036344076786687672191174826870e-1 a[18,14] = -.2806066613385549629370667652641771053798e-1 a[18,15] = -.2378121372710330174923994809623697280793e-1 a[18,16] = 0 a[18,17] = .2691804261928988988990035710208180655494e-1 # # ******************************************************** # # Coupling coefficients for j from 1 to 18 # for node c[19] = .749 # ---------------------------------------------------- a[19,1] = .1569705832522204237038521016179271074328e-1 a[19,2] = 0 a[19,3] = 0 a[19,4] = 0 a[19,5] = 0 a[19,6] = 0 a[19,7] = 0 a[19,8] = .4616075242202111936044636619654233606961 a[19,9] = .2113946516698113234089715994366982070261 a[19,10] = .1270330917167109232354671760560477406684 a[19,11] = .2318540550298708286361100159494550129997 a[19,12] = -.3385266406688372942274816228586594758312 a[19,13] = .5298251972194235542699858758068190007453e-1 a[19,14] = -.2750461365887187780272257875389040879148e-1 a[19,15] = -.2361906185395527177468217904834875318962e-1 a[19,16] = 0 a[19,17] = .2668458089504036186694391521145359328377e-1 a[19,18] = .1139683460285541525554621429934611232038e-1 # # ******************************************************** # # Coupling coefficients for j from 1 to 19 # for node c[20] = .65 # ---------------------------------------------------- a[20,1] = .1438964884291216402951708307100553583630e-1 a[20,2] = 0 a[20,3] = 0 a[20,4] = 0 a[20,5] = 0 a[20,6] = 0 a[20,7] = 0 a[20,8] = -1.206901219123788519767229935392578222389 a[20,9] = .2505628554639376294459724790398627595052 a[20,10] = .1303332915702127271150218433265955255131 a[20,11] = .2246717759263522068724895908997148014087 a[20,12] = 1.308419325781946480763221457213122900404 a[20,13] = .2589750180376235922034567811002287469857e-2 a[20,14] = .8070743254562856892787425005539711000810e-2 a[20,15] = -.1267568255392829384757055330207124946149e-1 a[20,16] = 0 a[20,17] = .1129158072373321675992773077865249355325e-1 a[20,18] = .3422056680709749787491027385482635813798e-1 a[20,19] = -.1149726368734142020610819623056729009783 # # ******************************************************** # # Coupling coefficients for j from 1 to 20 # for node c[21] = .487 # ---------------------------------------------------- a[21,1] = .1452348029801041972968851600082942318526e-1 a[21,2] = 0 a[21,3] = 0 a[21,4] = 0 a[21,5] = 0 a[21,6] = 0 a[21,7] = 0 a[21,8] = -.5214243102465818231380106106018875132310 a[21,9] = .1866698844204603947548300448071067931901 a[21,10] = .1299316354451273004594203171440036506875 a[21,11] = .2262141085765719298250459612902177827497 a[21,12] = .6104582639466710747030070634860296980845 a[21,13] = .1418715607022412453299982942636029187870e-1 a[21,14] = .1480061054412245756808618945399364705863e-1 a[21,15] = -.3711471609871774475268028540581604884370e-2 a[21,16] = 0 a[21,17] = .1393256979572558890317889514899316376009e-2 a[21,18] = 1.114731020624073270253135239909350910479 a[21,19] = -1.021208555757145859230933049243060207745 a[21,20] = -.2795650792912340738723193626472621878286 # # -------------------------------------------------------- # COEFFICIENTS FOR INTERPOLANT bi8 WITH 21 STAGES # for i from 1 to 21, bi8[i] = SUM_j=1^8 bi8[i,j]u^j # -------------------------------------------------------- # # COEFFICIENTS OF bi8[1] bi8[1,1] = 1 u bi8[1,2] = -11.009172273099548270487991762298741443 u^2 bi8[1,3] = 53.6037684418544552607732109075661671085 u^3 bi8[1,4] = -143.78186010517733940062488419099653948 u^4 bi8[1,5] = 227.780227272405967284751656763039821064 u^5 bi8[1,6] = -212.76677747453975152813132479196141561 u^6 bi8[1,7] = 108.476558503094350107230415191741302374 u^7 bi8[1,8] = -23.287737674388336205553419229966823592 u^8 # ----------------------------------------- # # COEFFICIENTS OF bi8[2] bi8[2,1] = 0 u bi8[2,2] = 0 u^2 bi8[2,3] = 0 u^3 bi8[2,4] = 0 u^4 bi8[2,5] = 0 u^5 bi8[2,6] = 0 u^6 bi8[2,7] = 0 u^7 bi8[2,8] = 0 u^8 # ----------------------------------------- # # COEFFICIENTS OF bi8[3] bi8[3,1] = 0 u bi8[3,2] = 0 u^2 bi8[3,3] = 0 u^3 bi8[3,4] = 0 u^4 bi8[3,5] = 0 u^5 bi8[3,6] = 0 u^6 bi8[3,7] = 0 u^7 bi8[3,8] = 0 u^8 # ----------------------------------------- # # COEFFICIENTS OF bi8[4] bi8[4,1] = 0 u bi8[4,2] = 0 u^2 bi8[4,3] = 0 u^3 bi8[4,4] = 0 u^4 bi8[4,5] = 0 u^5 bi8[4,6] = 0 u^6 bi8[4,7] = 0 u^7 bi8[4,8] = 0 u^8 # ----------------------------------------- # # COEFFICIENTS OF bi8[5] bi8[5,1] = 0 u bi8[5,2] = 0 u^2 bi8[5,3] = 0 u^3 bi8[5,4] = 0 u^4 bi8[5,5] = 0 u^5 bi8[5,6] = 0 u^6 bi8[5,7] = 0 u^7 bi8[5,8] = 0 u^8 # ----------------------------------------- # # COEFFICIENTS OF bi8[6] bi8[6,1] = 0 u bi8[6,2] = 0 u^2 bi8[6,3] = 0 u^3 bi8[6,4] = 0 u^4 bi8[6,5] = 0 u^5 bi8[6,6] = 0 u^6 bi8[6,7] = 0 u^7 bi8[6,8] = 0 u^8 # ----------------------------------------- # # COEFFICIENTS OF bi8[7] bi8[7,1] = 0 u bi8[7,2] = 0 u^2 bi8[7,3] = 0 u^3 bi8[7,4] = 0 u^4 bi8[7,5] = 0 u^5 bi8[7,6] = 0 u^6 bi8[7,7] = 0 u^7 bi8[7,8] = 0 u^8 # ----------------------------------------- # # COEFFICIENTS OF bi8[8] bi8[8,1] = 0 u bi8[8,2] = 179.087906217887237219763763517805603405 u^2 bi8[8,3] = -1728.4031360576313422200193755550186656 u^3 bi8[8,4] = 6564.46094056114241668520297149455937609 u^4 bi8[8,5] = -12902.910815841211056904979095488699523 u^5 bi8[8,6] = 13965.1797710026505924517219455017192980 u^6 bi8[8,7] = -7920.4240616873754433765592473279592986 u^7 bi8[8,8] = 1841.95421481179121486627465617251302548 u^8 # ----------------------------------------- # # COEFFICIENTS OF bi8[9] bi8[9,1] = 0 u bi8[9,2] = -22.883047147673905485214697445226699950 u^2 bi8[9,3] = 238.804050768529415536649476897078810139 u^3 bi8[9,4] = -936.85857822960737524933031174532417444 u^4 bi8[9,5] = 1874.04495776215135386789321723467192169 u^5 bi8[9,6] = -2050.5115593530093548177138133074309935 u^6 bi8[9,7] = 1171.50854519241625074722053452432895039 u^7 bi8[9,8] = -273.86587426642816628824059210232032907 u^8 # ----------------------------------------- # # COEFFICIENTS OF bi8[10] bi8[10,1] = 0 u bi8[10,2] = 12.1733154962812480100234994056280809726 u^2 bi8[10,3] = -82.750764764368702721462010574562168251 u^3 bi8[10,4] = 256.480443875130775557432334036545707434 u^4 bi8[10,5] = -441.09994297583109625949510679373047644 u^5 bi8[10,6] = 434.507983140121315864244979200535231052 u^6 bi8[10,7] = -229.89346058969002654886690501730424215 u^7 bi8[10,8] = 50.7112409957847852335854613316022755067 u^8 # ----------------------------------------- # # COEFFICIENTS OF bi8[11] bi8[11,1] = 0 u bi8[11,2] = 2.61320352013362275655103290153819291495 u^2 bi8[11,3] = 16.8800293726118444769201726246619058295 u^3 bi8[11,4] = -134.17279211552377899148845724211730771 u^4 bi8[11,5] = 340.120480169263766446505999236215647471 u^5 bi8[11,6] = -420.12588937768794008574360874572298723 u^6 bi8[11,7] = 258.323437214655710365707922536062173492 u^7 bi8[11,8] = -63.410806472348603406991569520971050106 u^8 # ----------------------------------------- # # COEFFICIENTS OF bi8[12] bi8[12,1] = 0 u bi8[12,2] = -212.31270059136659921947967158064013199 u^2 bi8[12,3] = 2045.35231639515436471224930279882472314 u^3 bi8[12,4] = -7762.0604136440909067904491207875065305 u^4 bi8[12,5] = 15250.1484145937813254342184379271259724 u^5 bi8[12,6] = -16501.080035445178597970138849040394763 u^6 bi8[12,7] = 9356.90736395730475810033295746054352009 u^7 bi8[12,8] = -2175.7254126781668268235727349627792128 u^8 # ----------------------------------------- # # COEFFICIENTS OF bi8[13] bi8[13,1] = 0 u bi8[13,2] = -11.501647085674297642838265513207120229 u^2 bi8[13,3] = 107.276970404572523556452634967780744097 u^3 bi8[13,4] = -401.23453329839476276779635630884895370 u^4 bi8[13,5] = 781.891698957293030032303464174078218081 u^5 bi8[13,6] = -841.65933710132544301884545345304835330 u^6 bi8[13,7] = 475.576739148558845765724142521617358352 u^7 bi8[13,8] = -110.30364125840179209012619330561933200 u^8 # ----------------------------------------- # # COEFFICIENTS OF bi8[14] bi8[14,1] = 0 u bi8[14,2] = -38.772396773403658004983716330925504164 u^2 bi8[14,3] = 358.639118015297013691293498291937710953 u^3 bi8[14,4] = -1336.2157321402269203215700860140131895 u^4 bi8[14,5] = 2598.19351531318703415004712031537458940 u^5 bi8[14,6] = -2792.8783300510167773490832509904683870 u^6 bi8[14,7] = 1576.58679034281004823064293851369170916 u^7 bi8[14,8] = -365.41434507471011100598367292347674440 u^8 # ----------------------------------------- # # COEFFICIENTS OF bi8[15] bi8[15,1] = 0 u bi8[15,2] = -10.004254699988526742118869843022452495 u^2 bi8[15,3] = 91.6775424156049820846038288550384217432 u^3 bi8[15,4] = -340.07776006528482087050428112358069131 u^4 bi8[15,5] = 659.600836997697836608703753625575467028 u^5 bi8[15,6] = -707.88262841129756715154773051891381715 u^6 bi8[15,7] = 399.158003026012712821111684268379597103 u^7 bi8[15,8] = -92.440939161061422396196349659724900862 u^8 # ----------------------------------------- # # COEFFICIENTS OF bi8[16] bi8[16,1] = 0 u bi8[16,2] = 0 u^2 bi8[16,3] = 0 u^3 bi8[16,4] = 0 u^4 bi8[16,5] = 0 u^5 bi8[16,6] = 0 u^6 bi8[16,7] = 0 u^7 bi8[16,8] = 0 u^8 # ----------------------------------------- # # COEFFICIENTS OF bi8[17] bi8[17,1] = 0 u bi8[17,2] = 16.2101449133756113052556293139075916278 u^2 bi8[17,3] = -150.44921310486628831227271221175946775 u^3 bi8[17,4] = 565.119663172215249030697931716027311294 u^4 bi8[17,5] = -1111.2102826515537053219989040609684290 u^5 bi8[17,6] = 1211.23697600084352654995897551778506557 u^6 bi8[17,7] = -695.33656069180915952688318054142230013 u^7 bi8[17,8] = 164.429272361794766275242260266430228456 u^8 # ----------------------------------------- # # COEFFICIENTS OF bi8[18] bi8[18,1] = 0 u bi8[18,2] = -1872.8050193184148301080651477238531464 u^2 bi8[18,3] = 16400.6663320704863438658761751771242437 u^3 bi8[18,4] = -57270.162601660532870808173776360907291 u^4 bi8[18,5] = 103330.591296066351625435679442988921099 u^5 bi8[18,6] = -102051.18881737274653810677122429441673 u^6 bi8[18,7] = 52424.7526087466269444582092355231970059 u^7 bi8[18,8] = -10961.853798531770674736754705310065172 u^8 # ----------------------------------------- # # COEFFICIENTS OF bi8[19] bi8[19,1] = 0 u bi8[19,2] = 1720.56212905796882726815179456740888746 u^2 bi8[19,3] = -15159.620464322547473818641003813529098 u^3 bi8[19,4] = 53314.5956193552349291147268699196667133 u^4 bi8[19,5] = -96988.374686915763514264143535872202761 u^5 bi8[19,6] = 96690.3918102368707328017014251049574348 u^6 bi8[19,7] = -50199.312489882466233785585470608781270 u^7 bi8[19,8] = 10621.7580824707027326837899207024800951 u^8 # ----------------------------------------- # # COEFFICIENTS OF bi8[20] bi8[20,1] = 0 u bi8[20,2] = 264.247048839608472575101341495756581748 u^2 bi8[20,3] = -2232.3380039056847315097415410801346265 u^3 bi8[20,4] = 7453.75380241914399865580881700002668831 u^4 bi8[20,5] = -12770.718814009167305619604409871514599 u^5 bi8[20,6] = 11903.2522676744622422517513883009525203 u^6 bi8[20,7] = -5733.1555765072257401698251585613343530 u^7 bi8[20,8] = 1114.95927548886306381650956271624778871 u^8 # ----------------------------------------- # # COEFFICIENTS OF bi8[21] bi8[21,1] = 0 u bi8[21,2] = -15.605510155633653661658701002871141419 u^2 bi8[21,3] = 40.6614542709875953973183427149912998699 u^3 bi8[21,4] = 170.153801875971406156068349606468881911 u^4 bi8[21,5] = -848.05688473860526088988204017788694648 u^5 bi8[21,6] = 1373.52456653185356010859654151640790661 u^6 bi8[21,7] = -993.16789677291301718845986848276015224 u^7 bi8[21,8] = 272.490468988339370078017375825650151758 u^8 # #******************************************************** # # FIVE ADDITIONAL STAGES FOR INTERPOLANT OF ORDER 9 # # Coupling coefficients for j from 1 to 21 # for node c[22] = .97e-2 # ---------------------------------------------------- a[22,1] = .8711816186418632397551552134391273160887e-2 a[22,2] = 0 a[22,3] = 0 a[22,4] = 0 a[22,5] = 0 a[22,6] = 0 a[22,7] = 0 a[22,8] = .1532993247326560354009297246756198757646e-1 a[22,9] = -.1943250606288015181452523235039064995331e-2 a[22,10] = .1072095950570478409673351209542485874475e-2 a[22,11] = .2601233036074381124605301648955068308034e-3 a[22,12] = -.1817718521410219268803888649170562657549e-1 a[22,13] = -.9877668338996712912794853387601975290298e-3 a[22,14] = -.3332383192417756458522769371199819990865e-2 a[22,15] = -.8605833352714280922890371244764817761631e-3 a[22,16] = 0 a[22,17] = .1392810143886650671629507962655268250424e-2 a[22,18] = -.1617419987897630058856682050471490927435 a[22,19] = .1485156561306044255987175208430691701911 a[22,20] = .2289051095253062003982633508518682297336e-1 a[22,21] = -.1429777169141779172700863258972231247118e-2 # # ******************************************************** # # Coupling coefficients for j from 1 to 22 # for node c[23] = .138 # ---------------------------------------------------- a[23,1] = .2710092628714179537584868171073201523422e-1 a[23,2] = 0 a[23,3] = 0 a[23,4] = 0 a[23,5] = 0 a[23,6] = 0 a[23,7] = 0 a[23,8] = .6923104986003885114772983986172611691902 a[23,9] = -.6725243132164494445073150422370697609664e-1 a[23,10] = .8808440659269460280192719417821072306037e-1 a[23,11] = .5982566312199630889836790781586188392095e-1 a[23,12] = -.8251290314815523858109567180972702723517 a[23,13] = -.4886453387508305661509501250020800230341e-1 a[23,14] = -.1682608371566511953625070499451460275789 a[23,15] = -.4443170503743608176862370639456174098455e-1 a[23,16] = 0 a[23,17] = .7037844639432780238885500484343585742713e-1 a[23,18] = -8.818686397504859095347055843377871115177 a[23,19] = 8.028821919603921347277176943010114087042 a[23,20] = 1.306591406499157862229641342220359676579 a[23,21] = -.1624883307224014710941456378572112779619 a[23,22] = 0 # # ******************************************************** # # Coupling coefficients for j from 1 to 23 # for node c[24] = .249 # ---------------------------------------------------- a[24,1] = .1466442426857961230207869561545353695312e-1 a[24,2] = 0 a[24,3] = 0 a[24,4] = 0 a[24,5] = 0 a[24,6] = 0 a[24,7] = 0 a[24,8] = .1898399561399124450873282966092465494738 a[24,9] = .3711530651907537053995301753121055816571e-1 a[24,10] = .1316247981199510429999180335475037001729 a[24,11] = .1466622928822685713088228930531811168096 a[24,12] = -.2377331162467838377764310139745632904067 a[24,13] = -.2493071232192949761907468241344412740688e-1 a[24,14] = -.9427777049221307692559326299260474943801e-1 a[24,15] = -.2726667627690941070300686079610559177653e-1 a[24,16] = 0 a[24,17] = .4097659404063362064133671142526181935526e-1 a[24,18] = -5.538725321977494545432278985839620727708 a[24,19] = 4.973442175247713118429917494679987165227 a[24,20] = .8626294449189858956229369969895335684221 a[24,21] = -.2250213948217893084759073334350395278429 a[24,22] = 0 a[24,23] = 0 # # ******************************************************** # # Coupling coefficients for j from 1 to 24 # for node c[25] = .439 # ---------------------------------------------------- a[25,1] = .1191252689920919922350166270042464848683e-1 a[25,2] = 0 a[25,3] = 0 a[25,4] = 0 a[25,5] = 0 a[25,6] = 0 a[25,7] = 0 a[25,8] = -.6723137858130876080443132301001699876552 a[25,9] = .1801284266825688933352476317750604591900 a[25,10] = .1365689447449347163335592887842694005196 a[25,11] = .2111038393798906615990851223053694372026 a[25,12] = .7775854726244728792738759030556073621755 a[25,13] = .2362117948505939633346199470372118351743e-1 a[25,14] = .6391325607581230059121991748062989083367e-1 a[25,15] = .1197654289246164466950805226357270574198e-1 a[25,16] = 0 a[25,17] = -.2188622145087031801816344581960035242509e-1 a[25,18] = 3.500218060218659348421859254969160434525 a[25,19] = -3.195765244251373581009030415350620096609 a[25,20] = -.5994401048600148264885439684062323011428 a[25,21] = .1137710737227729377873223163880721563872e-1 a[25,22] = 0 a[25,23] = 0 a[25,24] = 0 # # ******************************************************** # # Coupling coefficients for j from 1 to 25 # for node c[26] = .794 # ---------------------------------------------------- a[26,1] = .1184012014074604271397800922251586670558e-1 a[26,2] = 0 a[26,3] = 0 a[26,4] = 0 a[26,5] = 0 a[26,6] = 0 a[26,7] = 0 a[26,8] = -.6641263289768790103963252230448945748496 a[26,9] = .1788914061873194018510863179806909040625 a[26,10] = .1367562642287085016908037148859544845819 a[26,11] = .2107758119730286112356534765669813731420 a[26,12] = .7679184744806335059253231677825531412084 a[26,13] = .2313496869580297355994450248598930620113e-1 a[26,14] = .6230605705999706953753333719939087802061e-1 a[26,15] = .1157099029193866324646273674695267419004e-1 a[26,16] = 0 a[26,17] = -.2129441601042184344575923744216391881547e-1 a[26,18] = 3.163166251075345445930896810245694869610 a[26,19] = -2.819544872276817190650813744596266559171 a[26,20] = -.3978769046309758754403585402077121882020 a[26,21] = .1304821777615737042415746721743137433164 a[26,22] = 0 a[26,23] = 0 a[26,24] = 0 a[26,25] = 0 # # -------------------------------------------------------- # COEFFICIENTS FOR INTERPOLANT bi9 WITH 26 STAGES # for i from 1 to 26, bi9[i] = SUM_j=1^9 bi9[i,j]u^j # ----------------------------------------- # # COEFFICIENTS OF bi9[1] bi9[1,1] = 1 u bi9[1,2] = -60.671564990962758251925513636023967612 u^2 bi9[1,3] = 669.417333989096529852883401874298410509 u^3 bi9[1,4] = -3377.8789462251993515202266905133322562 u^4 bi9[1,5] = 9286.46896739104669222569524723782154199 u^5 bi9[1,6] = -14780.477136811020541163339406234475962 u^6 bi9[1,7] = 13604.9938632822472010264750780633453117 u^7 bi9[1,8] = -6724.7054433561144859867823236173963799 u^8 bi9[1,9] = 1381.86793341105651106517786971288707205 u^9 # ----------------------------------------- # # COEFFICIENTS OF bi9[2] bi9[2,1] = 0 u bi9[2,2] = 0 u^2 bi9[2,3] = 0 u^3 bi9[2,4] = 0 u^4 bi9[2,5] = 0 u^5 bi9[2,6] = 0 u^6 bi9[2,7] = 0 u^7 bi9[2,8] = 0 u^8 bi9[2,9] = 0 u^9 # ----------------------------------------- # # COEFFICIENTS OF bi9[3] bi9[3,1] = 0 u bi9[3,2] = 0 u^2 bi9[3,3] = 0 u^3 bi9[3,4] = 0 u^4 bi9[3,5] = 0 u^5 bi9[3,6] = 0 u^6 bi9[3,7] = 0 u^7 bi9[3,8] = 0 u^8 bi9[3,9] = 0 u^9 # ----------------------------------------- # # COEFFICIENTS OF bi9[4] bi9[4,1] = 0 u bi9[4,2] = 0 u^2 bi9[4,3] = 0 u^3 bi9[4,4] = 0 u^4 bi9[4,5] = 0 u^5 bi9[4,6] = 0 u^6 bi9[4,7] = 0 u^7 bi9[4,8] = 0 u^8 bi9[4,9] = 0 u^9 # ----------------------------------------- # # COEFFICIENTS OF bi9[5] bi9[5,1] = 0 u bi9[5,2] = 0 u^2 bi9[5,3] = 0 u^3 bi9[5,4] = 0 u^4 bi9[5,5] = 0 u^5 bi9[5,6] = 0 u^6 bi9[5,7] = 0 u^7 bi9[5,8] = 0 u^8 bi9[5,9] = 0 u^9 # ----------------------------------------- # # COEFFICIENTS OF bi9[6] bi9[6,1] = 0 u bi9[6,2] = 0 u^2 bi9[6,3] = 0 u^3 bi9[6,4] = 0 u^4 bi9[6,5] = 0 u^5 bi9[6,6] = 0 u^6 bi9[6,7] = 0 u^7 bi9[6,8] = 0 u^8 bi9[6,9] = 0 u^9 # ----------------------------------------- # # COEFFICIENTS OF bi9[7] bi9[7,1] = 0 u bi9[7,2] = 0 u^2 bi9[7,3] = 0 u^3 bi9[7,4] = 0 u^4 bi9[7,5] = 0 u^5 bi9[7,6] = 0 u^6 bi9[7,7] = 0 u^7 bi9[7,8] = 0 u^8 bi9[7,9] = 0 u^9 # ----------------------------------------- # # COEFFICIENTS OF bi9[8] bi9[8,1] = 0 u bi9[8,2] = -.59721676056684470801952994611021858852 u^2 bi9[8,3] = 47.9360483641316733888543826427761136649 u^3 bi9[8,4] = -565.71704170023321536474066239412313218 u^4 bi9[8,5] = 2799.62270720285856115169650129725180636 u^5 bi9[8,6] = -7032.1617949196787314919454503249805212 u^6 bi9[8,7] = 9404.15458349876626141738039605732324808 u^7 bi9[8,8] = -6374.6618043443040550263461557577033090 u^8 bi9[8,9] = 1720.36933766627996935452613674048582851 u^9 # ----------------------------------------- # # COEFFICIENTS OF bi9[9] bi9[9,1] = 0 u bi9[9,2] = .134984470796006971722056686682629898492 u^2 bi9[9,3] = -10.834629145944526812223590265084456754 u^3 bi9[9,4] = 127.864823186993429826665223280937285860 u^4 bi9[9,5] = -632.77793677721873323380079062006984876 u^5 bi9[9,6] = 1589.42732541225546480426354495668561928 u^6 bi9[9,7] = -2125.5512463055572075298710635247357933 u^7 bi9[9,8] = 1440.81748366575488132149000954073339262 u^8 bi9[9,9] = -388.84230978070109703698157599937134358 u^9 # ----------------------------------------- # # COEFFICIENTS OF bi9[10] bi9[10,1] = 0 u bi9[10,2] = .729074760675326573955388843418390629860 u^2 bi9[10,3] = -5.8519729010332163134273909972861513842 u^3 bi9[10,4] = 69.0620297387629935587576143813403566689 u^4 bi9[10,5] = -341.77444271620971355669272484345138739 u^5 bi9[10,6] = 858.477527120137300045841924638667044538 u^6 bi9[10,7] = -1148.0474435798640829627929492591018069 u^7 bi9[10,8] = 778.210749640921413718998790161059954495 u^8 bi9[10,9] = -210.02053960135392801261855137685544089 u^9 # ----------------------------------------- # # COEFFICIENTS OF bi9[11] bi9[11,1] = 0 u bi9[11,2] = .128853484734578556111620332930181436020 u^2 bi9[11,3] = -10.342520980591798291955500435733899743 u^3 bi9[11,4] = 122.057211066254163569747789887394785533 u^4 bi9[11,5] = -604.03720321370038943548927429014941953 u^5 bi9[11,6] = 1517.23563758112275748480130150120874085 u^6 bi9[11,7] = -2029.0088441529040679457612821944539176 u^7 bi9[11,8] = 1375.37564537632093724244189399979666341 u^8 bi9[11,9] = -371.18111685013155961843505701132655968 u^9 # ----------------------------------------- # # COEFFICIENTS OF bi9[12] bi9[12,1] = 0 u bi9[12,2] = .695897172076238736216479142518749276746 u^2 bi9[12,3] = -55.856705135748306493974842176201132618 u^3 bi9[12,4] = 659.192634079572821155723569480922153284 u^4 bi9[12,5] = -3262.2150841408463540688882313635567019 u^5 bi9[12,6] = 8194.11280758829012536721337359229150737 u^6 bi9[12,7] = -10958.039044672965515394222094269717056 u^7 bi9[12,8] = 7427.97157664347578998251574120831336997 u^8 bi9[12,9] = -2004.6325489464172818414236737993973128 u^9 # ----------------------------------------- # # COEFFICIENTS OF bi9[13] bi9[13,1] = 0 u bi9[13,2] = .261766805813262302830045074618681553995 u^2 bi9[13,3] = -2.1010907750946045040966004715281676569 u^3 bi9[13,4] = 24.7960125666003203913108259157166833632 u^4 bi9[13,5] = -122.71060391057876840039391579311663034 u^5 bi9[13,6] = 308.227540243681394314664842171282725434 u^6 bi9[13,7] = -412.19464510006695825973474777273784211 u^7 bi9[13,8] = 279.408578064554483146960831647219592946 u^8 bi9[13,9] = -75.405718003049089084120267121545668488 u^9 # ----------------------------------------- # # COEFFICIENTS OF bi9[14] bi9[14,1] = 0 u bi9[14,2] = .784566516126207461850266507371358159569 u^2 bi9[14,3] = -6.2973816116961953548351869709538775981 u^3 bi9[14,4] = 74.3185184720378264480315012555342585254 u^4 bi9[14,5] = -367.78777470565001975434105625062192278 u^5 bi9[14,6] = 923.818459990863967296803792781071063324 u^6 bi9[14,7] = -1235.4282876597389441773173089892267661 u^7 bi9[14,8] = 837.442371605615424660411371389701714778 u^8 bi9[14,9] = -226.00574311110805047457530900412142160 u^9 # ----------------------------------------- # # COEFFICIENTS OF bi9[15] bi9[15,1] = 0 u bi9[15,2] = .174323998241199708354266609961877656592 u^2 bi9[15,3] = -1.3992245634210807475661124260479057978 u^3 bi9[15,4] = 16.5129418820673116921688916038940504367 u^4 bi9[15,5] = -81.719311330651921118221550201203029441 u^5 bi9[15,6] = 205.264594249812809915291505118941991081 u^6 bi9[15,7] = -274.50164417988052816131283790668820527 u^7 bi9[15,8] = 186.072563019500146464809911604971039714 u^8 bi9[15,9] = -50.216551375567663661953198851112504438 u^9 # ----------------------------------------- # # COEFFICIENTS OF bi9[16] bi9[16,1] = 0 u bi9[16,2] = 0 u^2 bi9[16,3] = 0 u^3 bi9[16,4] = 0 u^4 bi9[16,5] = 0 u^5 bi9[16,6] = 0 u^6 bi9[16,7] = 0 u^7 bi9[16,8] = 0 u^8 bi9[16,9] = 0 u^9 # ----------------------------------------- # # COEFFICIENTS OF bi9[17] bi9[17,1] = 0 u bi9[17,2] = -.31202370347393561021490909104589729915 u^2 bi9[17,3] = 2.50532351385866542582114814692502855668 u^3 bi9[17,4] = -29.632099441822688910941870509585694850 u^4 bi9[17,5] = 147.173574899191829926941842865950673508 u^5 bi9[17,6] = -371.66826262970959497366149113587895821 u^6 bi9[17,7] = 500.803253762596738399023907244798218745 u^7 bi9[17,8] = -342.94904651443880465789781251565182585 u^8 bi9[17,9] = 93.7984587806712483517357668125471478402 u^9 # ----------------------------------------- # # COEFFICIENTS OF bi9[18] bi9[18,1] = 0 u bi9[18,2] = 1.00518918460055756612940887077202711807 u^2 bi9[18,3] = -80.145247628498561771052351691392668350 u^3 bi9[18,4] = 904.286287161854048273552227532129961260 u^4 bi9[18,5] = -4159.3779007972609160773661441354362080 u^5 bi9[18,6] = 9382.94844063873906978909606923086757913 u^6 bi9[18,7] = -10865.074362321895351322875097650589448 u^7 bi9[18,8] = 6171.21873158613438388357367454262300521 u^8 bi9[18,9] = -1354.8611378236732303410577866989742481 u^9 # ----------------------------------------- # # COEFFICIENTS OF bi9[19] bi9[19,1] = 0 u bi9[19,2] = -1.0187903048921490083392086220859123412 u^2 bi9[19,3] = 81.2296854343308670373230627559449087708 u^3 bi9[19,4] = -916.52209984084987508799139721191838178 u^4 bi9[19,5] = 4215.65805187092305730244216181520585548 u^5 bi9[19,6] = -9509.9082332688606953961331829842253224 u^6 bi9[19,7] = 11012.0886613642691804346667578915363039 u^7 bi9[19,8] = -6254.7209116731392876163095442126962197 u^8 bi9[19,9] = 1373.19363641821890233434135056823876813 u^9 # ----------------------------------------- # # COEFFICIENTS OF bi9[20] bi9[20,1] = 0 u bi9[20,2] = .148675638316259188843832285146684434829 u^2 bi9[20,3] = -1.1854132567011975185609013630414468662 u^3 bi9[20,4] = 13.3751280877394880365070372626692982663 u^4 bi9[20,5] = -61.520574820481602947755579141965921339 u^5 bi9[20,6] = 138.781422449831194026442067802994967214 u^6 bi9[20,7] = -160.70326769519984450233427661090825361 u^7 bi9[20,8] = 91.2773334775208822130042056459948017083 u^8 bi9[20,9] = -20.039495806540545226186936824258113811 u^9 # ----------------------------------------- # # COEFFICIENTS OF bi9[21] bi9[21,1] = 0 u bi9[21,2] = -.26119631066917488492992560550370701830 u^2 bi9[21,3] = .208255752438779775536051016631482507430 u^3 bi9[21,4] = -2.3497690346644753101563642562153051718 u^4 bi9[21,5] = 10.8080566226830373409785310179170015237 u^5 bi9[21,6] = -24.381395596370473593656809500224566606 u^6 bi9[21,7] = 28.2326688553901381296643680221162778846 u^7 bi9[21,8] = -16.035783011964464702761821641553971739 u^8 bi9[21,9] = 3.52057837559415010924534459738411867224 u^9 # ----------------------------------------- # # COEFFICIENTS OF bi9[22] bi9[22,1] = 0 u bi9[22,2] = 61.7513656850886811028625410193822166470 u^2 bi9[22,3] = -750.12689105528844734932299309395543253 u^3 bi9[22,4] = 3898.68558624484193098825600972584878432 u^4 bi9[22,5] = -10841.822728647457771543205430676141242 u^5 bi9[22,6] = 17321.5756646251661736613088744934995127 u^6 bi9[22,7] = -15937.704924918136911514034075812045455 u^7 bi9[22,8] = 7853.04762586250664965158337364697049736 u^8 bi9[22,9] = -1605.4056977967203049974482993035588813 u^9 # ----------------------------------------- # # COEFFICIENTS OF bi9[23] bi9[23,1] = 0 u bi9[23,2] = -1.9342980422240968080130439759208002892 u^2 bi9[23,3] = 147.588643166560364755815547685696411217 u^3 bi9[23,4] = -1181.4534621876819322477479227500066356 u^4 bi9[23,5] = 4124.81626831740245139843386194868654790 u^5 bi9[23,6] = -7708.5551065924472637822698176383416097 u^6 bi9[23,7] = 8052.34941989865008536988146615863307702 u^7 bi9[23,8] = -4443.0230550549520346267512917296872034 u^8 bi9[23,9] = 1010.21159049469242594065120030094021303 u^9 # ----------------------------------------- # # COEFFICIENTS OF bi9[24] bi9[24,1] = 0 u bi9[24,2] = .912642830786966650774488909524955581583 u^2 bi9[24,3] = -71.317106090194785769249593888661736245 u^3 bi9[24,4] = 695.456144607544639105995765905264616486 u^4 bi9[24,5] = -2620.8821934364843630835775612413187401 u^5 bi9[24,6] = 4957.00814115435023356389585992126611799 u^6 bi9[24,7] = -5038.2980855441275600980204854759417117 u^7 bi9[24,8] = 2633.33393505907862436876094759233267598 u^8 bi9[24,9] = -556.21347858095375473857942172246617795 u^9 # ----------------------------------------- # # COEFFICIENTS OF bi9[25] bi9[25,1] = 0 u bi9[25,2] = -.46544399707479916510507825569141210070 u^2 bi9[25,3] = 37.0909981044415418486565893300338539919 u^3 bi9[25,4] = -417.13822283398549912290149228106630232 u^4 bi9[25,5] = 1920.75879906118821707332162525764120727 u^5 bi9[25,6] = -4409.9601561495816960337202182818196327 u^6 bi9[25,7] = 5383.38350374976810011208128881165534148 u^7 bi9[25,8] = -3353.5185017727421417379849505934483413 u^8 bi9[25,9] = 839.849023837986277025652236012695285816 u^9 # ----------------------------------------- # # COEFFICIENTS OF bi9[26] bi9[26,1] = 0 u bi9[26,2] = -.11764517082552185612681029287052146921 u^2 bi9[26,3] = 9.48189481935429884137488032758066632784 u^3 bi9[26,4] = -114.91567582983193548201005631540452580 u^4 bi9[26,5] = 591.319329131246706800222487116556417765 u^5 bi9[26,6] = -1559.7654750865814938348967801088302956 u^6 bi9[26,7] = 2198.54584171864926697910295721673847716 u^7 bi9[26,8] = -1564.5620482737283422997168509115794570 u^8 bi9[26,9] = 440.013778691717020852050172967809238702 u^9 # #******************************************************** # Norms of low order INTERPOLANT coefficients on [0,2] # u Max norm 2-norm # ------------------------------------------------- .1000000000, .1227222739e-5, .3348508463e-5 .2000000000, .1615356496e-5, .3721626366e-5 .3000000000, .6667570096e-6, .1125720884e-5 .4000000000, .3385524146e-6, .1371991824e-5 .5000000000, -.4979316877e-6, .1866845810e-5 .6000000000, -.3855275900e-6, .1577140382e-5 .7000000000, -.5333350301e-6, .1597787028e-5 .8000000000, -.8981064587e-6, .1968485379e-5 .9000000000, -.6973183102e-6, .1411646512e-5 1.000000000, .1966258333e-56, .4681281745e-56 1.100000000, -.3763579804e-5, .8746517281e-5 1.200000000, -.2896131271e-4, .7960902000e-4 1.300000000, -.1153029259e-3, .3883837105e-3 1.400000000, .4370475809e-3, .1421009138e-2 1.500000000, .1370581071e-2, .4332655075e-2 1.600000000, .3719164661e-2, .1157394481e-1 1.700000000, .9047924268e-2, .2790864305e-1 1.800000000, .2019562961e-1, .6198056425e-1 1.900000000, .4203146751e-1, .1286244962 2.000000000, .8253193184e-1, .2521568348 # # ******************************************************** # Norms of high order INTERPOLANT error coefficients on [0,2] # u Max norm 2-norm # ------------------------------------------------------- .1000000000, .5235834963e-7, .1519257373e-6 .2000000000, .2095914343e-6, .5424950965e-6 .3000000000, .3292408958e-6, .7830811505e-6 .4000000000, .3445670411e-6, .7603468569e-6 .5000000000, .3009343118e-6, .5990805940e-6 .6000000000, .2568176444e-6, .4484576869e-6 .7000000000, .2143238993e-6, .3562383778e-6 .8000000000, .1418727695e-6, .2683013543e-6 .9000000000, .5417476349e-7, .2733305329e-6 1.000000000, .5849302357e-7, .3439902103e-6 1.100000000, .8714506430e-6, .3576174326e-5 1.200000000, .1131375469e-4, .4623422258e-4 1.300000000, .6872897170e-4, .2815480596e-3 1.400000000, .2923912562e-3, .1200914903e-2 1.500000000, .9964163854e-3, .4101897910e-2 1.600000000, .2905974139e-2, .1198624121e-1 1.700000000, .7539401031e-2, .3114936931e-1 1.800000000, .1784337214e-1, .7382582027e-1 1.900000000, .3919834562e-1, .1623811325 2.000000000, .8094769658e-1, .3356926326 # ********************************************************