原创力文档- 1、本文档共8页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
数值分析上机题作业
电器学院 交通信息工程及控制 罗宁20050290010110 第一章
源程序：
REALS
S=1
DO 1=1,71
S=S*29/I
END DO
WRITE(*,*) MY=M,S
END
结果：Y=79.5416
第三章
源程序：
REAL X,Y,N,K,DX,NEWTON,NE,S,XK
DIMENSION X(6),Y(6),NE(6,6),NEWTON(21)
DATA X/0.2,0.24,0.28,0.32,0.36,0.4/
DATA Y/0.1987,0.2377,02764,0.3146,0.3523,0.3894/
!计算差商表：
DX=0.04
DO 1=1,6
NE(I,1)=Y(I)
END DO
DO J =2,6
DO I=J,6
NE(IJ)=(NE(I,J-1)-NE(I-M-1))/(DX*(J-1))
END DO
END DO
DO 1=1,21
XK=0.2+(I-l)*0.01
NEWTON(I)=NE(1,1)
S=1
DO J=2,6
S=S*(XK-X(J-1)) NEWTON(I)=NEWTON(r)+S*NE(J,J)
END DO
END DO
WRITE(*,*) UNEWTON=H,NEWTON
END
NEWTON=0.1987000.2084510.2182090.2279
NEWTON=0.198700
0.208451
0.218209
0.2279
0.247417
0.257108
0.266770
0.276400
0.285998
0.305097
0.314600
0.324072
0.333513
0.342923
0.352300
0.370943
0.380199
0.389400
结果:
0.237700
0.2955640.361642
第四章
第四章
源程序:
REALA,B,ANS,D
DIMENSION A(4,4),B(4),ANS(4),D(4》)
DATA A/4197,6.8,88.6,1.45305,71.3,76.4,5.9,-206,
-47.4,-10.8,6.13,-840,52^0236.5/
DATA B/136,ll?5,25.7,6.6/
CALL LGAUS( A,B,4,D)
CALL GAUSQ(D,4,ANS)
WRITE(笃? ,ANS=,,ANS END
-1.684955E-02结果：ANS=3.226190E-02 0.322386 0.291290
-1.684955E-02
第五章 源程序：
RE AL A l ,B 1 ,A2,B2,E1 ,E2
Al=0;A2=0
B1=3.141592654/2;B2=1
E1=0.000001;E2=0.001
CALL SIN1PB1(A1,B1,E1,S)
WRITE?*)，第一题 simpson:
CALL SIMPB2(A2,B2,E1,S)
WRITE?*)，第二题 simpson:\S
CALL ROMBERG1(A1,B1,E2,R)
WRITER,*)，第一题 ROMBERG ,R
CALL ROMBERG2(A2,B2,E2,R)
WRITE(笃*)，第二题 ROMBERG ,R END
SUBROUTINE S1MPB1(A,B,E,S)
F(X)=SIN(X)/X
H=(B-A)/2
S2=0;N=l
SO=1+F(B)
S1=F(A+H)
S=(S0+4.0*Sl)*H/3.0
60 N=2*N;H=H/2.0
S2=S2+S1
SI =0.0
X=A+H
DO 1=1,N
S1=S1+F(X)
X二X+H+H
END DO
S2N=(S0+2?0*S2+4?0*Sl 严 H/3.0
IF(ABS(S2N-S)E)THEN S=S2N
GOTO 60
END EF
END SUBROUTINE
SUBROUTINE SIMPB2(A,B,E,S)
F(X)=LOG(1+X)/X
H=(B-A)/2
S2=0;N=l
SO=1+F(B)
S1=F(A+H)
S=(S0+4.0*Sl)*H/3.0
60 N=2*N;H=H/2.0
S2=S2+S1
Sl=0.0
X=A+H
DO 1=1,N
S1=S1+F(X)
X二X+H+H
END DO
S2N=(S0+2.0*S2+4.0*Sl 严 H/3.0
IF(ABS(S2N-S)E)THEN
S=S2N
GOTO 60
END IF
END SUBROUTINE
SUBROUTINE ROMBERG1(A,B,E,R)
F(X)=SIN(X)/X
REAL S,DS,T
INTEGER K,U
DIMENSION S(100,100)
S(l,l)=
文档评论（0）