О о и о
ГОСТ 25645.323-88 С. 35
СО Г0 11
\Ш М- (N2 - 3) / 2
DO 1? I * 1.М
17 С =С • FLOAT (2 4 I - I) / 2.
16 В = В /(С * 177245)
11 0 = 1. - А + В
18 Р = 1. —Q
СО ТОМ
4fW0 X * ( ( F * * d /3.) I*
(1
-
2./
(9. •
FLOAT (N2)
))
- (1. - 2.
/ (9. •
‘ FLOAT INI))»/ SORT \l. Jl9. * N1) + № *♦ (2.1 3.) ) • 2. / (9. ♦ N2) )
CAl.l PRN (P, X)
СОТОЗЛ
14 PRINT 20
20 FORMAT <I0X.-ОШИБКА. N X -=feИЛИ N2 < Ч>ИЛИГ < 0* )
STOP
30 CONTINUE
RETURN
END
SUBROUTINE OWHfX. P. N)
вычисление
p
квантиля
*
c iii
2*
р
асп
р
еделения c n степеням и
СВОБОДЫ
ir<N. vy.fi. OR. P. V f.fi.. OR.
P.
GT.
1.)
GO TO
10
CALL OWN (Z, P)
X ’N * <(1. ~ 2.1 <9. * I’LOAT (N)) + Z • SORT <2.1<9. ■FLOAT (N) ) ))
* * *3)
IHX. 1.Е.Д.) X -Я.5 * a * SORT (FLOAT (2 * N) - I)) ** 2
CALL PRH (PI, X. N)
IE(ABS (PI ~P>. ЕО.0.ЛСОТО1
H -11. 01
(I (PI CT. P)H = - !. • H
2 Z - X + H
11(2. LF.0.)GOTO 3006
H* H * 2.
CALL PRH (PI. 2, N)
IF((P1 - P) • H) 2,9. 3
3006 Z -fi.
3 II*(X - 2)4.9,5
4 XI * X
X2- Z
CO TO 6
5 XI -Z
X2 •*X
6 X“ (XI * X2J/2.
CALL PRH (PI. X, N)
IKABS (PI - P). LI. 0. fifil * P * (1. - P). \ND. ABS (XI X2). LT.. 0U01 *
* ABS(X1)/(ABS(X1) ♦ ABS (X2))> GO IО I
1F(P1 PI 7.1.8
7 XI -X
(ЮТО6
8 X2 - X
GO TO 6
10 PRINT 2Д
2fi FORMAT(10X, ОШИБКА. N < - jCTИЛИ
P
BHF. (fi, 1)*)
STOP
9 X »Z