. " -3 1" ( )

.. "122 1944 ." ( )

. " " ( )

.. " " ( )

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

UM *********************

IS=O
DO 3 1=1,15 DO 1 J=I,44 IF (CS(I).EQ.A(J))
* GO 2
1 CONTINUE GO 99
2 IF (IS.EQ.1) GO 3
IF (J.NE.33.AND.J.NE.31)
* GO 3 IS=I NK=I-I
3 NS(I)=J DO 4 J=I,5 NKJ)=O
4 N2(J)=0 ****
IF (NS(I).NE.29.OR.NS(2)
* .NE.30) GO 5
N1(1)=4 J
N2(1)=1 /
GO TO 19 J
5 N1(1)=2
IF (NS(2).NE.28.AND.NS(3)
* .NE.28) GO TO 6 IS= 3 I NI(D=I f
IF (NS(2).EQ.28) GO TO 6 IS= 4
IF (NS(I) NE.26 .AND.NS(D
* .NE.27) GO TO 6 N2(1)=NS(2)-35 N1(1)=3
IF (NS( D.EQ.27) N1( 1)=4
Q****
6 M=2
DO 18 J=IS9NK
7 IF (NS(J)-26) 898915
8 IF (N2(M)) 9,9,11
9 IF (N1(M)-26) 12,12,10
10 N2(M) = 1
11 M=M+1 GO TO 7
12 IF (NI(M)) 13,13,14
13 NI(M)=NS(J) GO TO 18
14 N1 (M)=NI (M)*26+NS( J) GO TO 18
15 IF (N2(M)) 16,16,17
16 N2(M)=NS(J)-35 GO TO 18
17 N2(M)=N2(M)*10+NS(J)-35 M=M+1
18 CONTINUE
IF (N2(M).EQ.O) N2(M)=1
Q****
19 DO 22 J=2,5 IS=O
IF (NKJ).EQ.O) GO TO 21
DO 20 N=1,94
IF (NI(J) .EQ.L(N)) IS=N
20 CONTINUE
250
IF (IS.BQ.O) GO TO 99 M N1(1)=N1(1)*100+IS *2 N2(1)=N2(1)*100+N2(J) JA=NKD JB=N2(D
\
RETURN
9? JA=7T7 RETURN END
SUBROUTINE POLIN (N1,N2) DOUBLE PRECISION A1S9TR1
* W,X1,X29X3,X4,X5,B COMMON //AA(60),ID(5),
* R( 599) 9A(1091D9C( 598)
* /DA/G(1100),F(IOO)9
* T(IOO)9CP(IOO)9H(IOO)1
* SS(100)9G1(30)9M2(73)
* /REZ/TR(10)9S(8) DIMENSION W(100),X1(100),
* X2(100)9X3(100)9X4(100)9
* X5(100)9B(110) EQUIVALENCE (A(I)9B(D)
Q*****************#***

C
DO 1 I=NI9N2 W(I)=T(I)*1D-4
XKI) = I.0/W(I) X2(I)=DLOG(W(I)) X3(I)=X1(I)**2 X4(I)=W(I)**2
1 X5(I)=X4(I)*W(I) DO 18 J=I,5
IF ((N2-N1).GT.J) GO TO 2
R(J,8)=1E+10
E2=1E+10
IF (J-I) 18,20,18
2 J1=J+2 J2=J+3
DO 3 11=1,7 TR(H)=O. DO 3 12=1,8
3 A(H ,12)=0. DO 11 I=N 1,N2
IF (ABS(F(I)-.77777)
* .LE.1E-4) GO TO 5 S(1)=1.0 S(2)=X1(I) S(3)=X2(I) SU)=W(I) S(5)=X3(I) S(6)=X4(I) S(7)=X5(I) S(J2)=F(I)
DO 4 12=1,J2 I3=(I2-1)*10 DO 4 11=1, 14=13+11
4 B(I4)=S(I1)*S(I2)+B(I4)
5 IF (ABS(SS(I)-.77777)
* .LE.1E-4) GO TO 7 S(I) = LO
S( 2)=0.
S(3)=1.0+X2(I)
S(4)=W(I)+W(I)
S(5)=-X3(I)
S(6)=3.0*X4(I)
S(7)=4.0X5(I)
S(J2)=SS(I)
DO 6 12=1,J2
I3=(I2-1)*10
DO 6 11=1,Jl
14=13+11
6 B(I4)=S(I1)*S(I2)+B(I4)
7 IF (ABS(CP(I)-.77777)
* .LE.1E-4) GO TO 9 S(I)=O.
S(2)=0.
S(3)=L0
S(4)=W(I)+W(I)
S(5)=X3(I)+X3(I)
S(6)=6.0X4(I)
251
S(7)=12.0*X5(I) S(J2)CP(I) LO 8 12*1,J2 ǻ(2-1)* DO 8 11I9JI 14=13+11
8 B( 14)=S(11)*S(12)+B(14)
9 IF (ABS(H(I)-.77777)
.LE.1E-4) GO TO 11 S(I)=O.
S(2)=-1.0
SO)=W(I)
S(4)=X4(I)
S(5)=-2.0Xl(I)
S(6)=2.0*X5(I)
S(7)=3.0*X4(I)*X4(I)
S(J2)=H(I)*1D-4
DO 10 12=19J2
13=(12-1)^10
DO 10 H = I9JI
14=13+11
10 B(I4)=S(Il)*S(I2)+B(I4)
11 CONTINUE IF(ABS(GK 13)-.77777)
.LE.1E-4.0R.G1O4)-
.77777.LE.1E-4) GO TO 12 X=G1(14)1E-4 A(I9D=LO
A(1,2)=0.
A(1,3)=ALOG(X)+t.O
A(1,4)=X+X
A(1,5)=-1.0/(1?)
AO,6)=3.0??
A(1,7)=4.0XXX
A(I ,J2) =01(13)
12 IF (ABS(GK 15)-.77777)
.LE.1E-4.0R.GK14)-
77777LEnEH) OO TO 13 X=G1(14)1E-4
A(2,1)=0. A( 2,2)=-1.(* A(2,3)=X A(2,4)=X*X
A( 2,5)=-2.OA I
A(2,6)=2. 0*X3
A(2,7)=3.0*X*4 j
A(2,J2)=G1(15)*1E-4 j
13 IF (A(I9I)EQ.0..OR./
* A(292).EQ.0.)GO TO 20 GALL STSTEM (J1) R(J9D=TR(I) R(J,2)=TR(3) R(J93)=TR(5) R(Ji4)sTR(2) R(J,5)=TR(4) R(J,6)=TR(6) R(J,7)=TR(7)
EI=O. E2=0. E3=0. E4=0.
DO 17 I=N1,N2
IF (T(I)-.77777 .LE.1E-4)
* GO TO 17 X=T(I)*1E-4
IF (ABS(F(I)-.77777)
* .LE.1E-4) GO TO 14 1=((R(J,7)*X+R(J96))X+
* R(J,5))*X+R(J,1)+ALOG(X)
* R(J,2)+(R(J,3)/X+ R(J,4))/X
T=ABS(F(I)-T)
IF (EI .GT .T) GO TO 14
Et=T
C(J,2)=T(I)
14 IF (ABS(SS(I)-.77777)
* .LS.12-4) QO TO 15 T=( (4*CNR(J97)X+3.a
* R4Jr6))**K?.(r*RtJ95))X+ R(J,1)+R(Jf2>**(1.0+
* ALOG(X))-R(J93)/(XX) T=ABS(SS(I)-T)
IF (E3.GT.T) GO TO 15 E3=T
C(J,4)=T(I)
15 IF (ABS(H( 1)-.777779)
252
[ * .LE.1E-4) GO 16 \ Y=((((3.0*R(J,7)*X+ \ * 2.0*R(J96))*XR(J,5))*X* \ * R(J,2))*X-R(J,4)-l * 2.0*R(J,3)/X)*1E+4 Y=ABS(H(I)-Y) IF (?4GT9Y) GO TO 16 E4=Y
C(J96)=T(I)
16 IF (ABS(CP(I)-.77777)
* .LE.1E-4) GO TO 17 Y=R(J,2)+2.*(((2.*R(J97)*
* X*R(J96))*3.0*X+R(J,5))*
* XR(J,3)/(X*X)) Y=ABS(CP(I)-Y)
IF (E2.GT.Y) GO TO 17 E2=Y
C(J,8)=T(I)
17 CONTINUE C(J1I)=EI
IF (E3.GT.E1) E1=E3
Y=E4*2.0/( T(N1) T (N2 ) )
IF (Y.GT.E1) EI=Y
C(Jf3)=E3
C(J,5)=E4
C(J,7)=E2
R(J98)=E1
18 R(J,9)=E2
19 RETURN
20 PRINT 100 M2(16)=777 RETURN
100 FO (' IHW1
* ' Ȳ 1,
* ' XBATiET ѫ,
* ' AAHHHXD END
SUBRDUTINE NIlRS(AA9B) DIMENSION (8) ,L(8),C(14) DATA /10,11,12,1,
* 14915,16,17,189199
* 1H-f1H.91H ,IH*/
****************** ************
40PMATHOE
**** ** ** ****** ****************
A=ABS(AA)
IF (ABS(A-.77777) LT.
* 1-4) GO TO 19
Ii (AA.LE.-1E+7.CR.A.GE.
* 1+8) GO 16
IF (A.LT.1E-6) GO 14
IA=A
A=A-IA
IB=A* 1+7
1=10000000
II (IB.LT.IC) GO 1
IA=IA-i-I
IB=IB-IC
1 N=1
IF (IA.EQ.O) GO 3 DO 2 1=1,8 =(1*4*1)/(1*4) IF (K.EQ.0*AND.N.EQ.1)
<< 1 .. 93 94 95 96 97 98 < 99 > 100 101 102 .. 103 >>
 
 
2010 PiroChem. .