# Вторая пачка часть 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, однако, как ука¬зы¬вается в ра¬боте Грин¬чишина и др

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