<div id='wikitext'>
<pre><b> File: DMEAN.FT of Tape: Various/ETH/f4</b></pre>
<pre><b>(Source file text) </b></pre>
<hr />
<pre>
      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 
</pre>
<hr/>
</div>