File: DMEAN.FT of Tape: Various/ETH/prog1
(Source file text)
SUBROUTINE MEAN C COMMON/INI/ ANF(10,10),DANF(10,10),PH(10,10),DPH(10,10) $,AF(10,10),DAF(10,10),ETA(10,10),DETA(10,10),EMU(10,10), $T(10),DT(10),ANFM(10),DANFM(10) CDC $,F(10,512) COMMON/TESPAR/ A0M,DA0M,AMU0,AF0(10),ETA0(10),A0(10) $,TDJ(10),DTDJ(10),TDM,DTDM COMMON/PARAM/IDENT,HANG,HI,HF,DK,NP,PI,SQHI,DKHI,I0DH,H0AL COMMON/STEER/MR,ADENT(10),MF,FREQ(10),I0 COMMON /LINE/ A,B,DA,DB,SA,SB,SAB,TD,AMU,RR,RMSD,CHI COMMON /REGRES/ XL(10),YL(10),DY(10) COMMON/FRERUN/ I,JR,INIT,LS COMMON/ACCUR/ DDMU1,DDMU2,DCHI,ERRMU C SP=0. SA=0. SB=0. DO 10 J=1,MR TAU=ETA0(J)*DK*NP EX=EXP(-TAU) A0(J)=AF0(J)*TAU*SQHI*EXP(ETA0(J)/HI)/(T(J)*(1.-EX)) PDETA0=(1.-TAU*EX/(1.-EX))/ETA0(J)+1./HI A0P=A0(J)**2*((DAF(I,J)/AF0(J))**2+(DETA(I,J)*PDETA0)**2) A0P=1./A0P SA=SA+A0(J)*A0P SB=SB+A0(J)**2*A0P SP=SP+A0P 10 CONTINUE A0M=SA/SP DA0M=SQRT((SB/SP-A0M**2)/(MR-1)) RETURN END