MINICERN/mathlib/gen/d/errorf.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1996/04/01 15:02:21  mclareni
   6 * Mathlib gen
   7 *
   8 *
   9 #include "gen/pilot.h"
  10       SUBROUTINE ERRORF (M)
  11 C-----FINAL ERROR ANALYSIS FOR CHISQ.
  12 #include "d510si.inc"
  13 #include "d510ui.inc"
  14 #include "d510uo.inc"
  15       EQUIVALENCE(NED(1),NS),(NED(2),NP)
  16       WRITE(6,3)
  17       IF (NP.GT.7) WRITE(6,4)
  18       DO 2 J=1,NS
  19       K1=NP*(J-1)+1
  20       FUN=EXDA(K1)
  21       SIG=EXDA(K1+1)
  22       NX=NP-2
  23       DO 1 K=1,NX
  24       K2=K1+1+K
  25       X(K)=EXDA(K2)
  26  1    CONTINUE
  27 C-----CALCULATE DY/DA"S FOR SCAL
  28       CALL ARITHM (Y)
  29       CALL SCAL (M,ER)
  30       XI=(FUN-Y)/SIG
  31       XI=0.5*XI*XI
  32       ER=SQRT(ER)
  33       ERROR(J)=ER
  34       WRITE(6,5) J,Y,ER,XI,(X(I),I=1,NX)
  35  2    CONTINUE
  36       RETURN
  37 C
  38  3    FORMAT('1','ERROR ANALYSIS FOR EACH DATA POINT USING FINAL PARA',
  39      +'METER VALUES (CHISQ ONLY)'////11X,'POINT',8X,'FITTED',9X,
  40      +'Y STANDARD',6X,'CONTRIBUTION',4X,'X CO-ORDINATES'/12X,'NO.',8X,
  41      +'Y  VALUE',8X,'DEVIATION',7X,'TO CHISQ/2',6X,'OF DATA POINT'//)
  42  4    FORMAT('+',86X,'(PRINTED IN ROW ORDER)')
  43  5    FORMAT(11X,I3,2X,7(5X,E12.5),(/67X,4(5X,E12.5)))
  44       END