SS=EXP(A) WRITE(6 IDA BRS SS 11 D1810/5XB=D18.10,5XR=,D18.10& D1 810,5X,k=D18.100 END SUBROUTINE CF(N,H,X,Y,Z) DIMENSION X(N),Y(N).Z(N) DOUBLEPRECISIONX.Y.Z.H.T.YA.YB T=X(1)+H CALL LGRG2(X.YN.TYAY T=X(1+2*H CALL LGRGZ(X YNTYBY Z(I)-(4-YA-3*Y(I-YBYH) D0101-2,N- T=X(I-H CALLLGRG2X YNTYAY T-XOIH CALLLGRG2(X.Y.N.T.YB) Z(D)-(YB-YA)(2*H 10 CONTINUE T-XONH CALL LGRG2CX.Y.N.TYAY T=X(N)2*H CALLLGRG2(X.Y.N.T.YB) Z(N(3*Y(N)4*YA+YB)(2*H) RETURN END SUBROUTINE PK(N.X.Y.A,B,R) DOUBLEPRECISIONX.Y.A.B.R DIMENSION X(N)Y(N) SX=0.0 8V=00 SXX=00 SYY-0.0 SXY-0.0 DO 10I=LN X1=X① Y1-Y0) SX-SX+XI SY-SY+Y1 SXX=SXX+XI*X SYY=SYY+Y1*YI SS=EXP(A) WRITE(6,11)A,B,R,S,SS 11 FORMAT(/5X,'A=',D18.10//5X,'B=',D18.10, 5X,'R=',D18.10// & 5X,'n=',D18.10,5X,'k=',D18.10/) END SUBROUTINE CF(N,H,X,Y,Z) DIMENSION X(N),Y(N),Z(N) DOUBLE PRECISION X,Y,Z,H,T,YA,YB T=X(1)+H CALL LGRG2(X,Y,N,T,YA) T=X(1)+2*H CALL LGRG2(X,Y,N,T,YB) Z(1)=(4*YA-3*Y(1)-YB)/(2*H) DO 10 I=2,N-1 T=X(I)-H CALL LGRG2(X,Y,N,T,YA) T=X(I)+H CALL LGRG2(X,Y,N,T,YB) Z(I)=(YB-YA)/(2*H) 10 CONTINUE T=X(N)-H CALL LGRG2(X,Y,N,T,YA) T=X(N)-2*H CALL LGRG2(X,Y,N,T,YB) Z(N)=(3*Y(N)-4*YA+YB)/(2*H) RETURN END SUBROUTINE PK(N,X,Y,A,B,R) DOUBLE PRECISION X,Y,A,B,R DIMENSION X(N),Y(N) SX=0.0 SY=0.0 SXX=0.0 SYY=0.0 SXY=0.0 DO 10 I=1,N X1=X(I) Y1=Y(I) SX=SX+X1 SY=SY+Y1 SXX=SXX+X1*X1 SYY=SYY+Y1*Y1