5 * Revision 1.1.1.1 1996/04/01 15:02:28 mclareni
10 SUBROUTINE LFITW(X,Y,W,L,KEY,A,B,E)
12 C TO PERFORM A WEIGHTED STRAIGHT LINE FIT
14 C FOR FORMULAE USED SEE MENZEL, FORMULAS OF PHYSICS P.116
16 C FIT IS OF Y=AX+B , WITH S**2 ESTIMATOR E. WEIGHTS ARE IN W.
17 C IF KEY=0, POINTS WITH Y=0 ARE IGNORED
20 DIMENSION X(1),Y(1),W(1)
32 IF(Y(J).EQ.0..AND.KEY.EQ.0) GO TO 2
46 A=(W2XY-W2X*W2Y/W2)/(W2X2-W2X**2/W2)
49 E=(W2Y2-W2Y**2/W2-(W2XY-W2X*W2Y/W2)**2/(W2X2-W2X**2/W2))/(ICNT-2)