Re: 選課
本篇為計程作業9低調版
1. Lagrange Interpolation:
PROGRAM Lagrange
REAL(8) XX(1:7),YY(1:7),X,Y,Z
INTEGER I,J
DATA XX/3,5,6,9,13,15,16/, YY/16,32,23,45,56,12,9/
READ(*,*)X
IF (X==0) STOP ! Eliminate this line
DO 777 I=1,7
Z = YY(I)
DO 888 J=1,7
IF (I.NE.J) THEN
Z = Z*(X-XX(J))/(XX(I)-XX(J))
END IF
888 CONTINUE
777 Y=Y+Z
WRITE(*,*)Y
END
2. 解聯立：
PROGRAM equations
INTEGER I,J,K
REAL N(1:4,1:5), TEMP
DATA N/1,2,4,-3,3,1,2,1,1,0,2,3,5,3,1,2,4,5,11,3/
DO 70 I=1,4
DO 70 J=1,4
IF(I.NE.J)THEN
60 TEMP=-(N(J,I)/N(I,I))
DO 90 K=1,5
90 N(J,K)=N(J,K)+N(I,K)*(TEMP)
END IF
70 CONTINUE
DO 55 I=1,4
55 WRITE(*,*)N(I,5)/N(I,I)
END
※ 編輯: snowangel 來自: 220.135.2.143 (12/10 17:52)
討論串 (同標題文章)