# Вторая пачка часть 152

Вычисление отношения неполной бета-функции (5.9) при x = z ? Re z, p ? Re p, q ? Re q ре¬а¬ли¬зо¬ва¬но в про¬цедуре-функции UBETF на основании раз¬ло¬же¬ния в степенной ряд, при этом в целях оп¬ти¬ми¬за¬ции про¬цесса сходимости ряда при x > 0.5 ис¬поль¬зу-ется со¬¬от¬¬ношение симметрии (5.10). В процедуре пред¬ва¬ри¬тель¬но про¬веряется соответствие вводимых значений ам¬плитуды и аргументов областям их оп¬ре¬де¬ле¬ния [см. вы-ражение (5.7)]; в случае если неравенства (5.7) не удов¬ле¬тво¬ря¬ют¬ся, осуществляется выход из процедуры к внешней про¬грам¬ме обработки ошибки. Точность вы¬чис¬ления задается при вызове процедуры-функции. Формальные параметры процедуры. Входные: x, p, q (тип real) - амплитуда и аргументы функции ; eps (тип real) - задаваемая точность. Выходной: ubetf (тип do¬ub-le) - вычисленное с точностью eps значение от-но¬шения не¬полной бета-функции . FUNCTION UBETF (X,P,Q,EPS : REAL ) : REAL; VAR W,T,T1,S1,S2,R,R1,Q1,I,U : REAL ; LABEL A2, A3, A4, A5, A6; {*** ВНАЧАЛЕ ОПРЕДЕЛЯЕТСЯ ПРОЦЕДУРА-ФУНКЦИЯ G, ИСПОЛЬЗУЕМАЯ В ОСНОВНОЙ ПРОЦЕДУРЕ-ФУНКЦИИ UBETF, НАЧИНАЮЩЕЙСЯ С МЕТКИ BEG; ***} FUNCTION G(U: REAL) : REAL; VAR GL : REAL; BEGIN G:=1; WHILE U>1 DO BEGIN U:=U-1; G:=G/U; END; G1:=(((((((((-0.00000018122*U+0.000001328554)*U -0.000002625721)*U-0.000017527917)*U+ 0.000145624324)*U-0.000360851496)*U- 0.000804341335)*U+0.008023278113)*U -0.017645242118)*U-0.024552490887)*U; G:=G*((((G1+0.191091101162)*U- 0.233093736365)*U-0.422784335092)*U+ 1)*U*(1+U); G:=1/G END { *** G *** }; { **** НАЧАЛО ОСНОВНОЙ ПРОГРАММЫ *** } BEGIN IF X>1 OR X<0 OR P<0 OR Q<0 THEN BEGIN WRITE (“ОШИБКА В ИСХОДНЫХ ДАННЫХ ‘); EXIT; END; IF X=0 OR X=1 THEN BEGIN UBETF:=X; EXIT; END; W:=0; IF X>0.5 THEN BEGIN T:=P; P:=Q; Q:=T; X:=X-1; W:=1 END; S2:=0; X0 := X; YR[1] := Y0; REPEAT F1 := H * FUNC (X0,Y0); F := H * FUNC (X0+H/2,Y0+F1/2); F1 := F1 + F*2; F := H * FUNC (X0+H/2,Y0+F/2); F1 := F1 + F*2; F := H * FUNC (X0+H,Y0+F); F1 := F1 + F; YR[I]:=YR[I-1] + F1 / 6; Y0 := YR [I]; X0 := X0 + H; INC (I); UNTIL I>N; END

Индекс
Элементарные функции Линейные уравнения Нелинейные уравнения Случайные числа