D=AB(S2-S1)/S2) END IF IF(D.LT.EPS)GOTO50 0 N-N+N H=0.5*H T1=T2 S1=S2 G0T010 RETURN END SUBROUTINE LGRG2(X,Y.N,T.Z) (5分) DIMENSION X(N).Y(N) DOUBLE PRECISION X.Y.TZ.S Z=0.0 IF(N.LE.0)RETURN IF(N.EQ.1)THEN Z=Y1) RETURN ENDIF IF(N.EQ.2)THEN Z=Y()(T-X(2)-Y(2)(T-X(D(X(1)-X(2)) RETURN END IF IF(T.LE.X(2))THEN K M=3 ELSE IF(T.GE.X(N-1))THEN K=N- M-N M=N 10 IF(IABS(K-M).NE.1)THEN L+MW2 IF(T.LT.X(L))THEN ELS S(T-X(K)).LT.ABS(T-X(M)))THEN ELSEK-I M+I END IE 8 D=ABS((S2-S1)/S2) END IF IF(D.LT.EPS)GOTO 50 40 N=N+N H=0.5*H T1=T2 S1=S2 GOTO 10 50 RETURN END SUBROUTINE LGRG2(X,Y,N,T,Z) （5 分） DIMENSION X(N),Y(N) DOUBLE PRECISION X,Y,T,Z,S Z=0.0 IF(N.LE.0)RETURN IF(N.EQ.1)THEN Z=Y(1) RETURN END IF IF(N.EQ.2)THEN Z=(Y(1)*(T-X(2))-Y(2)*(T-X(1)))/(X(1)-X(2)) RETURN END IF IF(T.LE.X(2))THEN K=1 M=3 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