@@ -84,9 +84,8 @@ function Uk_poly10(p, v, p2, even_odd_poly)
84
84
u8 = evalpoly (p2, (6.074042001273483 , - 493.91530477308805 , 7109.514302489364 , - 41192.65496889755 , 122200.46498301747 , - 203400.17728041555 , 192547.00123253153 , - 96980.59838863752 , 20204.29133096615 ))
85
85
u9 = evalpoly (p2, (24.380529699556064 , - 2499.8304818112097 , 45218.76898136273 , - 331645.17248456355 , 1.2683652733216248e6 , - 2.8135632265865337e6 , 3.763271297656404e6 , - 2.998015918538107e6 , 1.3117636146629772e6 , - 242919.18790055133 ))
86
86
u10 = evalpoly (p2, (110.01714026924674 , - 13886.08975371704 , 308186.40461266245 , - 2.7856181280864547e6 , 1.3288767166421818e7 , - 3.7567176660763346e7 , 6.634451227472903e7 , - 7.410514821153265e7 , 5.095260249266464e7 , - 1.9706819118432228e7 , 3.284469853072038e6 ))
87
- u11 = evalpoly (p2, (551.3358961220206 , - 84005.4336030241 , 2.2437681779224495e6 , - 2.4474062725738727e7 , 1.420629077975331e8 , - 4.9588978427503026e8 , 1.1068428168230145e9 , - 1.6210805521083372e9 , 1.5535968995705802e9 , - 9.394623596815784e8 , 3.2557307418576574e8 , - 4.932925366450996e7 ))
88
87
89
- Poly = (u0, u1, u2, u3, u4, u5, u6, u7, u8, u9, u10, u11 )
88
+ Poly = (u0, u1, u2, u3, u4, u5, u6, u7, u8, u9, u10)
90
89
return even_odd_poly (- p/ v, Poly)
91
90
end
92
91
@@ -112,9 +111,8 @@ function Uk_poly20(p, v, p2, even_odd_poly)
112
111
u18 = evalpoly (p2, (4.259392165047669e8 , - 1.722832387173505e11 , 1.2030115826419191e13 , - 3.4396530474307594e14 , 5.335106978708839e15 , - 5.1605093193485224e16 , 3.37667624979061e17 , - 1.5736434765189599e18 , 5.402894876715982e18 , - 1.3970803516443374e19 , 2.757282981650519e19 , - 4.178861444656839e19 , 4.859942729324836e19 , - 4.301555703831444e19 , 2.846521225167657e19 , - 1.3639420410571592e19 , 4.47020096401231e18 , - 8.966114215270463e17 , 8.30195760673191e16 ))
113
112
u19 = evalpoly (p2, (3.8362551802304335e9 , - 1.7277040123529995e12 , 1.3412416915180639e14 , - 4.2619355104268985e15 , 7.351663610930971e16 , - 7.921651119323832e17 , 5.789887667664653e18 , - 3.025566598990372e19 , 1.1707490535797259e20 , - 3.434621399768417e20 , 7.756704953461136e20 , - 1.360203777284994e21 , 1.8571089321463453e21 , - 1.9677247077053125e21 , 1.6016898573693598e21 , - 9.824438427689858e20 , 4.392792200888712e20 , - 1.351217503435996e20 , 2.5563802960529236e19 , - 2.242438856186775e18 ))
114
113
u20 = evalpoly (p2, (3.646840080706556e10 , - 1.818726203851104e13 , 1.5613123930484672e15 , - 5.48403360388329e16 , 1.0461721131134344e18 , - 1.2483700995047234e19 , 1.0126774169536592e20 , - 5.8917941350694964e20 , 2.548961114664972e21 , - 8.405915817108351e21 , 2.1487414815055883e22 , - 4.302534303482379e22 , 6.783661642951883e22 , - 8.423222750084323e22 , 8.19433100543513e22 , - 6.173206302884415e22 , 3.528435843903409e22 , - 1.4787743528433614e22 , 4.285296082829494e21 , - 7.671943936729004e20 , 6.393286613940837e19 ))
115
- u21 = evalpoly (p2, (3.6490108188498334e11 , - 2.0052440123627112e14 , 1.894406984252143e16 , - 7.319501491566134e17 , 1.5365025218443373e19 , - 2.0197335419300872e20 , 1.8081594057131945e21 , - 1.1640246461465369e22 , 5.591591380366263e22 , - 2.0566149136271542e23 , 5.8965434619782445e23 , - 1.3337178907798302e24 , 2.3967237744351682e24 , - 3.430872898515746e24 , 3.905264103536985e24 , - 3.511096528332644e24 , 2.461506085403875e24 , - 1.3170969618092387e24 , 5.194289094766812e23 , - 1.4228394823321413e23 , 2.417461500896379e22 , - 1.91862023880665e21 ))
116
114
117
- Poly = (u0, u1, u2, u3, u4, u5, u6, u7, u8, u9, u10, u11, u12, u13, u14, u15, u16, u17, u18, u19, u20, u21 )
115
+ Poly = (u0, u1, u2, u3, u4, u5, u6, u7, u8, u9, u10, u11, u12, u13, u14, u15, u16, u17, u18, u19, u20)
118
116
return even_odd_poly (- p/ v, Poly)
119
117
end
120
118
0 commit comments