Android- : play.google.com
->

0 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

. (= 1 =} . -J - (. 3.16, ).

!1 , " (. . 3.12) -


. 3.16

. - , \-. -. )! (3.88)-(3.91) .- ( --hi hi) (., , 1241).

.. . 24, KOTop(je . (3.88) (3.90) . , , , . , - .

- ;< . , z. , - penieHM . ., 1 iipn ], «» . , 1,

fjbiMH . , " )1 5 v« 3.

- - . , - 14, 241.

3 9. -

--! (. 3.17):

(3.92)

- \

CC(j ,In

ft/.

: ", , Z,

(3.93) (3.94)

(-, , Z, T)t .. ҄.

. . - {, ,, Zf,, );

j I)/i.v, i/n=-(m -1)/, - {k-~\)h.

- , , 2 " V;, t

2. . k ~

I .-

2 iT

! .m.k-

, hi

(. ,

1 3

21 V 3 - ft

cphx

--0,



-

DIMENSION ALF(6),QS(5),TV(ie),KW(iei), «A(2e).(20).C(20).D(2e),G(20),

#u(2a.20.2e).BY(2e).BZ(2e).w(2e}

1. DX.DY,D2 - !

AL.CR - , QV - ALFIS) -

(=. =, Y-e. Y=DY. 2=. 2=D2) QS(S) -

:~ - KX,f,N2 - TAU - -

JV - rV(JV) - KV , £ 2,

Ϩ KVV(KV) -

READ 1, DX,])Y,D2,AL,CR.QV

1 F0RMAT(5Fia.3)

PRINT 2.DX,DY,DZ.AL,CR,QV

2 EORMATC :.3Gie.3/

» Ab=,G10.3. CR=,G10.3, QV=\G10.3} READ 1. (ALF(JJ),N=1,5).(QS(N),N=1.S) PRINT 3,[ALF(N).N=1.8).IQS(N).N=1,5)

3 FORHATC : 750-10.3/ » ;75G10.3)

READ 1. ».TAU. PRINT 4,0,TAU,

4 FORMATS T0=-,FS.3, TAU=,G10.3, TMAX=,1.) READ 5. NX,i,NZ,JV,KV

5 FORHAT(10I5) S.NX,NY.NZ

6 FORMATC :.315) READ 5. (KVV(N),N=1.KV1 PRINT 7,(KVV(N),N=1,KV)

7 FORMATC :7(1015)) READ 1, (TV(N),N=1.JV)

PRINT 8,(TV(N).N=1,) S FORMATC : 7fSG10.3))

2. DO 10 K=1,NZ

DO 10 M=1.KY

DO 10 N=1,KX 10 U(K,M,K)T0

3.

. 3,17

55 56 57 58 59 60 61 62 63 64 65 66 67 58 66 70 71 72 73 74 75 75 77 78 79 80 81 82 83 84

87 86 88 90 81 92 93 94 95 86 97

101 102 108 104 105

KXl-HX-l m=ir:-i NZ1=KZ-1

Hx=ra/ci

HY=DY/NY1 HZ=DZ/KZ1

R1=CR»HX»HX/AL/TAU

R2eCR«HY«HY/AL/TAU

R3=CR#H2#H2/AL/TAU

R4=QV«HX*HX/AL/3

IV*HY»HY/AL,/3

R6=QV*HZ*HZ/AL/3

4. ,..

DO N=1.N«AX

A(N)=1. C(N)=1.

(1)=0,

DO 12 N=2.fiXl

12 BX(N)=-2.-Rl

(1)=-1.-R1/2.-ALF(1)«HX/AL BX (NX) =-1.-Rl/2. - ALF{ 2) t.HX/AL DO 13 H=2,KY1

13 BY(M)=-2,-R2 BY(l)=-l.-R2/2.-ALF(3)«HY/AL BY(KY)=-l.-R2/2.-ALF{4).HY/AL DO 14 K=2,N21

14 BZ(K)=-2.-R3 BZ(l)=~l.-R3/2.-ALF(5)#HZ/AL B2(N2)=-1--R3/2."ALF(8)»H2/AL TIME=0.

J=l Ul

15 TIME=J*TAU

5.

A(NX)=0.

DO 18 K=1.KZ

DO 18 M=I.KY 5,1 D{K>

DO 16 K=2.NX1

16 D(N)=R4tRl.U(N.H,K)

D(1)-QS(1)«HX/Ab+(R4+RI»U(1,H,K))/2.

D(NX)-QS(2).HX/AU(R4+R1»U(NX,H.K))/Z. 5.2

CALL SYSIWdH.A.BX.C.D.G.NX.IER) 5.3 1

DO 17 N-1,NX

17 U{K,H.K)-W(N)

18 COHTINUE A(C)-1.

0 6. Y

tKC. 3.17



A(NY)-e,

DO 21 =1 ,>JZ

DO 21 N=J .m

6.1 : m)

DO 19 M=2,KY1

19 D(H)-R5-R2-U(N,M,K)

D(l )=QS(3).Hl-/AL+CH5tR2.iJ(N.l ,K))/2.

D{N)=QS(4)»KY/AL-f(R5tR2#U(N,KY,K}}/2.

e 2 PFWKHE

CALL SVSTWX.A.Bf.C.D ..IER)

e,3 1

DO 20 H=1,N1

20 U(N.H,K)=W[M)

21 CONTltJUE

A(tJY)=l,

?, Z

ft(N2)=0.

24 M=l ,.41"

DO 24 N=1,NX

7.1 D(K)

22 H=2,N21

22 D(K)=R6+R3.U(N.M.K}

D(1)=QS(5).HZ/AL(RSR3*U(N.M.1)1/2.

D(NZ)=QS(6)«HZ/AL+(R6+R3»U(N,M,NZ))/2.

7,2

CALL SYSTRD(W,A,BZ,C,D.G,NZ.IER)

7.3

23 =1,2

24 CONTINUE

A(NZ)=1.

.

]FiTlME-LT.TV(L}} GO 29

14®

L=L+1

PRINT 25.J.TIME

FORMAT(/ HA :5,

DO 27 1=1,KV

K=KVV(I)

PRINT 26.

26 FORMAT(/ K=.I3)

M 27 N=1.KY

J7 PRIHT 28.*(,a4N..H,.Kt,N=l ,NX)

26 FORMATC H=. I3/(6G11,3.) i

.

IFdlME.GE.THAX) GO 30

J=J+S

GO TO 15

STOP

Piic. i.l

. ttJ, < 1 , / I V,,, ,, I

cphy

n. II - I . k

,, , /- .

Un.m.k- 1 2 4-

cphz

"7

peineHnu nil. . k « , . \-. V,,, ! « , (/», k . -f/eu, « , », ,-.. ,. Vn. ., . ;/,",!,, - < . ;, k ,, U ,\ , . v .

" § 3.5. ) - . """ . . . V,,. i. (-S ?~")- ( 106-121), ,, , (

137). - -- . .1 -"0 (.v, ) - HaV""" , ;1

;! . = " - , *Tcf( k liporoifKaMH no .v in 1, ) ;

st-e . , - .!



0 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



0.0109
.