Текст подпрограммы и версий ( Фортран )
sf17r.zip  sf17d.zip  		Тексты тестовых примеров ( Фортран )
tsf17r.zip  tsf17d.zip 
Текст подпрограммы и версий ( Си )
sf17r_c.zip  sf17d_c.zip  		Тексты тестовых примеров ( Си )
tsf17r_c.zip  tsf17d_c.zip 
Текст подпрограммы и версий ( Паскаль )
sf17r_p.zip , sf17e_p.zip 		Тексты тестовых примеров ( Паскаль )
tsf17r_p.zip , tsf17e_p.zip

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

Назначение

Вычисление производных функций Кельвина нулевого порядка.

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

SF17R вычисляет производные функций Кельвина первого рода (ber' (x) и bei' (x)) и второго рода (ker' (x) и kei' (x)) нулевого порядка.

J.F.Hart, E.W.Cheney, C.L.Lawson et al, Computer Approximations, Wiley, New York, 1968.

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

    SUBROUTINE  SF17R (X, BERD, BEID, AKERD, AKEID, IERR) 
   Параметры
X - заданное значение аргумента x (тип: вещественный);
           BERD -
           BEID  
         AKERD  
         AKEID   		вычисленные значения функций ber' (x), bei' (x), ker' (x), kei' (x) соответственно (тип: вещественный);
IERR - целая переменная, служащая для сообщения об ошибках, обнаруженных в ходе работы подпрограммы; при этом:
IERR= 1 - когда значение аргумента меньше 0.; значения AKERD и AKEID полагаются равными 3.4E38;
IERR=66 - когда абсолютное значение аргумента больше 57.58; значения BERD и BEID полагаются равными 0.; в случае, когда аргумент неотрицателен, AKERD и AKEID также полагаются равными 0.; иначе AKERD и AKEID полагаются равными 3.4E38. 
   Версии
SF17D - вычисление производных функций Кельвина нулевого порядка. 
   Вызываемые подпрограммы
       SF16R -
       SF16D   		подпрограммы вычисления функций Кельвина нулевого порядка с обычной точностью и повышенной точностью ( используются в SF17R и SF17D соответственно).
UTSF10 - подпрограмма выдачи диагностических сообщений при работе подпрограммы SF17R.
UTSF11 - подпрограмма выдачи диагностических сообщений при работе подпрограммы SF17D. 
   Замечания по использованию
 

Для подпрограммы SF17D параметры X, BERD, BEID, AKERD и AKEID имеют тип DOUBLE PRECISION и абсолютное значение аргумента X не должно превосходить 2.04E03. При нарушении этого условия значение IERR полагается равным 66 и для отрицательного X значения AKERD и AKEID полагаются равными 1.7D308.

При X = 0, BERD = 0., BEID = 0., AKEID = 0., AKERD = - 3.4E38 (для SF17R) или - 1.7D308 (для SF17D).

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

         X = 0.6
         CALL  SF17R (X, BERD, BEID, AKERD, AKEID, IERR)

Результаты:

       BERD     =  -0.134984812705, 
       BEID      =    0.299797506834, 
       AKERD  =  -1.45653855074, 
       AKEID   =    0.348164425113, 
       IERR      =    0