Текст подпрограммы и версий
sf15r_c.zip  sf15d_c.zip 
Тексты тестовых примеров
tsf15r_c.zip  tsf15d_c.zip 

Функция:  sf15r_c

Назначение

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

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

Функция sf15r_c производит вычисление следующих интегралов

                  π/2
 (1)      y =  ∫ (1 - x sin2 t) 1/2 dt             для    0 ≤ x ≤ 1
                  0
                  π/2
 (2)      y =  ∫ (1 - x2 sin2 t) 1/2 dt            для    -1 ≤ x ≤ 1
                  0
                  π/2
 (3)      y =  ∫ (1 - (1 - x) sin2 t) 1/2 dt       для     -1 ≤ x ≤ 1
                  0

Алгоритм реализованный в функции, изложен в статье

W.J.Cody, Chebyshev Approximations for the Complete Elliptic Integrals K and E, Math. Comp., 19, 1965.

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

    real sf15r_c (real *x, integer *l, integer *ierr)

Параметры

x - заданное значение аргумента x (тип: вещественный);
l - задает режим работы функции (тип: целый); при этом: если
l = 1 - вычисляется интеграл (1);
l = 2 - вычисляется интеграл (2);
l = 3 - вычисляется интеграл (3);
ierr - целая переменная, служащая для сообщения об ошибках, обнаруженных в ходе работы функции; при этом:
ierr=65 - когда значение l меньше 1 или больше 3; значение функции полагается равным 3.4e38;
ierr=66 - когда значение аргумента меньше 0. При l = 1; значение функции полагается pавным 3.4e38;
ierr=67 - когда значение аргумента больше 1 при l = 1, 2, 3 ; значение функции полагается равным 3.4e38.

Версии

sf15d_c - вычисление полных эллиптических интегралов второго рода с повышенной точностью.

Вызываемые подпрограммы

utsf10_c - подпрограмма выдачи диагностических сообщений при работе функции sf15r_c.
utsf11_c - подпрограмма выдачи диагностических сообщений при работе функции sf15d_c.

Замечания по использованию

 

Область допустимых значений аргумента для функций sf15r_c и sf15d_c указана выше при описании параметра ierr. Для функции sf15d_c аргумент  x  имеет тип двойной точности и при нарушении области допустимых значений функция sf15d_c полагается равной 1.7e308.

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

int main(void)
{
    /* Local variables */
    extern float sf15r_c(float *, int *, int *);
    static int ierr, l;
    static float x, y;

    l = 2;
    x = .33f;
    y = (float)sf15r_c(&x, &l, &ierr);

    printf("\n %16.7e %5i \n",x,l);
    printf("\n %5i %16.7e \n",ierr,y);
    return 0;
} /* main */


Результаты:

       y  =  1.52711603392,     ierr  =  0