Текст подпрограммы и версий ( Фортран )
sf85r.zip  sf85d.zip  		Тексты тестовых примеров ( Фортран )
tsf85r.zip  tsf85d.zip 
Текст подпрограммы и версий ( Си )
sf85r_c.zip  sf85d_c.zip  		Тексты тестовых примеров ( Си )
tsf85r_c.zip  tsf85d_c.zip 
Текст подпрограммы и версий ( Паскаль )
sf85r_p.zip , sf85e_p.zip 		Тексты тестовых примеров ( Паскаль )
tsf85r_p.zip , tsf85e_p.zip

Подпрограмма:  SF85R

Назначение

Вычисление эллиптических функций Якоби  sn (u, kc2),  cn (u, kc2) и  dn (u, kc2)

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

Подпрограмма SF85R вычисляет значения эллиптических функций Якоби  sn (u, kc2),  cn (u, kc2) и  dn (u, kc2), которые удовлетворяют следующим уравнениям: 

                           sn
                   u  =   ∫   dt  /  [ 1 - (1 - kc2) sin2 t ] 1/2
                           0                       
                 sn2 + cn2  =  1  ,      (1 - kc2) sn2 + dn2  =  1

Аргументы  u  и  kc принимают вещественные значения

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

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

    SUBROUTINE  SF85R (U, RKC2, SN, CN, DN) 
   Параметры
            U -
         RKC2   		заданные значения аргументов  u  и  kc2 соответственно (тип: вещественный);
       SN, CN -
         DN   		вещественные переменные, значения которых полагаются равными значениям функций  sn (u, kc2),  cn (u, kc2) и  dn (u, kc) 
   Версии
SF85D - вычисление эллиптических функций Якоби  sn (u, kc2),  cn (u, kc2) и  dn (u, kc2) в режиме удвоенной точности; при этом все параметры должны иметь тип DOUBLE PRESICION 
   Вызываемые подпрограммы:  нет 
   Замечания по использованию:  нет

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

       U = 0.61802 
       RKC2 = 0.5 
       CALL  SF85R (U, RKC2, SN, CN, DN) 

Результаты:

       SN = 0.564576 
       CN = 0.825381 
       DN = 0.916857