This commit was generated by cvs2svn to compensate for changes in r2,
[u/mrichter/AliRoot.git] / MINICERN / mathlib / gen / d / errorf.F
CommitLineData
fe4da5cc 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)
11C-----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
27C-----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
37C
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