資源描述：
《數(shù)值分析上機題作業(yè).doc》由會員上傳分享，免費在線閱讀，更多相關(guān)內(nèi)容在工程資料-天天文庫。
1、數(shù)值分析上機題作業(yè)電器學(xué)院交通信息工程及控制羅寧20050290010110第一章源程序：REALSS=1DO1=1,71S=S*29/IENDDOWRITE(*,*)MY=M,SEND結(jié)果：Y=79.5416第三章源程序：REALX,Y,N,K,DX,NEWTON,NE,S,XKDIMENSIONX(6),Y(6),NE(6,6),NEWTON(21)DATAX/0.2,0.24,0.28,0.32,0.36,0.4/DATAY/0.1987,0.2377,02764,0.3146,0.3523,0.3894/!計算差商表：DX=0.0
2、4DO1=1,6NE(I,1)=Y(I)ENDDODOJ=2,6DOI=J,6NE(IJ)=(NE(I,J-1)-NE(I-M-1))/(DX*(J-1))ENDDOENDDODO1=1,21XK=0.2+(I-l)*0.01NEWTON(I)=NE(1,1)S=1DOJ=2,6S=S*(XK-X(J-1))NEWTON(I)=NEWTON(r)+S*NE(J,J)ENDDOENDDOWRITE(*,*)UNEWTON=H,NEWTONENDNEWTON=0.1987000.2084510.2182090.22790.2474170.25
3、71080.2667700.2764000.2859980.3050970.3146000.3240720.3335130.3429230.3523000.3709430.3801990.389400結(jié)果:0.2377000.2955640.361642第四章源程序:REALA,B,ANS,DDIMENSIONA(4,4),B(4),ANS(4),D(4》)DATAA/4197,6.8,88.6,1.45305,71.3,76.4,5.9,-206,&&-47.4,-10.8,6.13,-840,52^0236.5/DATAB/136,l
4、l?5,25.7,6.6/CALLLGAUS(A,B,4,D)CALLGAUSQ(D,4,ANS)WRITE(篤?,ANS=,,ANSEND-1.684955E-02結(jié)果：ANS=3.226190E-020.3223860.291290第五章源程序：REALAl,B1,A2,B2,E1,E2Al=0;A2=0B1=3.141592654/2;B2=1E1=0.000001;E2=0.001CALLSIN1PB1(A1,B1,E1,S)WRITE?*)，第一題simpson:CALLSIMPB2(A2,B2,E1,S)WRITE?*)，第二
5、題simpson:SCALLROMBERG1(A1,B1,E2,R)WRITER,*)，第一題ROMBERG,RCALLROMBERG2(A2,B2,E2,R)WRITE(篤*)，第二題ROMBERG,RENDSUBROUTINES1MPB1(A,B,E,S)F(X)=SIN(X)/XH=(B-A)/2S2=0;N=lSO=1+F(B)S1=F(A+H)S=(S0+4.0*Sl)*H/3.060N=2*N;H=H/2.0S2=S2+S1SI=0.0X=A+HDO1=1,NS1=S1+F(X)X二X+H+HENDDOS2N=(S0+2?0
6、*S2+4?0*Sl嚴H/3.0IF(ABS(S2N-S)>E)THENS=S2NGOTO60ENDEFENDSUBROUTINESUBROUTINESIMPB2(A,B,E,S)F(X)=LOG(1+X)/XH=(B-A)/2S2=0;N=lSO=1+F(B)S1=F(A+H)S=(S0+4.0*Sl)*H/3.060N=2*N;H=H/2.0S2=S2+S1Sl=0.0X=A+HDO1=1,NS1=S1+F(X)X二X+H+HENDDOS2N=(S0+2.0*S2+4.0*Sl嚴H/3.0IF(ABS(S2N-S)>E)THENS=S
7、2NGOTO60ENDIFENDSUBROUTINESUBROUTINEROMBERG1(A,B,E,R)F(X)=SIN(X)/XREALS,DS,TINTEGERK,UDIMENSIONS(100,100)S(l,l)=(B-A)*(F(B)+l)/2K=1T=0100K=K+1S(K,1)=1/2*S((K-1),1)U=2**(K-2)-lDOJ=0,UT=T+F(A+(2*J+1)*(B-A)/(2**(K-1)))ENDDOS(K,1)=S(K,1)+T*(B?A)/2**(K?1)DOJ=2,KS((K-J+l)J)=(4*
8、*(J-l)*S((K-J+2),(J4))+S((K-J+l),(J-l)))/&&(4**(J-1)-1)ENDDODS=S(1,K)-S(1J<-1)IF(ABS(DS).GT.E)TH
