C FIT LINE TO EQUALLY SPACED POINTS... C Y=BX+A SUBROUTINE LINFIT(X,Y,N,A,B,DEL) DIMENSION X(1),Y(1) REAL*8 X,Y,A,B,XY,X2,SX,SY,SY2,EN,DEL INTEGER*2 N XY=0. X2=0. SX=0. SY=0. SY2=0. EN=N DO 10 I=1,N XY=XY+X(I)*Y(I) C FORM SUMS NEEDED TO FIT LINE... X2=X2+X(I)*X(I) SX=SX+X(I) SY=SY+Y(I) SY2=SY2+Y(I)*Y(I) 10 CONTINUE C NOW GET SLOPE B=(XY-(SX*SY)/EN)/(SX2-(SX*SX)/EN) C THEN INTERCEPT A=(SY/EN)-B*(SX/EN) DEL=DSQRT((SY2-(A*SY+B*XY))/EN) RETURN END