Текст подпрограммы и версий ( Фортран ) sf19r.zip sf19d.zip |
Тексты тестовых примеров ( Фортран ) tsf19r.zip tsf19d.zip |
Текст подпрограммы и версий ( Си ) sf19r_c.zip sf19d_c.zip |
Тексты тестовых примеров ( Си ) tsf19r_c.zip tsf19d_c.zip |
Текст подпрограммы и версий ( Паскаль ) sf19r_p.zip , sf19e_p.zip |
Тексты тестовых примеров ( Паскаль ) tsf19r_p.zip , tsf19e_p.zip |
Вычисление неполных эллиптических интегралов первого рода.
Функция SF19R производит вычисление интегралов вида:
T F (K, T) = ∫ (1 - K2 sin2 t) -1/2 dt для 0 ≤ K2 ≤ 1 и 0 ≤ T ≤ π/2 0
J.F.Hart, E.W.Cheney, C.L.Lawson et al, Computer Approximations, Wiley, New York, 1968.
FUNCTION SF19R (RK2, T, IERR)
Параметры
RK2, T - | заданные значения аргументов K2 и T (тип: вещественный); |
IERR - | целая переменная, служащая для сообщения об ошибках, обнаруженных в ходе работы функции; при этом: |
IERR=65 - | когда значение RK2 меньше 0. или больше 1.; значение функции полагается равным 3.4E38; |
IERR=66 - | когда значение T меньше 0. или больше π / 2; значение функции полагается равным 3.4E38; |
IERR=67 - | когда RK2 = 1 и T = π / 2 (значение подинтегральной функции pавно 0); значение функции полагается равным 3.4E38. |
Версии
SF19D - | вычисление неполных эллиптических интегралов первого рода с повышенной точностью. |
Вызываемые подпрограммы
UTSF10 - | подпрограмма выдачи диагностических сообщений при работе функции SF19R. |
UTSF11 - | подпрограмма выдачи диагностических сообщений при работе функции SF19D. |
Замечания по использованию
Для функции SF19D параметры RK2 и T имеют тип DOUBLE PRECISION. Значение SF19D при IERR ≠ 0 полагается pавным 1.7D308. |
RK2 = 0.5 T = 3.14159265358 / 4 Y = SF19R (RK2, T, IERR) Результаты: Y = 0.826017876245, IERR = 0