26
X3=YY
IF(RON.LT.0.66D0.OR.RON.GT.1D0)Z=0D0
IF(X2.LT.0D0.OR.X2.GT.0.2D0)Z=0D0
IF(X3.LT.0D0.OR.X3.GT.0.15D0) Z=0D0
IF(Z.EQ.0D0) GO TO 133
X1=1D0-X2-X3
Х11=Х1*Х1
X12=X1*X2
X13=X1*X3
X22=X2*X2
X23=X2*X3
X33=X3*X3
Z=1D0-(.0741*RON-.006D0-.063*YA-.0575*YY)**2
BMNG=24.05525*Z*RON
Y1=1D0-YA-YY
BMY=(BMNG-28.0135*YA-44.01*YY)/Y1
С Расчет теплоты сгорания эквивалентного углеводорода (Н)
H=47.479*BMY+128.64D0
RETURN
3 Т=Т1
ТС=Т1-Т0
P=PRESS
IF(PRESS.LE.0D0.OR.PRESS.GT.12D0)Z=0D0
IF(T1.LT.250D0.OR.T1.GT.340D0)Z=0D0
IF(Z.EQ.0D0) GO TO 133
CALL B11BER(T,H,B11)
CALL BBER(T,B11,B,Z)
IF(Z.EQ.0D0) GO TO 133
CALL CBER(T,H,C,Z)
IF(Z.EQ.0D0) GO TO 133
CALL ITER2(P,T,B,C,Z)
133 RETURN
END
SUBROUTINE B11BER(T,H,B11)
IMPLICIT REAL*8(A-H,O-Z)
COMMON/BBLOK/BR11H0(3),BR11H1(3),BR11H2(3),BR22(3),BR23(3),BR33(3)
T2=T*T
B11=BR11H0(1)+BR11H0(2)*T+BR11H0(3)*T2+
*(BR11H1(1)+BR11H1(2)*T+BR11H1(3)*T2)*H+
*(BR11H2(1)+BR11H2(2)*T+BR11H2(3)*T2)*H*H
END
SUBROUTINE BBER(T,B11,BEFF,Z)
IMPLICIT REAL*8(A-H,O-Z)
COMMON/BBLOK/BR11H0(3),BR11H1(3),BR11H2(3),BR22(3),BR23(3),BR33(3)
COMMON/ZETA/Z12,Z13,Y12,Y13,Y123
COMMON/XBLOK/X1,X2,X3,X11,X12,X13,X22,X23,X33
T2=T*T
B22=BR22(1)+BR22(2)*T+BR22(3)*T2
B23=BR23(1)+BR23(2)*T+BR23(3)*T2
B33=BR33(1)+BR33(2)*T+BR33(3)*T2
BA13=B11*B33
IF(BA13.LT.0D0) THEN