Вторая пачка часть 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
Индекс
Элементарные функции
Линейные уравнения
Нелинейные уравнения
Случайные числа
