БЧА НИВЦ МГУ. SF68R. Цилиндрические функции

Документ взят из кэша поисковой машины. Адрес оригинального документа : http://num-anal.srcc.msu.su/lib_na/cat/sf_htm_c/sf68r_c.htm
Дата изменения: Thu Mar 12 10:01:26 2015
Дата индексирования: Sun Apr 10 01:49:43 2016
Кодировка: Windows-1251

Текст подпрограммы и версий
sf68r_c.zip

Тексты тестовых примеров
tsf68r_c.zip

Функция: sf68r_c

Назначение

Вычисление вещественной и мнимой части модифицированной функции Бесселя второго рода мнимого индекса и вещественного аргумента K_{1/2 + i β} (x).

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

sf68r_c вычисляет вещественную и мнимую части модифицированной функции Бесселя K_{1/2 + i β} (x).

B зависимости от значений аргумента и индекса используется метод, основанный на степенных разложениях, на интегральных представлениях или τ - метод Ланцоша.

Раппопорт Ю.М. Таблицы модифицированных функций Бесселя K_{1/2 + i β} (x), Hаука, 1979.

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

    int sf68r_c (real *x, real *r, real *t1, real *t2)

   Параметры

x -	заданное значение аргумента x (тип: вещественный);
r -	заданное значение мнимой части индекса β (тип: вещественный);
t1 -	вещественная переменная, значение которой полагается равным результату вычисления вещественной части K_{1/2 + i β} (x);
t2 -	вещественная переменная, значение которой полагается равным результату вычисления мнимой части K_{1/2 + i β} (x).

   Версии:  нет

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

sf23c_c -	подпрограмма вычисления значений гамма - функции и ее натурального логарифма от комплексного аpгумента.
qsf3r_c -	подпрограмма вычисления определенных интегралов от осциллирующих функций по формулам интерполяционного типа пятой степени точности с автоматическим выбором шага для больших отрезков интегрирования.

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

Для подпрограммы sf68r_c областями допустимых значений параметров x и β являются 0.1 ≤ x ≤ 10., 0.01 ≤ β ≤ 10.

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

int main(void)
{
    /* Initialized data */
    static float x = 1.f;
    static float r__ = 1.f;

    /* Local variables */
    extern int sf68r_c(float *, float *, float *, float *);
    static float t1, t2;

    sf68r_c(&x, &r__, &t1, &t2);

    printf("\n %16.7e %16.7e \n",t1,t2);
    return 0;
} /* main */


Результаты:

       t1  =  0.298824989 , 
       t2  =  0.118944694