Вторая пачка часть 142
2346331+i*0.1580554
7/4 1.517608e-2-i*3.706224 1.517607e-2-i*3.706224
-7/4 -1.517608e-2-i*3.706224 -1.517607e-2+i*3.706224
R:=1;
T:=1-X;
Q1:=Q;
I:=Q;.
REPEAT
DEC (I);
IF I>0 THEN
BEGIN
Q1:=I;
R:=R*(Q1+1)/T/(P+Q1);
S2:=S2+R;
END
UNTIL I = 0;
R:=1;
S1:=1;
I:=0;
REPEAT
I:=I+1;
IF R>=EPS*S1 THEN
BEGIN
R:=R*X*(I-Q1)*(P+I-1)/I/(P+I);
S1:=S1+R;
END;
UNTIL R >= EPS*S1;
U:=Q1;
T:=G(U);
T1:=T;
U:=Q1+P;
R:=G(U);
R1:=R;
I:=Q1;
WHILE I<= Q-0.5 DO
BEGIN
T1:=T1*I;
R1:=R1*(I+P);
INC (I);
END;
T:=X P*(S1*R/P/T+S2*R1*(1-X) Q/Q/T1);
U:=P;
T:=T/G(U);
IF W=1 THEN `
BEGIN
UBETF:=1-T;
EXIT
END;
UBETF:=T
END { *** UBETF *** }.
Процедура-функция UBETF получена путем перера¬бот¬¬¬ки и перевода вначале на язык FORTRAN [Белашов, 1997], а затем на язык PASCAL программы вы¬чис¬ле¬ния , опубликованной в работе Гринчишинa и др. (1988), и протестирована на IBM PC/AT-286 для тех же, что и в указанной последней работе, значений x, p и q. Р嬬зультаты вычислений функций I0.7(0.5; 0.5) = 0.63098988, I0.2(2;1.5) = =0.06979572 при заданной точности ? = 1e-8 с точ¬ностью до вось¬ми десятичных знаков совпадают с по¬лученными в ра¬ботах Белашова (1997), Гринчиши¬на и др. (1988).
§ 2. НЕКОТОРЫЕ ИНТЕГРАЛЬНЫЕ ФУНКЦИИ
Интегральная показательная функция Еk(z) опре¬де¬ля¬ется интегралом [Справочник ..., 1979]
0, 1, 2, ...; z = x + iy, x > 0 .
Наряду с функцией Еk(z) часто используется свя¬зан¬ная с ней ком¬плексная интегральная функция
Wk(z) = U + iV = zk ez Ek (z), (5.11)
представляемая в виде непрерывной дроби [Спра¬воч¬ник ..., 1979]
вычислить которую можно итерационным методом в соответствии с формулами [Библиотека ..., 1975]
;
; ;
; ;
Итерационный процесс завершается, когда вы¬пол¬ня¬¬ет¬ся , где ? - заданная точность вы¬чис¬¬ле¬ний.
Данный алгоритм в соответствии с областью оп¬ре¬де¬ле¬ния z позволяет получить результат на ве-щественной по¬лу¬плос¬кости Re z > 0, однако, как ука¬зы¬вается в ра¬боте Грин¬чишина и др
Индекс
Элементарные функции
Линейные уравнения
Нелинейные уравнения
Случайные числа
