北航数值分析大作业第三题(fortran).pdf

《北航数值分析大作业第三题(fortran).pdf》由会员分享，可在线阅读，更多相关《北航数值分析大作业第三题(fortran).pdf（24页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、“数值分析数值分析“计算实习大作业第三题计算实习大作业第三题SY1415215孔维鹏一、计算说明一、计算说明1、将=0.08，=0.5+0.05j分别代入方程组（A.3）得到关于 t，u，v，w 的的方程组，调用离散牛顿迭代子函数求出与，对应的，。2、调用分片二次代数插值子函数在点(,)处插值得到(,)=(,)，得到数表(,(,)。3、对于=1,2,3,4，分别调用最小二乘拟合子函数计算系数矩阵及误差，直到满足精度，即求得最小的 k 值及系数矩阵。4、将=0.1，=0.5+0.2j分别代入方程组（A.3）得到关于，的的方程组，调用离散牛顿迭代子函数求出与，对应的，调用分片二次代数插值子函数在点

2、(,)处插值得到(,)=(,)；调用步骤 3 中求得的系数矩阵求得(,)，打印数表(,(,),(,)。二、源程序（二、源程序（FORTRANFORTRAN）PROGRAMPROGRAM SY1415215DIMENSIONDIMENSION X(11),Y(21),T(6),U(6),Z(6,6),UX(11,21),TY(11,21),FXY(11,21),C(6,6)DIMENSIONDIMENSION X1(8),Y1(5),FXY1(8,5),PXY1(8,5),UX1(8,5),TY1(8,5)REAL(8)REAL(8)X,Y,T,U,Z,FXY,UX,TY,C,E,X1,Y1,F

3、XY1,PXY1,UX1,TY1OPENOPEN(1,FILE=第三题计算结果.TXT)DODO I=1,11 X(I)=0.08*(I-1)ENDDOENDDODODO I=1,21 Y(I)=0.5+0.05*(I-1)ENDDOENDDO!*求解非线性方程组，得到z=f（t，u）的函数*DODO I=1,11DODO J=1,21CALLCALL DISNEWTON_NONLINEAR(X(I),Y(J),UX(I,J),TY(I,J)ENDDOENDDOENDDOENDDO!*分片二次插值得到z=f（x，y）*DODO I=1,11DODO J=1,21CALLCALL INTERPO

4、LATION(UX(I,J),TY(I,J),FXY(I,J)ENDDOENDDOENDDOENDDOWRITEWRITE(1,(数表(x,y,f(x,y):)WRITEWRITE(1,(3X,X,7X,Y,10X,F(X,Y)DODO I=1,11DODO J=1,21WRITEWRITE(1,(1X,F5.2,2X,F5.3,2X,E20.13)X(I),Y(J),FXY(I,J)ENDDOENDDOWRITEWRITE(1,()ENDDOENDDO!*最小二乘拟合得到P（x，y）*N=11M=21WRITEWRITE(1,(,K和 分别为：)DODO K=1,20CALLCALL LSF

5、ITTING(X,Y,FXY,C,N,M,K,K,E)WRITEWRITE(1,(I3,2X,E20.13)K-1,EIFIF(E10E-7)EXITEXITENDDOENDDOWRITEWRITE(1,()WRITEWRITE(1,(系数矩阵Crs（按行）为：)DODO I=1,KDODO J=1,KWRITEWRITE(1,(E20.13,2X,)C(I,J)ENDDOENDDOWRITEWRITE(1,()WRITEWRITE(*,()ENDDOENDDODODO I=1,8 X1(I)=0.1*IENDDOENDDODODO J=1,5 Y1(J)=0.5+0.2*JENDDOENDD

6、ODODO I=1,8DODO J=1,5CALLCALL DISNEWTON_NONLINEAR(X1(I),Y1(J),UX1(I,J),TY1(I,J)ENDDOENDDOENDDOENDDODODO I=1,8DODO J=1,5CALLCALL INTERPOLATION(UX1(I,J),TY1(I,J),FXY1(I,J)ENDDOENDDOENDDOENDDOPXY1=0DODO I=1,8DODO J=1,5DODO II=1,KDODO JJ=1,K PXY1(I,J)=PXY1(I,J)+C(II,JJ)*(X1(I)*(II-1)*(Y1(J)*(JJ-1)ENDDOE

7、NDDOENDDOENDDOENDDOENDDOENDDOENDDOWRITEWRITE(1,()WRITEWRITE(1,(数表（x,y,f(x,y),p(x,y):)WRITEWRITE(1,(2X,X,6X,Y,12X,F(X,Y),14X,P(X,Y)DODO I=1,8DODO J=1,5WRITEWRITE(1,(F5.3,2X,F5.3,2X,E20.13,2X,E20.13)X1(I),Y1(J),FXY1(I,J),PXY1(I,J)ENDDOENDDOWRITEWRITE(1,()ENDDOENDDOCLOSECLOSE(1)ENDEND!*用离散牛顿法求解非线性方程组*S

8、UBROUTINESUBROUTINE DISNEWTON_NONLINEAR(X1,Y1,U,T)PARAMETERPARAMETER(N=4)REALREAL EPS!EPS为迭代精度，M为最大迭代次数DIMENSIONDIMENSION X(N),H(N),Y(N),JA(N,N),E(N),XK(N)REAL(8)REAL(8)JA,X,H,Y,E,XK,U,T,V,W,X1,Y1,E1,E2F1(T,U,V,W)=0.5*COS(T)+U+V+W-X1-2.67F2(T,U,V,W)=T+0.5*SIN(U)+V+W-Y1-1.07F3(T,U,V,W)=0.5*T+U+COS(V)

9、+W-X1-3.74F4(T,U,V,W)=T+0.5*U+V+SIN(W)-Y1-0.79EPS=10E-12M=100X=1.0DODO K=1,M H=1!计算Y=F（x）Y(1)=F1(X(1),X(2),X(3),X(4)Y(2)=F2(X(1),X(2),X(3),X(4)Y(3)=F3(X(1),X(2),X(3),X(4)Y(4)=F4(X(1),X(2),X(3),X(4)!计算JA（N,N）E=0.0DODO I=1,NDODO J=1,NDODO JJ=1,NIFIF(JJ=J)THENTHEN E(JJ)=X(JJ)+H(JJ)ELSEELSE E(JJ)=X(JJ)E

10、NDIFENDIFENDDOENDDOIFIF(I=1)THENTHEN JA(I,J)=(F1(E(1),E(2),E(3),E(4)-F1(X(1),X(2),X(3),X(4)/H(J)ELSEIFELSEIF(I=2)THENTHEN JA(I,J)=(F2(E(1),E(2),E(3),E(4)-F2(X(1),X(2),X(3),X(4)/H(J)ELSEIFELSEIF(I=3)THENTHEN JA(I,J)=(F3(E(1),E(2),E(3),E(4)-F3(X(1),X(2),X(3),X(4)/H(J)ELSEIFELSEIF(I=4)THENTHEN JA(I,J)=

11、(F4(E(1),E(2),E(3),E(4)-F4(X(1),X(2),X(3),X(4)/H(J)ENDIFENDIFENDDOENDDOENDDOENDDO!求解线性方程组CALLCALL GAUSS(JA,XK,-Y,N)!判断精度CALLCALL NORM(XK,N,E1)CALLCALL NORM(X,N,E2)IFIF(E1/E2TA).OR.(A(L,K)=TA)THENTHENTA=A(L,K)TL=LDODO J=K,NT(K,J)=A(K,J)A(K,J)=A(TL,J)A(TL,J)=T(K,J)ENDDOENDDOTB(K)=B(K)B(K)=B(TL)B(TL)=T

12、B(K)ENDIFENDIFENDDOENDDODODO I=K+1,NM(I,K)=A(I,K)/A(K,K)A(I,K)=0DODO J=K+1,NA(I,J)=A(I,J)-M(I,K)*A(K,J)ENDDOENDDOB(I)=B(I)-M(I,K)*B(K)ENDDOENDDOENDDOENDDO!回代过程X(N)=B(N)/A(N,N)DODO K=N-1,1,-1S=0.0DODO J=K+1,NS=S+A(K,J)*X(J)ENDDOENDDOX(K)=(B(K)-S)/A(K,K)ENDDOENDDORETURNRETURNENDEND!*求向量的无穷范数*SUBROUTIN

13、ESUBROUTINE NORM(X,N,A)DIMENSIONDIMENSION X(N)REAL(8)REAL(8)X,AA=ABS(X(1)DODO I=2,NIFIF(ABS(X(I)ABS(X(I-1)THENTHEN A=ABS(X(I)ENDIFENDIFENDDOENDDORETURNRETURNENDEND!*分片二次代数插值*SUBROUTINESUBROUTINE INTERPOLATION(U,V,W)PARAMETERPARAMETER(N=6,M=6)DIMENSIONDIMENSION X(N),Y(M),Z(M,N),LK(3),LR(3)REAL(8)REAL

14、(8)X,Y,Z,H,TREAL(8)REAL(8)U,V,W,LK,LR!U,V分别为插值点处的坐标，W为插值结果INTEGERINTEGER R!*数据赋值*DATADATA Y/0.0,0.2,0.4,0.6,0.8,1.0/DATADATA X/0.0,0.4,0.8,1.2,1.6,2.0/DATADATA Z/-0.5,-0.42,-0.18,0.22,0.78,1.5,&-0.34,-0.5,-0.5,-0.34,-0.02,0.46,&0.14,-0.26,-0.5,-0.58,-0.5,-0.26,&0.94,0.3,-0.18,-0.5,-0.66,-0.66,&2.06,

15、1.18,0.46,-0.1,-0.5,-0.74,&3.5,2.38,1.42,0.62,-0.02,-0.5/H=0.4T=0.2!*计算K,R*IFIF(UX(N-1)-H/2)THENTHEN K=N-1ELSEELSEDODO I=3,N-2IFIF(UX(I)-H/2).AND.(U=X(I)+H/2)THENTHEN K=IENDIFENDIFENDDOENDDOENDIFENDIFIFIF(VY(M-1)-T/2)THENTHEN R=M-1ELSEELSEDODO J=3,M-2IFIF(VY(J)-T/2).AND.(VN)P=NIFIF(P20)P=20IFIF(QM)Q

16、=MIFIF(Q20)Q=20XX=0YY=0D1=NAPX(1)=0.0DODO I=1,N APX(1)=APX(1)+X(I)ENDDOENDDOAPX(1)=APX(1)/D1DODO J=1,M V(1,J)=0.0DODO I=1,N V(1,J)=V(1,J)+Z(I,J)ENDDOENDDO V(1,J)=V(1,J)/D1ENDDOENDDOIFIF(P1)THENTHEN D2=0.0 APX(2)=0.0DODO I=1,N G=X(I)-APX(1)D2=D2+G*G APX(2)=APX(2)+(X(I)-XX)*G*GENDDOENDDO APX(2)=APX(2)

17、/D2 BX(2)=D2/D1DODO J=1,M V(2,J)=0.0DODO I=1,N G=X(I)-APX(1)V(2,J)=V(2,J)+Z(I,J)*GENDDOENDDO V(2,J)=V(2,J)/D2ENDDOENDDO D1=D2ENDIFENDIFDODO K=3,P D2=0.0 APX(K)=0.0DODO J=1,M V(K,J)=0.0ENDDOENDDODODO I=1,N G1=1.0 G2=X(I)-APX(1)DODO J=3,K G=(X(I)-APX(J-1)*G2-BX(J-1)*G1 G1=G2 G2=GENDDOENDDO D2=D2+G*G A

18、PX(K)=APX(K)+X(I)*G*GDODO J=1,M V(K,J)=V(K,J)+Z(I,J)*GENDDOENDDOENDDOENDDODODO J=1,M V(K,J)=V(K,J)/D2ENDDOENDDO APX(K)=APX(K)/D2 BX(K)=D2/D1 D1=D2ENDDOENDDOD1=MAPY(1)=0.0DODO I=1,M APY(1)=APY(1)+Y(I)ENDDOENDDOAPY(1)=APY(1)/D1DODO J=1,P U(J,1)=0.0DODO I=1,M U(J,1)=U(J,1)+V(J,I)ENDDOENDDO U(J,1)=U(J,1

19、)/D1ENDDOENDDOIFIF(Q1)THENTHEN D2=0.0 APY(2)=0.0DODO I=1,M G=Y(I)-APY(1)D2=D2+G*G APY(2)=APY(2)+(Y(I)*G*GENDDOENDDO APY(2)=APY(2)/D2 BY(2)=D2/D1DODO J=1,P U(J,2)=0.0DODO I=1,M G=Y(I)-APY(1)U(J,2)=U(J,2)+V(J,I)*GENDDOENDDO U(J,2)=U(J,2)/D2ENDDOENDDO D1=D2ENDIFENDIFDODO K=3,Q D2=0.0 APY(K)=0.0DODO J=1

20、,P U(J,K)=0.0ENDDOENDDODODO I=1,M G1=1.0 G2=Y(I)-APY(1)DODO J=3,K G=(Y(I)-APY(J-1)*G2-BY(J-1)*G1 G1=G2 G2=GENDDOENDDO D2=D2+G*G APY(K)=APY(K)+Y(I)*G*GDODO J=1,P U(J,K)=U(J,K)+V(J,I)*GENDDOENDDOENDDOENDDODODO J=1,P U(J,K)=U(J,K)/D2ENDDOENDDO APY(K)=APY(K)/D2 BY(K)=D2/D1 D1=D2ENDDOENDDOV(1,1)=1.0V(2,1

21、)=-APY(1)V(2,2)=1.0DODO I=1,PDODO J=1,Q A(I,J)=0.0ENDDOENDDOENDDOENDDODODO I=3,Q V(I,I)=V(I-1,I-1)V(I,I-1)=-APY(I-1)*V(I-1,I-1)+V(I-1,I-2)IFIF(I=4)THENTHENDODO K=I-2,2,-1 V(I,K)=-APY(I-1)*V(I-1,K)+V(I-1,K-1)-BY(I-1)*V(I-2,K)ENDDOENDDOENDIFENDIF V(I,1)=-APY(I-1)*V(I-1,1)-BY(I-1)*V(I-2,1)ENDDOENDDODOD

22、O I=1,PIFIF(I=1)THENTHEN T(1)=1.0 T1(1)=1.0ELSEIFELSEIF(I=2)THENTHEN T(1)=-APX(1)T(2)=1.0 T2(1)=T(1)T2(2)=T(2)ELSEELSE T(I)=T2(I-1)T(I-1)=-APX(I-1)*T2(I-1)+T2(I-2)IFIF(I=4)THENTHENDODO K=I-2,2,-1 T(K)=-APX(I-1)*T2(K)+T2(K-1)-BX(I-1)*T1(K)ENDDOENDDOENDIFENDIF T(1)=-APX(I-1)*T2(1)-BX(I-1)*T1(1)T2(I)=T

23、(I)DODO K=I-1,1,-1 T1(K)=T2(K)T2(K)=T(K)ENDDOENDDOENDIFENDIFDODO J=1,QDODO K=I,1,-1DODO L=J,1,-1 A(K,L)=A(K,L)+U(I,J)*T(K)*V(J,L)ENDDOENDDOENDDOENDDOENDDOENDDOENDDOENDDODT1=0.0DODO I=1,N X1=X(I)DODO J=1,M Y1=Y(J)X2=1.0 DD=0.0DODO K=1,P G=A(K,Q)DODO KK=Q-1,1,-1 G=G*Y1+A(K,KK)ENDDOENDDO G=G*X2 DD=DD+G

24、 X2=X2*X1ENDDOENDDO DT=DD-Z(I,J)DT1=DT1+DT*DTENDDOENDDOENDDOENDDORETURNRETURNENDEND三、计算结果三、计算结果数表(x,y,f(x,y):XYUXTYF(X,Y)0.000.5001.3450.2430.4465040260217E+000.000.5501.3220.2690.3246832410866E+000.000.6001.2990.2950.2101596387135E+000.000.6501.2770.3220.1030436532994E+000.000.7001.2550.3500.340191

25、2404030E-020.000.7501.2350.377-0.8873582055387E-010.000.8001.2150.406-0.1733716589358E+000.000.8501.1960.434-0.2505346523772E+000.000.9001.1770.463-0.3202764787154E+000.000.9501.1590.492-0.3826681228086E+000.001.0001.1420.521-0.4377957495101E+000.001.0501.1250.550-0.4857588795274E+000.001.1001.1090.

26、580-0.5266672522606E+000.001.1501.0930.609-0.5606385181100E+000.001.2001.0790.639-0.5877965434618E+000.001.2501.0640.669-0.6082697644452E+000.001.3000.001.3500.001.4000.001.4500.001.5000.080.5000.080.5500.080.6000.080.6500.080.7000.080.7500.080.8000.080.8500.080.9000.080.9500.081.0000.081.0500.081.1

27、000.081.1500.081.2000.081.2500.081.3000.081.3501.0500.6991.0370.7291.0240.7591.0110.7901.0000.8201.4150.2281.3910.2531.3680.2791.3460.3061.3250.3331.3040.3601.2840.3881.2650.4161.2460.4441.2290.4731.2110.5021.1940.5311.1780.5601.1630.5891.1480.6191.1330.6491.1190.6791.1060.709-0.6221894327719E+00-0.

28、6296883895048E+00-0.6308997691168E+00-0.6259561598452E+00-0.6149885531129E+000.6380152050067E+000.5066117295708E+000.3821763157102E+000.2648635435247E+000.1547802220057E+000.5199267859948E-01-0.4346806749373E-01-0.1316011021416E+00-0.2124310552929E+00-0.2860046192736E+00-0.3523860630008E+00-0.411655

29、3841229E+00-0.4639049136278E+00-0.5092367829793E+00-0.5477611307544E+00-0.5795943699592E+00-0.6048572233071E+00-0.6236734340898E+000.081.4001.0930.739-0.6361682760837E+000.081.4501.0800.769-0.6424676727883E+000.081.5001.0680.799-0.6426971176492E+000.160.5001.4830.2140.8400813766831E+000.160.5501.460

30、0.2390.6997641125864E+000.160.6000.160.6500.160.7000.160.7500.160.8000.160.8500.160.9000.160.9500.161.0000.161.0500.161.1000.161.1500.161.2000.161.2500.161.3000.161.3500.161.4000.161.4500.161.5000.240.5000.240.5500.240.6000.240.6501.4370.2641.4140.2901.3930.3161.3720.3431.3520.3701.3330.3981.3150.42

31、61.2970.4541.2790.4831.2620.5121.2460.5411.2310.5701.2160.6001.2010.6291.1870.6591.1740.6891.1610.7181.1480.7481.1360.7781.5510.2011.5270.2251.5040.2491.4820.2750.5660613612391E+000.4391716517306E+000.3192421655770E+000.2063761909359E+000.1006384973112E+000.2060692397577E-02-0.8935398464083E-01-0.17

32、36270461058E+00-0.2507999539520E+00-0.3209321867511E+00-0.3840977416274E+00-0.4403822530309E+00-0.4898811730659E+00-0.5326979436866E+00-0.5689418288271E+00-0.5987265475032E+00-0.6221686739756E+00-0.6393865589031E+00-0.6504993634875E+000.1051515083266E+010.9029273821703E+000.7605801844086E+000.624715

33、2461364E+000.240.7001.4600.3010.4955197895447E+000.240.7501.4390.3270.3731340505334E+000.240.8001.4190.3540.2576567336524E+000.240.8501.4000.3810.1491505238469E+000.240.9001.3810.4090.4764703563104E-010.240.9501.3630.437-0.4684940702542E-010.241.0000.241.0500.241.1000.241.1500.241.2000.241.2500.241.

34、3000.241.3500.241.4000.241.4500.241.5000.320.5000.320.5500.320.6000.320.6500.320.7000.320.7500.320.8000.320.8500.320.9000.320.9500.321.0000.321.0501.3460.4651.3290.4941.3130.5231.2970.5521.2820.5811.2670.6101.2530.6401.2400.6691.2270.6991.2140.7291.2020.7591.6170.1881.5930.2121.5700.2361.5470.2611.5

35、260.2861.5050.3121.4850.3391.4660.3651.4470.3931.4290.4201.4110.4481.3950.477-0.1343567558006E+00-0.2149132642759E+00-0.2885737112683E+00-0.3554064608815E+00-0.4154914249019E+00-0.4689182248284E+00-0.5157838183466E+00-0.5561910731630E+00-0.5902469746186E+00-0.6180615769184E+00-0.6397468483377E+000.1

36、271246761728E+010.1115001981949E+010.9646076509368E+000.8203474277775E+000.6824477015560E+000.5510852181777E+000.4263923704805E+000.3084629582599E+000.1973572015127E+000.9310554775472E-01-0.4285967343901E-02-0.9483383158541E-010.321.1001.3780.505-0.1785730041818E+000.321.1501.3630.534-0.255553892512

37、5E+000.321.2001.3470.563-0.3258401857129E+000.321.2501.3330.592-0.3895069599390E+000.321.3001.3190.621-0.4466381259700E+000.321.3501.3050.650-0.4973249497540E+000.321.4000.321.4500.321.5000.400.5000.400.5500.400.6000.400.6500.400.7000.400.7500.400.8000.400.8500.400.9000.400.9500.401.0000.401.0500.40

38、1.1000.401.1500.401.2000.401.2500.401.3000.401.3500.401.4000.401.4501.2920.6801.2790.7101.2670.7391.6810.1771.6580.1991.6340.2231.6120.2471.5910.2721.5700.2981.5500.3241.5300.3501.5120.3771.4930.4051.4760.4321.4590.4601.4430.4881.4270.5171.4120.5451.3970.5741.3830.6031.3690.6321.3560.6621.3430.691-0

39、.5416640906154E+00-0.5797565124179E+00-0.6117062922991E+000.1498321046591E+010.1334998574800E+010.1177125021283E+010.1025024091602E+010.8789600202794E+000.7391450727049E+000.6057448220594E+000.4788837887740E+000.3586506733933E+000.2451021277699E+000.1382683519668E+000.3815496242408E-01-0.55252856396

40、84E-01-0.1419870293898E+00-0.2220945000227E+00-0.2956352193106E+00-0.3626794362747E+00-0.4233061773066E+00-0.4776011207937E+00-0.5256554594343E+000.401.5001.3310.721-0.5675647513312E+000.480.5001.7450.1660.1731892811769E+010.480.5501.7210.1880.1562034590102E+010.480.6001.6980.2110.1397216878274E+010

41、.480.6501.6760.2350.1237801106740E+010.480.7000.480.7500.480.8000.480.8500.480.9000.480.9500.481.0000.481.0500.481.1000.481.1500.481.2000.481.2500.481.3000.481.3500.481.4000.481.4500.481.5000.560.5000.560.5500.560.6000.560.6500.560.7000.560.7501.6540.2591.6340.2841.6130.3101.5940.3361.5750.3631.5570

42、.3901.5390.4171.5220.4441.5060.4721.4900.5001.4750.5291.4600.5571.4460.5861.4320.6151.4190.6441.4060.6741.3940.7031.8080.1561.7840.1781.7610.2001.7390.2231.7170.2471.6960.2720.1084087585323E+010.9363227839056E+000.7947044362728E+000.6593871570749E+000.5304876715431E+000.4080885997266E+000.2922442321

43、330E+000.1829823414334E+000.8030848215007E-01-0.1579055318162E-01-0.1053445967045E+00-0.1883980538686E+00-0.2650070438339E+00-0.3352378374822E+00-0.3991645893867E+00-0.4568681687555E+00-0.5084349841214E+000.1971221944248E+010.1795329692430E+010.1624067146210E+010.1457830756268E+010.1296954768194E+01

44、0.1141718174733E+010.560.8001.6760.2970.9923495597011E+000.560.8501.6570.3230.8490326592663E+000.560.9001.6380.3490.7119114782397E+000.560.9501.6200.3750.5810940997652E+000.561.0001.6020.4020.4566585768956E+000.561.0501.5850.4290.3386546704747E+000.561.1000.561.1500.561.2000.561.2500.561.3000.561.35

45、00.561.4000.561.4500.561.5000.640.5000.640.5500.640.6000.640.6500.640.7000.640.7500.640.8000.640.8500.640.9000.640.9500.641.0000.641.0500.641.1000.641.1501.5680.4571.5520.4851.5370.5131.5220.5411.5080.5701.4940.5991.4810.6271.4680.6571.4550.6861.8700.1471.8460.1681.8230.1901.8010.2131.7790.2361.7580

46、.2601.7380.2851.7180.3101.6990.3361.6810.3621.6640.3881.6460.4151.6300.4431.6140.4700.2271082699420E+000.1220249611213E+000.2339219534309E-01-0.6881863698047E-01-0.1546492029899E+00-0.2341526439382E+00-0.3073911710726E+00-0.3744348739671E+00-0.4353605132598E+000.2215668024274E+010.2034201225610E+010

47、.1856955171354E+010.1684358335342E+010.1516776463456E+010.1354519097603E+010.1197844094442E+010.1046963013441E+010.9020461834743E+000.7632263799936E+000.6306048550718E+000.5042529674728E+000.3842166978607E+000.2705203064766E+000.641.2001.5980.4980.1631685158709E+000.641.2501.5840.5260.6214860470250E

48、-010.641.3001.5690.554-0.3256648583588E-010.641.3501.5550.583-0.1210165306376E+000.641.4001.5420.611-0.2032515079566E+000.641.4501.5290.640-0.2793303897633E+000.641.5000.720.5000.720.5500.720.6000.720.6500.720.7000.720.7500.720.8000.720.8500.720.9000.720.9500.721.0000.721.0500.721.1000.721.1500.721.

49、2000.721.2500.721.3000.721.3500.721.4000.721.4500.721.5001.5160.6691.9310.1391.9070.1591.8840.1811.8620.2031.8400.2261.8190.2491.7990.2731.7790.2981.7600.3231.7420.3491.7240.3751.7070.4021.6910.4291.6750.4561.6590.4841.6440.5111.6300.5391.6160.5681.6020.5961.5890.6251.5760.653-0.3493199222056E+000.2

50、464684390094E+010.2278059075084E+010.2095251278336E+010.1916718301640E+010.1742854737447E+010.1573998475956E+010.1410434829219E+010.1252401695837E+010.1100094482886E+010.9536697141176E+000.8132510552324E+000.6789308703755E+000.5507747946397E+000.4288254602727E+000.3131046786931E+000.2036155380751E+0