ELSE IF(T.GE.X(N-1))THEN K=N.2 M-N ELSE K=1 M=N 10 IF(IABS(K-M).NE.1)THEN L=K+M0/2 IF(T.LTX(L))THEN ELSE K=L EnDIf G0T010 END IF IF(ABS(T-X(K))LT.ABS(T-X(M)))THEN K=K-1 ELSE M-M+1 END IF END IF Z=0.0 DO30I=K.M S=10 DO20J=K.M IF(J.NE.I)YTHEN S=S*(T-X()V(X(1-X()) END IF CONTINUE Z=Z+S*Y(I) 30 CONTINUE RETURN END 4.运行结果。(20分) A -1479000221D+01 B= 2303009515D+0 .9756013465D+00 n 2303009515D+01 k= .2278653895D+00 动力学方程式为： vs、de 2303009515D+01 2278653895D+00 ELSE IF(T.GE.X(N-1))THEN K=N-2 M=N ELSE K=1 M=N 10 IF(IABS(K-M).NE.1)THEN L=(K+M)/2 IF(T.LT.X(L))THEN M=L ELSE K=L ENDIF GOTO 10 END IF IF(ABS(T-X(K)).LT.ABS(T-X(M))) THEN K=K-1 ELSE M=M+1 END IF END IF Z=0.0 DO 30 I=K,M S=1.0 DO 20 J=K,M IF(J.NE.I)THEN S=S*(T-X(J))/(X(I)-X(J)) END IF 20 CONTINUE Z=Z+S*Y(I) 30 CONTINUE RETURN END 4．运行结果。（20 分） A= -.1479000221D+01 B= .2303009515D+01 R= .9756013465D+00 n= .2303009515D+01 k= .2278653895D+00 动力学方程式为： n kc t c = − = d d ν （n= .2303009515D+01 k= .2278653895D+00）