Текст подпрограммы и версий
sf85r_c.zip  sf85d_c.zip  		Тексты тестовых примеров
tsf85r_c.zip  tsf85d_c.zip 

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

Назначение

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

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

Подпрограмма sf85r_c вычисляет значения эллиптических функций Якоби  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.

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

    int sf85r_c (real *u, real *rkc2, real *sn, real *cn, real *dn)

   Параметры
u -
rkc2   		заданные значения аргументов  u  и  kc2 соответственно (тип: вещественный);
sn, cn -
dn   		вещественные переменные, значения которых полагаются равными значениям функций  sn (u, kc2),  cn (u, kc2) и  dn (u, kc) 
   Версии
sf85d_c - вычисление эллиптических функций Якоби  sn (u, kc2),  cn (u, kc2) и  dn (u, kc2) в режиме удвоенной точности; при этом все параметры должны иметь тип double. 
   Вызываемые подпрограммы:  нет 
   Замечания по использованию:  нет

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

int main(void)
{
    /* Local variables */
    extern int sf85r_c(float *, float *, float *, float *, float *);
    static float u, cn, dn, sn, rkc2;

    u = .61802f;
    rkc2 = .5f;
    sf85r_c(&u, &rkc2, &sn, &cn, &dn);

    printf("\n %16.7e %16.7e %16.7e \n",sn,cn,dn);
    return 0;
} /* main */


Результаты:

       sn = 0.564576 
       cn = 0.825381 
       dn = 0.916857