С. 32 ГОСТ25645.323-88
3000 Z=0.
3IKK - Z) 4. 9.5
4 XI = F
X2 =Z
CO TO 6
5 XI =Z
X 2-F
6F "(XI * \2 )/2 .
CALL PRI (P1.F.N1.N2)
1F(ABS <P1 - P) . LT. 0. 001 » P *
(J.
- P). AND. ABS. (XI - X2>. LT..
•8ВШ * ABS (X I)/ 1(ABS (XI) ♦ ABS <X2))) GOTO!
IF (PI - P) 7.1.8
7 XI - F
CO TO 6
8 X2 * F
CO TO 6
102 PP1 = P
N3-N2
103 PP *(1.+ PPl)/ 2.
CALL QWN (YY. PP)
AI.F-(YY •• 2 - 3.)/6.
HU - FLOAT (N3-1)
WW - YY/ HH • SQRT(HH + ALF)
XP- !./(!. + EXP (2. • WW))
F - ( (1. - 2. • XP) •• 2) • N3 l (4. ♦ XP - 4. * (XP •• 2))
И (N2 EQ. 1 ) 1 ’ l./l
GOTO 150
101 P P I - I - P
N3 - N1
СОТО 103
YY 3 3. 141592 ♦ Р/ 2.
F a (SIN ( YY)/ COS (YY)) •• 2
СОТО 1
20 PRINT 30
30 FORMAT (1ЙХ. ’ОБШИВКА. P < ff. ИЛИ P > 1.ИЛИМ1><- Й’)
STOP
9 F ■ Z
I RETURN
END
SUBROUTINE OWN (X, P)
С ВЫЧИСЛЕНИЕ P КВАНТИЛЯ СТАНДАРТНОГО NUT. I) - НОРМАЛЬНОГО
С РАСПРЕДЕЛЕНИЯ
с*..................................
..
.................................................
II (Р. o r I„OR. LT. (Г.) GOТО 20
X■0.99999F* 12
D -X
IF(P) 4.4.2
2 И ГР - 1.0)7, 12,4
4 X - - #.9999991-*T2
СО TO 12
7 D -P
II (D -0 .5 ) 9.9.8
8 D M .0 -D