ГОСТ 25645.115—84 С. 25
GOTO 4
3 КР= АР
4 DKP = KP*KP
A=Q(22, J)+Q(2’3, J)*KP+Q(24, J)*DKP
K4= Q(18, J)+Q(1I9, J)*H+Q(20, J)*HD+Q(21f J)*HD*H
K4= l.D0+K4*A
RO= K4*RO
GOTO 8
5 CONTINUE
DO 6 1=1, 9
IF (H.GE.U(10—I, 1)) GOTO 7
6 CONTINUE
PRINT 120
20 FORMAT (5X, ’H <0’)
RETURN
7 1=10—1
DH= H—U(I, 1)
DHD= DH*DH
RO=U(I, 2)*DEXP(■—U(I, 3)*DH+U(I, 4)*DHD)
8 RU= RO/GR
RETURN
END
SUBROUTINE VKM(FSH)
IMPLICIT REAL*8(A—H, O—Z)
DIMENSION XF(7)
DIMENSION Rl,(27,7), R2(27,7), R3>(27,7)
DIMENSION Rll (54), R12(81), R13(54)
DIMENSION R21l(54i), R2’2(81), R23i(64)
DIMENSION R31‘(54), R32(81), R33(54)
EQUIVALENCE (R1 (1,1), Rll/(1))
EQUIVALENCE (RI(1,3), Rlf2(l))
EQUIVALENCE (Rl.(l,6), R13(l))
EQUIVALENCE (R(2l(Ll), R2l(l))
EQUIVALENCE (R2(l,3), R22(l))
EQUIVALENCE (R2(i,6), R23(l))
EQUIVALENCE (R3i(l,l), R31/(l))
EQUIVALENCE (R3(l,3), R32f(l))
EQUIVALENCE (R3(l,6), R33(1))
COMMON/UKOEF/U(8,4), GR, WZ
COMMON/QKOEF/Q(27,8)/YEAR/P (38)
DATA XF/87.5D0, l[I2.5lDO, 187.5D0, 162.5DQ, 187.5D0, 225.D0, 40O.DO/
DATA R11/75.DO, —18.2991(D0, .70O9DO, 115.3430D0,
* —.6828D0, .5576D—2, .95238D—6,
* —4.384D0, .080629DO, —.4925D—3, .1I042D—5,
* L5D0, J0O6DO, .541O52D0,
* —5.ШШХ), J06258D0, —.1672D—3,