о п п п
ГОСТ 25645.323-88 С 31
CALL QWH (XR, AL1, N1)
РР8 = N1/ XR
CALL 0WS(T.AL.N1)
CALL OWF (F, AL, N2. N1)
A5 “ Y - T • SORT (SY/(N1 ♦ 1))
A6 » Z • SORT (SE)
A? = (FP8 • FF - 1. ♦ 1./ (F ♦ FF) ♦ (1. - FP8/ F)) • • 0 .5
y
p
g
=■
аз
A6 • A7
117 FORMAT (SX. ’LG (A) » \ F7. 3.5X. ’B - *. F7. 3)
PRINT 4 ’
PRINT 117.AK.BK
A8 = 10 •• AK
13fl FORMAT (7X. *A•=*. E10. 3)
PRINT 130, A8
118 FORMAT (SX. ’НИЖНИЙ ДОВЕРИТ. ПРЕДЕЛ - ’ . F7. 3)
PRINT 118. YP
119 FORMAT (5X.‘НИЖНИЙ ТОЛЕРАНТНЫЙ ПРЕДЕЛ » ’ , F7. 3)
PRINT I19.YPG
B5 - 10 •• YP
B6 = 10 •• YPG
PRINT 200. B5. B6
1010 CONTINUE
1000 CONTINUE
GOTO 160
130 PRINT 3
PRINT 2.G. RR.QU.RN
160 STOP
END
SUBROUTINE OWF (F. P. N1. N2)
ВЫЧИСЛЕНИЕ P- КВАНТИЛЯ F РАСПРЕДЕЛЕНИЯ C NI.N2 СТЕПЕНЯМИ
СВОБОДЫ
1F(N2. LF-0. OR. N1. LE.0. OR. P. LT.0. .OR. P. GT. 1.) GO TO 20
IF(N2. EQ. LAND. Nl.EQ. 1) GOTO 100
II (N2. EO- 1) GOTO 101
IF(N1. FO. I) GOTO 102
CALL QWN (X. P)
A * (X • X - 3.) / 6.
Ml -N1 - I
М2 - N2 - 1
H - 2./ (1. / FLOAT (Ml) ♦ 1. / FLOAT (М2) )
W - X / 11 • SORT (H ♦ A) - (1. I FLOAT (Ml) - 1. / FLOAT (М2) )M A +
• + 5. / 6. - 2. / (3. eH))
fF " EXP (2. • W)
- 130 CALL PRГ (PI, F, N1, N2)
’II (ABS (PI - P). EQ. 0. 0) GO TO 1
11=0.01
If (PI. GT p) II - - 1.* II
2 Z « P * H
1F(Z. Lli.0.)GOTO 3000
H -H «2.
CALL PRF (P1.Z.NI.N2)
II {(PI - P) • H) 2.9, 3