|
Текст подпрограммы и версий sf65r_p.zip |
Тексты тестовых примеров tsf65r_p.zip |
Вычисление логарифмической производной ψ (1 + ix) гамма - функции.
SF65R вычисляет логарифмическую производную гамма - функции (пси - функцию) ψ (1 + ix).
Подпрограмма - функция SF65R представлена в Библиотеку Вычислительным центром AH CCCP (группа Л.Н.Кармазиной) в виде процедуры на АЛГОЛе - 60. Перевод этой процедуры на стандарт языка ФОРТРАН осуществлен в НИВЦ МГУ. Алгоритм взят из кандидатской диссертации:
Г.Неметх, "Разложение обобщенных гипергеометрических функций по полиномам Чебышева", ОИЯИ, Дубна, 1971.
function SF65R(X :Real): Complex;
Параметр
| X - | заданное значение аргумента x (тип: вещественный). |
Версии: нет
Вызываемые подпрограммы: нет
Замечания по использованию: нет
Unit tsf65r_p;
interface
uses
SysUtils, Math, { Delphi }
Lstruct, Lfunc, UtRes_p, SF65R_p;
function tsf65r: String;
implementation
function tsf65r: String;
var
X :Real;
A :Complex;
label
_1;
begin
Result := ''; { результат функции }
X := 2.0;
A := SF65R(X);
Result := Result + Format('%s',[' X=']);
Result := Result + Format('%20.16f ',[X]);
Result := Result + Format('%s',[' A=']);
Result := Result + Format('%20.16f %20.16f ',[A.re,A.im]) + #$0D#$0A;
_1:
UtRes('tsf65r',Result); { вывод результатов в файл tsf65r.res }
exit;
end;
end.
Результат: Z = (0.7145915153687, 1.320807282642)