|
Документ взят из кэша поисковой машины. Адрес
оригинального документа
: http://www.srcc.msu.ru/num_anal/lib_na/cat/sf/sf50r.htm
Дата изменения: Mon Mar 14 11:07:48 2011 Дата индексирования: Mon Feb 4 02:15:28 2013 Кодировка: Windows-1251 |
|
Текст подпрограммы и версий ( Фортран ) sf50r.zip sf50d.zip |
Тексты тестовых примеров ( Фортран ) tsf50r.zip tsf50d.zip |
|
Текст подпрограммы и версий ( Си ) sf50r_c.zip sf50d_c.zip |
Тексты тестовых примеров ( Си ) tsf50r_c.zip tsf50d_c.zip |
|
Текст подпрограммы и версий ( Паскаль ) sf50r_p.zip , sf50e_p.zip |
Тексты тестовых примеров ( Паскаль ) tsf50r_p.zip , tsf50e_p.zip |
Вычисление полных эллиптических интегралов третьего рода П (x, k).
Функция SF50R производит вычисление интегралов вида:
p/2
П (x, k) = т dt / (1 + x sin 2 t) (1 - k 2 sin 2 t) 1 / 2
0
для 0 ? k 2 ? 1 и x ? - 1
J.F.Hart, E.W.Cheney, C.L.Lawson et al, Computer Approximations, Wiley, New York, 1968.
REAL FUNCTION SF50R (X, RK2, IERR)
Параметры
| X, RK2 - | заданные значения аргументов x и k2 соответственно (тип: вещественный); |
| IERR - | целая переменная, служащая для сообщения об ошибках, обнаруженных в ходе работы функции; при этом: |
| IERR=65 - | когда значение X pавно - 1; значение функции полагается равным 3.4E38; |
| IERR=66 - | когда значение RK2 меньше 0 или больше или pавно 1; значение функции полагается равным 3.4E38. |
Версии
| SF50D - | вычисление полных эллиптических интегралов третьего рода с повышенной точностью. |
Вызываемые подпрограммы
| SF14R - SF14D | вычисление полных эллиптических интегралов первого рода с обычной и повышенной точностью (используются в SF50R и SF50D соответственно). |
| SF15R - SF15D | вычисление полных эллиптических интегралов второго рода с обычной и повышенной точностью (используются в SF50R и SF50D соответственно). |
| SF19R - SF19D | вычисление неполных эллиптических интегралов первого рода с обычной и повышенной точностью (используются в SF50R и SF50D соответственно). |
| SF20R - SF20D | вычисление неполных эллиптических интегралов второго рода с обычной и повышенной точностью (используются в SF50R и SF50D соответственно). |
| UTSF14 - | подпрограмма выдачи диагностических сообщений при работе функции SF50R. |
| UTSF15 - | подпрограмма выдачи диагностических сообщений при работе функции SF50D. |
Замечания по использованию
| Для функции SF50D параметры X и RK2 имеют тип DOUBLE PRECISION. Значение SF50D при IERR ? 0 полагается pавным 1.7D308. |
X = 1.
RK2 = 0.5
Y = SF50R (X, RK2, IERR)
Результаты: Y = 1.27312736674 , IERR = 0