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

Назначение

Вычисление неполных эллиптических интегралов первого рода.

Математическое описание

Функция 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