IFON.EO.DGOTO 40 30 D=ABS(S2-S1) IF(ABS(S2).GT.1.0)THEN D=ABS((S2-S1)S2) END IF IF(D.LT.EPS)GOTO 50 40 n=nin H=05*H T172 s1=S2 G0T010 50 Return END SUBROUTINE LGRG2(X.Y.N.T.Z) DIMENSION X(N)Y(N) DOUBLE PRECISION X.Y.T.Z.S 7=00 IF(N.LE.O)RETURN IF(N.EQ.1)THEN Z=Y1) RETURN END IF IF(N.EQ.2)THEN ZY1)*(TX2Y(2)*(T-X1)MX1)X2) RETURN END IF IF(T.LE.X(2))THEN K=1 M=3 ELSE IF(T.GE.X(N-1)YTHEN K=N-2 M=N ELSE 10 IF(IABS(K-M).NE.1)THEN L-K+M0/2 IF(TLTXOLDTHEN M-I ELSE K=L ENDIF G0T010IF(N.EQ.1)GOTO 40 30 D=ABS(S2-S1) IF(ABS(S2).GT.1.0)THEN 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) 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