Текст подпрограммы и версий 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)