|
Текст подпрограммы и версий ( Фортран ) sf20r.zip sf20d.zip |
Тексты тестовых примеров ( Фортран ) tsf20r.zip tsf20d.zip |
|
Текст подпрограммы и версий ( Си ) sf20r_c.zip sf20d_c.zip |
Тексты тестовых примеров ( Си ) tsf20r_c.zip tsf20d_c.zip |
|
Текст подпрограммы и версий ( Паскаль ) sf20r_p.zip , sf20e_p.zip |
Тексты тестовых примеров ( Паскаль ) tsf20r_p.zip , tsf20e_p.zip |
Вычисление неполных эллиптических интегралов второго рода.
Функция SF20R производит вычисление интегралов вида:
T
E (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 SF20R (RK2, T, IERR)
Параметры
| RK2, T - | заданные значения аргументов K2 и T (тип: вещественный); |
| IERR - | целая переменная, служащая для сообщения об ошибках, обнаруженных в ходе работы функции; при этом: |
| IERR=65 - | когда значение RK2 меньше 0. или больше 1.; значение функции полагается равным 3.4E38; |
| IERR=66 - | когда T меньше 0. или больше π / 2; значение функции полагается равным 3.4E38. |
Версии
| SF20D - | вычисление неполных эллиптических интегралов второго рода с повышенной точностью. |
Вызываемые подпрограммы
| UTSF10 - | подпрограмма выдачи диагностических сообщений при работе функции SF20R. |
| UTSF11 - | подпрограмма выдачи диагностических сообщений при работе функции SF20D. |
Замечания по использованию
| Для функции SF20D параметры RK2, T имеют тип DOUBLE PRECISION. Значение функции SF20D при IERR ≠ 0 полагается равным 1.7D308. |
RK2 = 0.75
T = 3.14159265358 / 3
Y = SF20R (RK2, T, IERR)
Результаты:
Y = 0.918393294308, IERR = 0