. " -3 1" ( )

.. "122 1944 ." ( )

. " " ( )

.. " " ( )

.. " " ( )
- ..
.., .., .., .. . .. .: , 1982. 267 c.
( ): primenenevm1982.djvu
<< 1 .. 92 93 94 95 96 97 < 98 > 99 100 101 102 .. 103 >>

/A2**2+A5)*R0*T(I)+H(I) S(I)=(A1+A1+D/A2+A4*D/B+
A5)*R0+S(I) DT=-11604.5/T(I) A=O.
B=O. C=O.
DO 32 J=I1NI
D=G(J)*EXP(DT*E(J))
A=A+D
D=D*E(J)
B=B+D
32 C=D*E(J)+C W1=B/A
CP(I)=CP(I)+(DT*DT*(C/A-
W1**2))*R0
H(I)=H(I)-DT*W1*T(I)*R0 S(I)=(ALOG(A)-DT*W1)*R0+
S(I)
33 F(I)=S(I)-H(I)/T(I) GO TO 500
246
JX= 1 \ NH= 2 \ NK=N
\ NX=2 508
501\ EF=R(26) ****
EC=EF5.0 600
502 CALL POLIN (NH1NK) IF (KIEQ.777) RETURN DO 503 1=1,5 IF (RR(I,8).LE.EF.AND.
* RR(I,9).LE.EC) GO TO 505
503 CONTINUE IF (NK-NH .GE.6) GO TO 504 EF=EF*2.0 EC=EC*2.0 NH=NX GO TO 502
504 NH=NH+1 IF (NH.LE.6) GO TO 504 GO TO 502
505 IF (EF.LT.999.) GO TO 510 IF=1 601 ICP=I
DO 509 1=2,5 IF (RR(1,8).NE.0. .AND.
* RR(I,8).LT.RR(IF,8))IF=I IF (RR(I,9).LT.RR(ICP,9)
* .AND.RR(I,9).NE.0.)ICP=I
509 CONTINUE I=IF
IF (ICP.GT.I) I=ICP
510 NP=NP+1 IF (NP.LE.5) GO TO 506 PRINT 104 GO TO 508
506 CO(NP1I)=T(NH) CO(NP,2)=T(NK)
DO 507 J=1,9 602
507 CO(NP, J+2)=RR(I, J) 603 DO 511 J=I,8
511 CC(NP,J)=AC(I,J) NK=NH
NH=NX C****
IF (NK.EQ.NH) GO TO 508 R( 14)=.77777 GO TO 502
GO TO (700,616,609), JX
JX=2 NH= 1 NK=N MO=O NX=I
IF (K1.GE.200000000)
* GO TO 501 J=O
DO 601 1=18,24,2 IF (R(I) .LE.TO).OR.R(I)
* .GE.T(N)) GO TO 601 J=J+1
G(J)=R(I+1) IF (ABS(G(J)-.77777)
* .LT.1E-4) G(J)=O. E(J)=R(I) CONTINUE
IF (J.EQ.0) GO TO 501 EF=1E+3 EC=1E+3 MO=J
IF (J .EQH) GO TO 603 DO 60 2 M0=2,J DO 602 1=2,J IF (E(I).GT.E(I-I))
* GO TO 602 A=E(I-I) E(I-I)=E(I) E(I)=A A=G(I-D G(I-I)=G(I) G(I)=A CONTINUE NH=1 NX=I M1 = 1 JX= 3
247
604 DO 605 I=MH1N
IF (T(I)-E(M1).GT.0.1)
* GO TO 606
605 CONTINUE NK=N NKK=NK JX= 2
GO TO 607
606 NK=I-I
IF (NK-NH.I/T.1) GO TO 608 IF (ABS(T(NK-I)-E(M1))
* .LT.0.D NK=NK-I NKK=NK
607 IF (NK-NHGEI ) GO TO 502
608 PRINT 105 GO TO 999
0***
609 A=E(MI )*1E-4 B=CO(NP,3)+CO(NP,4)+
* CO(NP,4)AIiOG(A)-
* CO(NP,5)*A**(-2)+((4.0*
* CO(NP,9)*A+3.0*CO(NP,8))*
* A+2.0*CO(NP,7))*A R(13)=B+G(M1)/E(M1)
B=((((3.0*CO(NP,9)A+
* 2.*CO(NP,8))*A+CO(NP,7))*
* A+CO(NP,4))*A-CO(NP,6)
* -2.0*CO(NP,5)/A)*1E+4 R(15)=B+G(M1) R(M)=E(MI)
M1=M1+1 NH=NKK+1 NX=NH
IF (Ml .LE.MO) GO TO 604
NK=N
JX= 2
GO TO 607
616 IF (MO.LT.D GO TO 610 DO 617 J=2,5 IF (ABS(CO(J,1)-.77777)
* .LT.1E-4) GO TO 610 IF (CO(J,D-CO(J-1,2)
* . .0) GO TO 617 DO 615 1=1,5
IF (E(I).GT.CO(J,1)
* .OR.E(I).LT.CO(J-I ,2))
* GO TO 615 CO(J1I)=E(I) CO(J-1,2)=E(I)
615 CONTINUE 617 CONTINUE
610 DO 614 I=I9N DO 611 J= 1 ,NP
IF ((T(I)-CO(J1D) .GT.
* -.1.AND.(CO(J,2)-T(I)) .GT.-.1) GO TO 612
611 CONTINUE T(I)=.77777 GO TO 614
612 IF (J.EQ.NP.OR.I.EQ.1)
* GO TO 613
IF (ABS(T(I)-T(I-D)
* .LT.0.1) J=J+1
613 A=T(I)*1E-4 S(I)=CO(J,3)+00(J,4)+
* CO(J,4)*ALOG(A)-CO(J,5)
* *A**(-2)+((4.0*CO(J,9)
* *A+3.0*CO(J,8))*A+
* 2.0*CO(J,7))*A H(I)=((((3.0*CO(J,9)*A+
* 2.0*CO(J,8))*A+CO(J,7)) A+CO(J,4))A-CO(J,6)-
* 2.0*CO(J,5)/A)*1E+4
CP(I)=CO(J,4)+2.0*(((2.*
* CO(J,9)*A+CO(J,8))*
* 3.*A+CO(J,7))*A+CO(J,5)
* /(A*A)) F(I)=S(I)-H(I)/T(I)
614 CONTINUE C***
700 IF (ABS(T(D-298.)
* .GT.0.5) GO TO 999 IF (ABS(RU)-.77777)
* .LE. 1 E-4) R(4)=H(1) 999 RETURN
248
101
102
109 FORMAT (' ..
* * .1X1
* ', \* ',
', 4 '') FORMAT (//3X1 ' ' ,
* X
* ' T '.
* (' ,
* T) CP(T) H(T)'
* S(T)V/) FORMAT (//3X1 B * ,
* X
* () [CP(T),
* '(T)1S(T)] '
* T'/)
103 FORMAT ( EHG'
* ' ,
* '
* ')
104 FORMAT (,'',
* ECKHE ',
* ' 1,
* ' '/)
105 FORMAT (,'
* ' *
* * ,
* ' '/)
106 FORMAT (//.'1a
* ' * ,
* ' ,
* " .. '//) END
SUBROUTINE SYSTEM (N) DOUBLE PRECISION T1A1S1
* W,Z
COMMON /HEZ/T(10),RZ(16)
* //Y(IIO)1A(IO1II) DIMENSION M(IO)
Q*****************************

Q******************************
N1=N+1
DO 1 I=I1N
1 M(I)=I NN=N-1
Q***
DO 8 K=I9NN
KK=K-H
W=O.
DO 3 J=K9N DO 3 I=K9N Z=DABS(A(I9J)) IF (W-Z) 2,3,3
2 W=Z 1=1 K2=J
3 CONTINUE
IF (K.EQ.K1) GO 5 DO 4 J=K9NI S=A(K9J) A(K9J)=A(KI9J)
4 A(KI9J)=S
5 IF (K.EQ.K2) GO TO 7 DO 6 I=I9N S=A(I9K) A(I9K)=A(I9K2)
6 A(I9K2)=S K3=M(K) M(K)=M(K 2) M(K2)=K3
7 DO 8 I=KK9N Z=A(19K)/A(K9K) DO 8 J=K9NI
8 A(I ,J)=A(I9J) -Z*A(K9J)
?***
DO 11 J= 1,N
NN=J-I
I=NI-J
.. 1554
249
S=O.
IF (N.EQ.I) GO TO 10 DO 9 K=I9MN KI=NI-K 9 S=A(KI9N1)*A(I9KI)+S
10 K3=M(I)
A(I9N1) = (A(I9N1)-S)/A(I9 * I)
11 T(K3)=A( I9NI) RETTJBN
END
SUBROUTINE FLAGS (CS9JA9
* JB)
COMMON /SYMB/ L(94),A(44)
* // AZ(60) DIMENSION CS(15)9N1(5),
* N2(5)9NS(15)
<< 1 .. 92 93 94 95 96 97 < 98 > 99 100 101 102 .. 103 >>
 
 
2010 PiroChem. .