Текст подпрограммы и версий ( Фортран )
sf84r.zip  sf84d.zip  		Тексты тестовых примеров ( Фортран )
tsf84r.zip  tsf84d.zip 
Текст подпрограммы и версий ( Си )
sf84r_c.zip  sf84d_c.zip  		Тексты тестовых примеров ( Си )
tsf84r_c.zip  tsf84d_c.zip 
Текст подпрограммы и версий ( Паскаль )
sf84r_p.zip , sf84e_p.zip 		Тексты тестовых примеров ( Паскаль )
tsf84r_p.zip , tsf84e_p.zip

Функция:  SF84R

Назначение

Вычисление полного эллиптического интеграла  cel (kc, p, a, b) второго рода общего вида

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

Подпрограмма - функция SF84R вычисляет полный эллиптический интеграл второго рода общего вида: 

                               ∞
cel( kc, p, a, b)  =   ∫   [(a + b t 2) dt] / [ (1 + p t 2)  ( (1 + t 2)(1 + kc2 t 2) ) 1/2 ] ,
                              0

где  kc,  p,  a,  b - вещественные числа

Справочник по специальным функциям. Под ред. Абрамовица М. и Стиган И., Изд - во "Наука", 1979.

Использование 

    FUNCTION  SF84R (RKC, P, A, B) 
   Параметры
       RKC, P -
         A, B   		заданные значения аргументов  kc,  p,  a  и  b соответственно (тип: вещественный) 
   Версии
SF84D - вычисление полного эллиптического интеграла  cel (kc, p, a, b) второго рода общего вида в режиме удвоенной точности; при этом все параметры функции SF84D должны иметь тип DOUBLE PRECISION, а сама функция должна быть описана как DOUBLE PRECISION  FUNCTION 
   Вызываемые подпрограммы:  нет 
   Замечания по использованию
 

Если RKC = 0 и  B = 0, то значение функции полагается равным 1.

Если RKC = 0 и  B ≠ 0, то значение функции полагается равным наибольшему представимому на используемой машине числу (это означает, что при данном условии интеграл расходится)

Пример использования 

       RKC = 0.5 
       P = 1.0 
       A = 1.0 
       B = 1.0 
       RES = SF84R (RKC, P, A, B) 

Результат:    RES = 2.15652