5 * Revision 1.1.1.1 1996/04/01 15:02:28 mclareni
10 #if defined(CERNLIB_FORTRAN)||!defined(CERNLIB_CDC)
11 SUBROUTINE LFIT(X,Y,L,KEY,A,B,E)
13 C TO FIT A STRAIGHT LINE Y=A*X+B TO L POINTS WITH ERROR E
14 C SEE MENZEL , FORMULAS OF PHYSICS P.116
15 C POINTS WITH Y=0 ARE IGNOERD IF KEY=0
29 IF(Y(J).EQ.0..AND.KEY.EQ.0) GO TO 10
34 IF(COUNT.LE.1.) GO TO 25
38 IF(Y(J).EQ.0..AND.KEY.EQ.0) GO TO 20
41 SUMXY=SUMXY+SCARTX *SCARTY
42 SUMXX=SUMXX+SCARTX *SCARTX
43 SUMYY=SUMYY+SCARTY *SCARTY
47 IF(SUMXX.EQ.0.) GO TO 25
50 IF(COUNT.LT.3.) GO TO 101
51 E=(SUMYY-SUMXY*A )/(COUNT-2.0)