matlab解微分方程实例1

韩鼎 020080031 ————————————————————————————————————————————— 1.对初值问题 ( )⎪⎩ ⎪⎨ ⎧ = ≤<−= 10 1020 u xu dx du 真解为 ，分别取步长 和 用四阶 Runge-Kutta 法计算。 xexu 20)( −= 2.0=h 1.0=h 四阶龙格-库塔法的绝对稳定区间是（-2.785，0）， 当取 时，2.0=h h 的值为 4−=hμ ，不再稳定区间。故结果是错误的。 当取 时，1.0=h h 的值为 2−=hμ ，在稳定区间，故结果是正确的，但偏差较大。 1u 2u 3u 4u 5u 6u 2.0=h 1 5 25 1.0=h 1.00000000 0.33333333 0.11111111 0.03703704 0.01234568 0.00411523 真值 1.00000000 0.13533528 0.01831564 0.00247875 0.00033546 0.00004540 7u 8u 9u 10u 11u 2.0=h 125 625 3125 1.0=h 0.00137174 0.00045725 0.00015242 0.00005081 0.00001694 真值 0.00000614 0.00000083 0.00000011 0.00000002 0.00000000 四阶 Runge-Kutta 1.0=h ───────────── x=[0:0.1:1]; u=ones(1,11); h=0.1; for i=1:10 k1=-20*u(i); k2=-20*(u(i)+h*k1/2); k3=-20*(u(i)+h*k2/2); k4=-20*(u(i)+h*k3); u(i+1)=u(i)+h*(k1+2*k2+2*k3+k4)/6; end 四阶 Runge-Kutta 2.0=h ───────────── x=[0:0.2:1]; u=ones(1,6); h=0.2; for i=1:5 k1=-20*u(i); k2=-20*(u(i)+h*k1/2); k3=-20*(u(i)+h*k2/2); k4=-20*(u(i)+h*k3); u(i+1)=u(i)+h*(k1+2*k2+2*k3+k4)/6; end 再取 和 用四阶的 Adams 外插法计算，并比较误差，说明为什么会产生这样的结果。 02.0=h 01.0=h 四阶的 Adams 外插法 02.0=h ───────────── format long； x=[0:0.02:1]; u=ones(1,51); f=ones(1,50) h=0.02; for m=1:4 k1=-20*u(m); k2=-20*(u(m)+h/2*k1); k3=-20*(u(m)+h/2*k2); k4=-20*(u(m)+h*k3); u(m+1)=u(m)+h*(k1+2*k2+2*k2+k4)/6; %龙格-库塔求初值 end for m=5:50 f(m)=-20*u(m); f(m-1)=-20*u(m-1); f(m-2)=-20*u(m-2); f(m-3)=-20*u(m-3); u(m+1)=u(m)+h*(55*f(m)-59*f(m-1) +37*f(m-2)-9*f(m-3))/24; %预测 end 0 0.2 0.4 0.6 0.8 1 −0.6 −0.4 −0.2 0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 四阶的 Adams 外插法 01.0=h ───────────── format long； x=[0:0.01:1]; u=ones(1,101); f=ones(1,100) h=0.01; for m=1:4 k1=-20*u(m); k2=-20*(u(m)+h/2*k1); k3=-20*(u(m)+h/2*k2); k4=-20*(u(m)+h*k3); u(m+1)=u(m)+h*(k1+2*k2+2*k2+k4)/6; %龙格-库塔求初值 end for m=5:100 f(m)=-20*u(m); f(m-1)=-20*u(m-1); f(m-2)=-20*u(m-2); f(m-3)=-20*u(m-3); u(m+1)=u(m)+h*(55*f(m)-59*f(m-1) +37*f(m-2)-9*f(m-3))/24; %预测 end 四阶 Adams 外插 方法 快递客服问题件处理详细方法山木方法pdf计算方法pdf华与华方法下载八字理论方法下载 的绝对稳定区间是（-0.3，0）， 当取 时，02.0=h h 的值为 4.0−=hμ ，不在稳定区间。故结果是错误的。 当取 时，1.0=h h 的值为 2.0−=hμ ，在稳定区间，故结果是正确的。 X h=0.2 h=0.1 真值 h=0.2 误差 h=0.1 误差 1 1.0000E+00 1.0000E+00 1.0000E+00 0.0000E+00 0.0000E+00 2 8.1940E-01 8.1873E-01 6.6925E-04 3 6.7573E-01 6.7142E-01 6.7032E-01 5.4133E-03 1.0963E-03 4 5.5016E-01 5.4881E-01 1.3469E-03 5 4.5662E-01 4.5080E-01 4.4933E-01 7.2866E-03 1.4710E-03 6 3.6911E-01 3.6788E-01 1.2334E-03 7 3.0855E-01 3.0230E-01 3.0119E-01 7.3561E-03 1.1092E-03 8 2.4749E-01 2.4660E-01 8.9648E-04 9 2.0850E-01 2.0269E-01 2.0190E-01 6.6012E-03 7.9515E-04 10 1.6595E-01 1.6530E-01 6.4995E-04 11 1.4056E-01 1.3591E-01 1.3534E-01 5.2278E-03 5.7268E-04 12 1.1127E-01 1.1080E-01 4.7058E-04 13 9.4956E-02 9.1129E-02 9.0718E-02 4.2380E-03 4.1106E-04 14 7.4612E-02 7.4274E-02 3.3876E-04 15 6.3842E-02 6.1104E-02 6.0810E-02 3.0322E-03 2.9375E-04 16 5.0030E-02 4.9787E-02 2.4265E-04 17 4.3288E-02 4.0971E-02 4.0762E-02 2.5255E-03 2.0914E-04 18 3.3546E-02 3.3373E-02 1.7307E-04 19 2.8914E-02 2.7472E-02 2.7324E-02 1.5901E-03 1.4841E-04 20 2.2494E-02 2.2371E-02 1.2298E-04 21 1.9850E-02 1.8421E-02 1.8316E-02 1.5341E-03 1.0500E-04 22 1.5083E-02 1.4996E-02 8.7102E-05 23 1.2968E-02 1.2351E-02 1.2277E-02 6.9099E-04 7.4097E-05 24 1.0113E-02 1.0052E-02 6.1509E-05 25 9.2625E-03 8.2819E-03 8.2297E-03 1.0328E-03 5.2160E-05 26 6.7813E-03 6.7379E-03 4.3322E-05 27 5.6205E-03 5.5532E-03 5.5166E-03 1.0392E-04 3.6637E-05 28 4.5470E-03 4.5166E-03 3.0440E-05 29 4.5569E-03 3.7235E-03 3.6979E-03 8.5899E-04 2.5682E-05 30 3.0489E-03 3.0276E-03 2.1342E-05 31 2.1399E-03 2.4967E-03 2.4788E-03 -3.3885E-04 1.7969E-05 32 2.0444E-03 2.0294E-03 1.4934E-05 33 2.5826E-03 1.6741E-03 1.6616E-03 9.2109E-04 1.2551E-05 34 1.3708E-03 1.3604E-03 1.0432E-05 35 3.5247E-04 1.1225E-03 1.1138E-03 -7.6131E-04 8.7528E-06 36 9.1916E-04 9.1188E-04 7.2745E-06 37 1.9329E-03 7.5268E-04 7.4659E-04 1.1863E-03 6.0951E-06 38 6.1632E-04 6.1125E-04 5.0652E-06 39 -7.6398E-04 5.0469E-04 5.0045E-04 -1.2644E-03 4.2386E-06 40 4.1326E-04 4.0973E-04 3.5219E-06 41 2.0071E-03 3.3841E-04 3.3546E-04 1.6716E-03 2.9438E-06 42 2.7710E-04 2.7465E-04 2.4456E-06 43 -1.7231E-03 2.2691E-04 2.2487E-04 -1.9479E-03 2.0421E-06 44 1.8580E-04 1.8411E-04 1.6962E-06 45 2.5911E-03 1.5215E-04 1.5073E-04 2.4403E-03 1.4150E-06 46 1.2458E-04 1.2341E-04 1.1751E-06 47 -2.8307E-03 1.0202E-04 1.0104E-04 -2.9318E-03 9.7951E-07 48 8.3537E-05 8.2724E-05 8.1324E-07 49 3.6756E-03 6.8406E-05 6.7729E-05 3.6079E-03 6.7736E-07 50 5.6014E-05 5.5452E-05 5.6225E-07 51 -4.3335E-03 4.5868E-05 4.5400E-05 -4.3789E-03 4.6799E-07 52 3.7559E-05 3.7170E-05 3.8836E-07 53 5.3875E-03 3.0756E-05 3.0432E-05 5.3571E-03 3.2305E-07 54 2.5184E-05 2.4916E-05 2.6802E-07 55 -6.5036E-03 2.0622E-05 2.0400E-05 -6.5240E-03 2.2282E-07 56 1.6887E-05 1.6702E-05 1.8482E-07 57 7.9795E-03 1.3828E-05 1.3674E-05 7.9658E-03 1.5356E-07 58 1.1323E-05 1.1195E-05 1.2734E-07 59 -9.7026E-03 9.2718E-06 9.1661E-06 -9.7118E-03 1.0576E-07 60 7.5922E-06 7.5046E-06 8.7678E-08 61 1.1857E-02 6.2170E-06 6.1442E-06 1.1850E-02 7.2782E-08 62 5.0908E-06 5.0305E-06 6.0326E-08 63 -1.4449E-02 4.1686E-06 4.1186E-06 -1.4453E-02 5.0056E-08 64 3.4135E-06 3.3720E-06 4.1479E-08 65 1.7635E-02 2.7952E-06 2.7608E-06 1.7632E-02 3.4404E-08 66 2.2888E-06 2.2603E-06 2.8502E-08 67 -2.1506E-02 1.8742E-06 1.8506E-06 -2.1507E-02 2.3632E-08 68 1.5347E-06 1.5151E-06 1.9574E-08 69 2.6238E-02 1.2567E-06 1.2405E-06 2.6237E-02 1.6224E-08 70 1.0291E-06 1.0156E-06 1.3435E-08 71 -3.2003E-02 8.4266E-07 8.3153E-07 -3.2004E-02 1.1132E-08 72 6.9001E-07 6.8080E-07 9.2159E-09 73 3.9041E-02 5.6502E-07 5.5739E-07 3.9040E-02 7.6338E-09 74 4.6267E-07 4.5635E-07 6.3187E-09 75 -4.7622E-02 3.7886E-07 3.7363E-07 -4.7622E-02 5.2324E-09 76 3.1023E-07 3.0590E-07 4.3302E-09 77 5.8092E-02 2.5404E-07 2.5045E-07 5.8092E-02 3.5848E-09 78 2.0802E-07 2.0505E-07 2.9660E-09 79 -7.0863E-02 1.7034E-07 1.6788E-07 -7.0863E-02 2.4548E-09 80 1.3948E-07 1.3745E-07 2.0307E-09 81 8.6442E-02 1.1422E-07 1.1254E-07 8.6442E-02 1.6803E-09 82 9.3526E-08 9.2136E-08 1.3897E-09 83 -1.0545E-01 7.6584E-08 7.5435E-08 -1.0545E-01 1.1497E-09 84 6.2711E-08 6.1761E-08 9.5068E-10 85 1.2863E-01 5.1352E-08 5.0565E-08 1.2863E-01 7.8628E-10 86 4.2049E-08 4.1399E-08 6.5008E-10 87 -1.5690E-01 3.4432E-08 3.3895E-08 -1.5690E-01 5.3755E-10 88 2.8195E-08 2.7751E-08 4.4436E-10 89 1.9140E-01 2.3088E-08 2.2720E-08 1.9140E-01 3.6736E-10 90 1.8906E-08 1.8602E-08 3.0362E-10 91 -2.3348E-01 1.5481E-08 1.5230E-08 -2.3348E-01 2.5097E-10 92 1.2677E-08 1.2469E-08 2.0739E-10 93 2.8480E-01 1.0380E-08 1.0209E-08 2.8480E-01 1.7139E-10 94 8.5000E-09 8.3584E-09 1.4161E-10 95 -3.4742E-01 6.9603E-09 6.8433E-09 -3.4742E-01 1.1701E-10 96 5.6995E-09 5.6028E-09 9.6659E-11 97 4.2379E-01 4.6670E-09 4.5872E-09 4.2379E-01 7.9854E-11 98 3.8216E-09 3.7557E-09 6.5957E-11 99 -5.1696E-01 3.1294E-09 3.0749E-09 -5.1696E-01 5.4481E-11 100 2.5625E-09 2.5175E-09 4.4994E-11 101 6.3061E-01 2.0983E-09 2.0612E-09 6.3061E-01 3.7159E-11 2.对初值 ( )⎪⎩ ⎪⎨ ⎧ = ≤<= 10 204 u xux dx du ，取步长 05.0=h 进行计算（真解为） 分别用对应于 四阶 adams 方法的预测-校正法和修正的预测-校正法计算并比较误差。 22 )1()( xxu += 四阶 Adams 方法的预测-校正法 ──────────────────────────── x=[0:0.05:2]; format long； u=ones(1,41); f=ones(1,41) h=0.05; for m=1:4 k1=4*x(m)*sqrt(u(m)); k2=4*(x(m)+h/2)*sqrt(u(m)+k1*h/2); k3=4*(x(m)+h/2)*sqrt(u(m)+k2*h/2); k4=4*(x(m)+h)*sqrt(u(m)+k3*h); u(m+1)=u(m)+h*(k1+2*k2+2*k2+k4)/6; %龙格-库塔求初值 end for m=5:40 f(m)=4*x(m)*sqrt(u(m)); f(m-1)=4*x(m-1)*sqrt(u(m-1)); f(m-2)=4*x(m-2)*sqrt(u(m-2)); f(m-3)=4*x(m-3)*sqrt(u(m-3)); u(m+1)=u(m)+h*(55*f(m)-59*f(m-1)+37*f(m-2)-9*f(m-3))/24; % 预 测 f(m+1)=4*x(m+1)*sqrt(u(m+1)); u(m+1)=u(m)+h*(9*f(m+1)+19*f(m)-5*f(m-1)+f(m-2))/24;%校正 f(m+1)=4*x(m+1)*sqrt(u(m+1)); end 四阶 Adams 方法的修正预测-校正法 ──────────────────────────── x=[0:0.05:2]; format long； u=ones(1,41); up=ones(1,41); %预测值 uc=ones(1,41); %修正值 f=ones(1,41); n=ones(1,41); h=0.05; for m=1:4 k1=4*x(m)*sqrt(u(m)); k2=4*(x(m)+h/2)*sqrt(u(m)+k1*h/2); k3=4*(x(m)+h/2)*sqrt(u(m)+k2*h/2); k4=4*(x(m)+h)*sqrt(u(m)+k3*h); u(m+1)=u(m)+h*(k1+2*k2+2*k2+k4)/6;%龙格-库塔求初值 end for m=5:40 f(m)=4*x(m)*sqrt(u(m)); f(m-1)=4*x(m-1)*sqrt(u(m-1)); f(m-2)=4*x(m-2)*sqrt(u(m-2)); f(m-3)=4*x(m-3)*sqrt(u(m-3)); up(m+1)=u(m)+h*(55*f(m)-59*f(m-1)+37*f(m-2)-9*f(m-3))/24;%预测 n(m+1)=up(m+1)+(uc(m-1)-up(m-1)); f(m+1)=4*x(m+1)*sqrt(n(m+1)); uc(m+1)=u(m)+h*(9*f(m+1)+19*f(m)-5*f(m-1)+f(m-2))/24;%校正 u(m+1)=uc(m+1)-19*(uc(m+1)-up(m+1)); f(m+1)=4*x(m+1)*sqrt(u(m+1)); end 预测校正 修正预测校正 真值 预测校正误差 修正预测校正误差 1 1.0000E+00 1.0000E+00 1.0000E+00 0.0000E+00 0.0000E+00 2 1.0050E+00 1.0050E+00 1.0050E+00 -2.0833E-06 -2.0833E-06 3 1.0201E+00 1.0201E+00 1.0201E+00 -8.3800E-06 -8.3800E-06 4 1.0455E+00 1.0455E+00 1.0455E+00 -1.9105E-05 -1.9105E-05 5 1.0816E+00 1.0816E+00 1.0816E+00 -3.4649E-05 -3.4649E-05 6 1.1289E+00 1.1289E+00 1.1289E+00 -3.5436E-05 -3.8766E-05 7 1.1881E+00 1.1881E+00 1.1881E+00 -3.6342E-05 -3.6907E-05 8 1.2600E+00 1.2600E+00 1.2600E+00 -3.7426E-05 -3.6556E-05 9 1.3456E+00 1.3456E+00 1.3456E+00 -3.8677E-05 -3.8753E-05 10 1.4460E+00 1.4460E+00 1.4460E+00 -4.0094E-05 -4.1241E-05 11 1.5625E+00 1.5625E+00 1.5625E+00 -4.1677E-05 -4.2777E-05 12 1.6965E+00 1.6965E+00 1.6965E+00 -4.3428E-05 -4.3938E-05 13 1.8496E+00 1.8496E+00 1.8496E+00 -4.5345E-05 -4.5650E-05 14 2.0235E+00 2.0235E+00 2.0235E+00 -4.7429E-05 -4.8037E-05 15 2.2201E+00 2.2200E+00 2.2201E+00 -4.9680E-05 -5.0614E-05 16 2.4414E+00 2.4414E+00 2.4414E+00 -5.2097E-05 -5.3042E-05 17 2.6895E+00 2.6895E+00 2.6896E+00 -5.4681E-05 -5.5448E-05 18 2.9669E+00 2.9669E+00 2.9670E+00 -5.7432E-05 -5.8136E-05 19 3.2760E+00 3.2760E+00 3.2761E+00 -6.0349E-05 -6.1190E-05 20 3.6194E+00 3.6194E+00 3.6195E+00 -6.3433E-05 -6.4450E-05 21 3.9999E+00 3.9999E+00 4.0000E+00 -6.6684E-05 -6.7757E-05 22 4.4204E+00 4.4204E+00 4.4205E+00 -7.0102E-05 -7.1135E-05 23 4.8840E+00 4.8840E+00 4.8841E+00 -7.3686E-05 -7.4713E-05 24 5.3939E+00 5.3939E+00 5.3940E+00 -7.7437E-05 -7.8552E-05 25 5.9535E+00 5.9535E+00 5.9536E+00 -8.1355E-05 -8.2590E-05 26 6.5663E+00 6.5663E+00 6.5664E+00 -8.5439E-05 -8.6746E-05 27 7.2360E+00 7.2360E+00 7.2361E+00 -8.9690E-05 -9.1022E-05 28 7.9664E+00 7.9664E+00 7.9665E+00 -9.4108E-05 -9.5476E-05 29 8.7615E+00 8.7615E+00 8.7616E+00 -9.8693E-05 -1.0014E-04 30 9.6254E+00 9.6254E+00 9.6255E+00 -1.0344E-04 -1.0499E-04 31 1.0562E+01 1.0562E+01 1.0563E+01 -1.0836E-04 -1.0999E-04 32 1.1577E+01 1.1577E+01 1.1577E+01 -1.1345E-04 -1.1513E-04 33 1.2673E+01 1.2673E+01 1.2674E+01 -1.1870E-04 -1.2045E-04 34 1.3857E+01 1.3857E+01 1.3857E+01 -1.2412E-04 -1.2595E-04 35 1.5132E+01 1.5132E+01 1.5132E+01 -1.2970E-04 -1.3163E-04 36 1.6504E+01 1.6504E+01 1.6504E+01 -1.3545E-04 -1.3747E-04 37 1.7977E+01 1.7977E+01 1.7978E+01 -1.4137E-04 -1.4347E-04 38 1.9558E+01 1.9558E+01 1.9559E+01 -1.4746E-04 -1.4964E-04 39 2.1252E+01 2.1252E+01 2.1252E+01 -1.5371E-04 -1.5598E-04 40 2.3064E+01 2.3064E+01 2.3064E+01 -1.6013E-04 -1.6250E-04 41 2.5000E+01 2.5000E+01 2.5000E+01 -1.6671E-04 -1.6918E-04 3.对线性二步法 })1()3{( 2 1)1( 111 −−+ +−−=++− mmmmm fafahauuau ①当 和 时各为几阶方法 0=a 5−=a ②当 和 带入 公式 小学单位换算公式大全免费下载公式下载行测公式大全下载excel公式下载逻辑回归公式下载 分别计算，结果如何， 分析 定性数据统计分析pdf销售业绩分析模板建筑结构震害分析销售进度分析表京东商城竞争战略分析 产生的原因。 0=a 5−=a 线性二步法a 0= ────────────── x=[0:0.05:2]; u=ones(1,41); f=ones(1,41); h=0.05; k1=4*x(1)*sqrt(u(1)); k2=4*(x(1)+h/2)*sqrt(u(1)+k1*h/2); k3=4*(x(1)+h/2)*sqrt(u(1)+k2*h/2); k4=4*(x(1)+h)*sqrt(u(1)+k3*h); u(2)=u(1)+h*(k1+2*k2+2*k2+k4)/6; for i=2:40 f(i)=4*x(i)*sqrt(u(i)); f(i-1)=4*x(i-1)*sqrt(u(i-1)); u(i+1)=u(i)+h*(3*f(i)-f(i-1))/2; end 线性二步法 5−=a ────────────── x=[0:0.05:2]; u=ones(1,41); f=ones(1,41); h=0.05; k1=4*x(1)*sqrt(u(1)); k2=4*(x(1)+h/2)*sqrt(u(1)+k1*h/2); k3=4*(x(1)+h/2)*sqrt(u(1)+k2*h/2); k4=4*(x(1)+h)*sqrt(u(1)+k3*h); u(2)=u(1)+h*(k1+2*k2+2*k2+k4)/6; for i=2:40 f(i)=4*x(i)*sqrt(u(i)); f(i-1)=4*x(i-1)*sqrt(u(i-1)); u(i+1)=-4*u(i)+5*u(i-1)+h*(4*f(i)+ 2*f(i-1)); a=0 a=-5 真值 a=0 误差 a=-5 误差 1 1.0000E+00 1.00000000000000 + 0.00000000000000i 1.0000E+ 00 0.0000E+0 0 * 2 1.0050E+00 1.00500416666342 + -0.00000000000000i 1.0050E+ 00 -2.0833E- 06 * 3 1.0200E+00 1.02008329178346 + -0.00000000000000i 1.0201E+ 00 -5.8349E- 05 * 4 1.0453E+00 1.04553698368760 + -0.00000000000000i 1.0455E+ 00 -1.7796E- 04 * 5 1.0812E+00 1.08136999674027 + -0.00000000000000i 1.0816E+ 00 -3.6284E- 04 * 6 1.1283E+00 1.12993813969308 + -0.00000000000000i 1.1289E+ 00 -6.1575E- 04 * 7 1.1872E+00 1.18288567502636 + -0.00000000000000i 1.1881E+ 00 -9.4025E- 04 * 8 1.2587E+00 1.28547186142071 + -0.00000000000000i 1.2600E+ 00 -1.3406E- 03 * 9 1.3438E+00 1.22051381094603 + -0.00000000000000i 1.3456E+ 00 -1.8216E- 03 * 1 0 1.4436E+ 00 2.05756014933536 + -0.00000000000000i 1.4460E+ 00 -2.3887E- 03 * 1 1 1.5595E+ 00 -1.43451742931644 + -0.00000000000000i 1.5625E+ 00 -3.0477E- 03 * 1 2 1.6927E+ 00 16.2840660271083 - 0.479085366809122i 1.6965E+ 00 -3.8047E- 03 * 1 3 1.8449E+ 00 -70.5331043051500 + 1.65068275220330i 1.8496E+ 00 -4.6660E- 03 * 1 4 2.0179E+ 00 364.487788957862 - 4.97970844033798i 2.0235E+ 00 -5.6384E- 03 * 1 5 2.2134E+ 00 -1800.66524876411 + 30.1201862512429i 2.2201E+ 00 -6.7284E- 03 * 1 6 2.4335E+ 00 9030.26260152458 - 121.649185789069i 2.4414E+ 00 -7.9427E- 03 * 1 7 2.6803E+ 00 -45067.2593705767 + 648.695644208090i 2.6896E+ 00 -9.2882E- 03 * 1 8 2.9562E+ 00 225449.837242735 - 3067.35108015850i 2.9670E+ 00 -1.0771E- 02 * 1 9 3.2637E+ 00 -1126812.27485840 + 15578.6208835926i 3.2761E+ 00 -1.2399E- 02 * 2 0 3.6053E+ 00 5634665.00985976 - 76888.0286824452i 3.6195E+ 00 -1.4177E- 02 * 2 1 3.9839E+ 00 -28170914.6840050 + 385815.065066401i 4.0000E+ 00 -1.6113E- 02 * 2 2 4.4023E+ 00 140857914.904976 - 1923460.35590415i 4.4205E+ 00 -1.8212E- 02 * 2 3 4.8636E+ 00 -704276248.849292 + 9624971.78366803i 4.8841E+ 00 -2.0481E- 02 * 2 4 5.3711E+ 00 3521399714.32554 - 48093868.7829759i 5.3940E+ 00 -2.2925E- 02 * 2 5 5.9280E+ 00 -17606925426.4150 + 240511638.329105i 5.9536E+ 00 -2.5551E- 02 * 2 6 6.5380E+ 00 88034728444.9727 - 1202388697.21533i 6.5664E+ 00 -2.8363E- 02 * 2 7 7.2047E+ 00 -440173243763.573 + 6012174647.53694i 7.2361E+ 00 -3.1368E- 02 * 2 8 7.9319E+ 00 2200866770347.93 - 30059953079.4702i 7.9665E+ 00 -3.4570E- 02 * 2 9 8.7236E+ 00 -11004331695601.8 + 150301019619.033i 8.7616E+ 00 -3.7976E- 02 * 3 0 9.5839E+ 00 55021661460644.9 - 751500133906.493i 9.6255E+ 00 -4.1589E- 02 * 3 1 1.0517E+ 01 -275108295703218 + 3757507432680.23i 1.0563E+ 01 -4.5415E- 02 * 3 2 1.1528E+ 01 1.37554149455436e+15 - 18787510525502.5i 1.1577E+ 01 -4.9458E- 02 * 3 3 1.2620E+ 01 -6.87770741067501e+15 + 93937588903415.6i 1.2674E+ 01 -5.3724E- 02 * 3 4 1.3799E+ 01 3.43885371391920e+16 - 469687802242842i 1.3857E+ 01 -5.8215E- 02 * 3 5 1.5069E+ 01 -1.71942685364992e+17 + 2.34843920489451e+15i 1.5132E+ 01 -6.2938E- 02 * 3 6 1.6436E+ 01 8.59713427282174e+17 - 1.17421952676777e+16i 1.6504E+ 01 -6.7895E- 02 * 3 7 1.7905E+ 01 -4.29856713465361e+18 + 5.87109773682939e+16i 1.7978E+ 01 -7.3091E- 02 * 3 8 1.9480E+ 01 2.14928356756948e+19 - 2.93554882830377e+17i 1.9559E+ 01 -7.8530E- 02 * 3 9 2.1168E+ 01 -1.07464178369175e+20 + 1.46777441960893e+18i 2.1252E+ 01 -8.4216E- 02 * 4 0 2.2974E+ 01 5.37320891858714e+20 - 7.33887207685360e+18i 2.3064E+ 01 -9.0152E- 02 * 4 1 2.4904E+ 01 -2.68660445924452e+21 + 3.66943604130908e+19i 2.5000E+ 01 -9.6342E- 02 *