COMMON /FIT/X,Y,SIGMAY,NPTS,MODE,A,SIGMAA,B,SIGMAB,R
      DOUBLE PRECISION X,Y,SIGMAY,A,SIGMAA,B,SIGMAB,R
      INTEGER NPTS,MODE,I
      DIMENSION X(50),Y(50),SIGMAY(50)
1     TYPE 100
      ACCEPT 101,NPTS
      IF (NPTS.EQ.0) GOTO 999
      TYPE 102
      ACCEPT 103,MODE
      TYPE 104
      DO 10 I=1,NPTS
10    ACCEPT 105,X(I),Y(I),SIGMAY(I)
      CALL LINFIT
      TYPE 200,A,SIGMAA
      TYPE 201,B,SIGMAB
      TYPE 202,R
      TYPE 203
      GOTO 1
999   STOP
100   FORMAT(' NUMBER OF POINTS='$)
101   FORMAT(I)
102   FORMAT(' MODE='$)
103   FORMAT(I)
104   FORMAT('0INPUT DATA: INDEPENDENT/DEPENDENT[/SIGMA(DEPENDENT)]'/)
105   FORMAT(3D)
200   FORMAT(' A= ',1PD20.5,9X,'SIGMAA=',1PD20.5)
201   FORMAT(' B= ',1PD20.5,9X,'SIGMAB=',1PD20.5)
202   FORMAT(' CORRELATION= ',F15.7)
203   FORMAT('1')
      END