File: DFMU.FT of Tape: Various/ETH/prog1
(Source file text)
FUNCTION FMU(MU) 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/ IF,JR,INIT,LS COMMON/ACCUR/ DDMU1,DDMU2,DCHI,ERRMU C 1000 FORMAT(5X'MU ='E17.10' FMU ='E17.10' AMU ='E17.10) REAL MU LOGICAL ERRMU AMU0=MU DO 10 J=1,MR JR=J IF(LS.EQ.0) CALL LSTAU IF(LS.EQ.1) CALL LSAFTAU 10 CONTINUE CALL LOGAMP CALL LINREG(MR) FMU=AMU-MU C WRITE(4,1000) MU,FMU,AMU RETURN END