This commit was generated by cvs2svn to compensate for changes in r2,
[u/mrichter/AliRoot.git] / 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