perm filename PH3.F4[NET,GUE]1 blob
sn#026509 filedate 1973-02-27 generic text, type T, neo UTF8
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