-  

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 [ 81 ] 82 83

f *

0 1 +

t

. 8.

-

/1

- -

, 7

-

-

. 9. -

. .

..=

(2>

1= - - ; - ( ); 1/ - , - .

(. . 8)

-. ( Vk, U . .) ;



SUBROUTINE AWttVl.VZ.T.EM.R ,... W, ? . ¥ 1, M2 Ml 1 ,M6 ,M7 , , n ,TE , TBI I Dl,M8,H9y MlOJl cOMMON/fll /1 ,P2/J\2/Pl-PJ,P7,pe, P9/J1J/Z ,Wl We . H . AW, HG .RA UrnENSlON S3{30)CRt*l K.N

.oOubUE precision ., ET .fnMEnMl.WS.WSl,WAF.WMf 1,,,,1/, HJPiMl, l,WA,WAl.Tn.Tni.Wl.w2,vO, Nl>Vl.V?,S,yol<VN2.TI,EnT 51.S3, 2uAT,vl,V2,H.Pl.P2.P3,PiiP5-P6-P7.S.P -WH ..flW ,VF,vFl,Z2.V3,V4,T6 J,Dl,PM,PMT

r 1-2= 2 M1 J-M3<jiM1 HKlioMf 5rw5<,M16iM*e. 17= 7 tRUS . iSi HUiSofM . -/ ll

pl = Vl-V20f = Vl*V2 P3 = Vl.v2cP<, = vl. 2cP5=V2 2 PeP*-P5<>P7=P**P3 p8=P* P5 P9=P6..2CZ=K EM.Pl/r/P3 Z2:(K.EM> .2/n/R WH:2 , )>0 = 2 , n 1 =-/ .eWJ = . n 10W2-2.>.l.p) ET = 05eRT(R-*n* 2 TBTK.N.fl ./,ei/a)6lFDl.<(iHj2 2,}

.2 = . 7/ /< ....- .2) = 6/2 . PM = UAM><2/ ./RA EnMl=K fiif.EH/(2,.n 2 v2 *H..2)cUftHl = SnM141G/2.

wS = Z. 2/(2. n. . 2 R *1 3)<>WSlsK il2.EM. -2/ Zil . . R V2 . 2-! 3 ) =2 2 (1., 2*1. 2)/(-,2. (2,**..2- -2> 3

WA=S3 -Wflf.(loWAF 1=1 X.EM) . 2 . / (.. 5 * <.) OWA 1 = S WPf 1

=:2/1.. 2 - 2)=£/1

pRlNT<.,Ml,ET,tnM,enMl.UAM,UftMl> PM.WS .WS l .WAF ,WAF 1 .WA.WAl , To .T3

=.5/cnf=.5/

CALL ( Ml 2 ,Ml , Ml . 1 cC (LL I M 1 5 . Ml 6 , M 1 ? . 1 ) To 17 ) H=(W2-W1)/M. 2N=M*1

3 dO1J = 1J0C = 2-?4JOM=1W2-W1) /01 =/ Vl=0.t)oc.v2 = C.D0oVO=V-i.wl]cVNj=vf(W2) 006I = l.KlcVl;iH.i H.VF(Wi (i l.ll Hl

V2* 2, VF (Wl 2 1 )

5HH.<V0-SINl)*yj y2)3.OS3 J)=S 1 F J . EO. 1 iGOTO 1

iFiDaBS(S3CJ)-S3(J-1)1.6E.CA6S(S3(J1/S4)G0t01

WUF = Z2.SeWA= s3 .WAf OOPRINT7,M1,J SJ(J-1).S3(J),WAF,WAeboTO 15 1 CONTINUE 5 dOBK2Me.M9.MlCCTI=.125 2/

D016J = 1 ,3COM=2 2 . J H = (H2- l j/MOtCiiM/Z

VJ=0.0OCV =0.D0 VOl=VFi(Wl)uVN2=Vfl(W2)

00SiKJ=l,Ki.OV3 = VJ.<i. H VFl(wl.l2.K3-ll Hl S 2 . HVF 1 Wi < 1

SlslM lV0l VN2) V3 i .)/l.eS3lJ) = Sl

IF(J. , 1 jcOl0 16

1F(04BS(S3(J).S3(J-1)>.GE,DabScS3<J)/S )) OtO 16

= .5 1/= = 2/*./

pRINT 1.Tl,J.S3(J-n.S3ij)... 8 ;le cONTlNyE e CONTINUE .17 CONTINUE

4 FORMAT ( . ,l2/l2x,ET=,Dl2.5, . = . DI 2 . 5 . X . E 1 = . 1 D12 . 5/2ex. UAM5 , DI 2 . 5. >: . ! .Dl2.5.2x. PM= D12.5/12XWS= .01 22.?.X>WSl=,Dl2.5.2x.WAF=,0l2.5.X,WAFl=-.01Z.3/*7)l,-wA=.Dl2.3 3,2X,WAji.,Dj2.5/l2.T0=.Ol2,5,2X,-T3=.iDl2 .-SZ/leX , .. 4,ex, .- ,9X.EOMT , 10> . AT , 1lx, < pMT j

7 FORM ATII4X, - .I2/12X,J=.I3,X,S3<J-l)=.Dl2.5,x.S3(Jl = . l,Dl2.5/lBX,W F=..Dl2.5.x.Vs=.D12.5/i7x.TJ.8x.J.5x.S3(J-l) 2,eA.S3(J),j0X,.EnT,iiii,UAT.,j0l, BMT. 1

10 FOR Ti nx.ojz.5,.13. i2x,012.511

pETuRN END

SUBROUTINE (Ml,M2.M3,T1I i OMMON/i/Z.Wl.W2,WH.flW,HGA OOueLE PRECISION Z . X , T, Tl . Wl , WH flW. f ,. , E OMT , uA T . pMT DOlIrMl,:M2.M3 !. I 0 1.2.

EnT=7 (DSlN(W2 T )-OSlN(Wi,T ) ,/XeEDnT-2. Z DsIKJ( AW.T/2 )/X e0MT=2..Z/X UAT=H6.EnT/2. pMT=:UAT..2/t /.

1 pRlN72,T,EnT,E0nT,6OMT,UAT PMT

Z f0RmaT(I1X.6(X.Ol2.51)

rEturncend

DOuBlE PRECISION FUNCTION VF(X)

C0M CN/JI2/PJ , P3, P7,Pe,P cp0tiBLE pREClblON P J , P3 . P 7 , Pe , P9 , X , X J Xl=X. 20VF:((Pl.(l.-Pj.Xl) )..2♦X1 P9)/( (J.♦P7 >1♦X1 2 2) rETURN END

DOUBLE PRECISION PWCTION VFJIX) COMMON//1,,7.8.9/l/Tl,P2

DOUBLE PRECISION n .X . XI ,PJ,p2,P3.P7,Pe.P9CXISX 2

vFlsDSeRT(iPj.(i .p3.,Xl)) 2,Xl.F l.DCDSX TI-DATAN(X P2/H.-P3 xl I 1 ) 1/ 1 . .P7 Xi Pe.X 1 2) RETURNOEND



PROGRAM.

DATA TL.Tl.TOiTA.UM.SC.n.kzn.DUT a/4.E-7,1.2F-6,S-E-7,3.E-7,l.,l.,3-i4l59 3,Z.506 28,4.E-7/ T01=T0/2.oT02=T0/10.CA=(TO-TA)/2. TO=THOTH=DUT-T01 Sl=0, PRINT 2,

rOl K=l,lfl 0SUSX+.4 <X1.=-Z.S1 PH=1./SG/K2n

PK = PM EXP(-(Xl*.2)/2./(SG..2)> THl=TOa.PK K 1=:*01 PRINT 2,

8 CALL STROB(T1,TL,TH1,A,TA,A.UM,PUT1,0.0,0) PRINT 2,

DO 3 Kl=l,10 =0.01 1 1= 1-EZ PRINT 4,1,1

3 CALL STROB(Tl,TL.TH.A,,A,UM,OUT,0,TO,) PRINT 2, DUT2 = DUT 0tH2=02-T01

DO 5 K2=l,10 CDTK=0.Oi K2 PRINT 4,K2,DTK

5 Call StRob(TI.tl,tH2,a,,a.uM,2,dtk,,O)

PRINT 2, *SJ=0 D06 K3 = l,10 t,S3 = S3+2.JE-2 TL1-TL*S3 TL OPRJNT 7,K3,TL1

6 CALL STROB(1,TLl.TH.A.,,,,0,0,0) 2 F0RMATfy/5X, /)

4 F0RMAT(3X,3H N=,I3,6H =,Ell.S)

7 rORMAT(5X,3H N=,I3,6H TLl.Ell.S) END

SUBROUTINE STROB(TI.TL,TH1,A1,B1,C1,UM1.,DTK,01,) TH=THl 0A=A1 C8-B1 C=C1 CUM=UM1 0T0=Tol OU=-UM/(A+B) PRINT 3,

001 Kl = 6,26,4 oDT=<T1-TL)/(K1-1) S2=0 53=0

00 2K2=l,KlCTK-fK2-i;-DTOTyO=U K2oUK2=UK(TK,,,,,,.)- = ( 1*DTK.(k2-1)/K1-1)).uK2 S2 = S2*nP S3=S3-.UK2

2 CONTINUE

TU=S2/5J COLJAABS ( -) oT3 = TL*TU 6 0= 1 . E A/T 3 PRINT 4,1,0,,.,,,0 1 CONTINUE

3 FORMaTISx, 98(-)/5,I ! DUT 1 ! TH T3 !

Z our !/5. 98(-))

A rORMATSX, IHJ , 13, 3H I ,7(E11.5,2H !M RETURNcEND

FUNCTION (,,.,,,,) X=TK-TH ClF(X) 7,7,1

1 y=V-A 4rF(V) 2,4,3

2 X = UM >!/A OGDTD 8

3 \ = - IF(y) 4, 4. 5

4 x=UM

5 =- 61() 7,7,6

6 =./

7 . 6G0TO 9

8 =. (->

9 = ORETURNoENO

h, fi, to, ta, Um, , ( );

t), 1/ -

/ (- U!

( , ) ii. , , Uk ( STROB).

, STROB UK . 249.



1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 [ 81 ] 82 83