|
Текст подпрограммы и версий ( Фортран ) sf83r.zip sf83d.zip |
Тексты тестовых примеров ( Фортран ) tsf83r.zip tsf83d.zip |
|
Текст подпрограммы и версий ( Си ) sf83r_c.zip sf83d_c.zip |
Тексты тестовых примеров ( Си ) tsf83r_c.zip tsf83d_c.zip |
|
Текст подпрограммы и версий ( Паскаль ) sf83r_p.zip , sf83e_p.zip |
Тексты тестовых примеров ( Паскаль ) tsf83r_p.zip , tsf83e_p.zip |
Вычисление неполного эллиптического интеграла el2 (x, kc, a, b) второго рода общего вида
Подпрограмма - функция SF83R вычисляет неполный эллиптический интеграл второго рода общего вида:
x
el2(x, kc, a, b) = ∫ [(a + b t 2) dt] / [ (1 + t 2) ( (1 + t 2)(1 + kc2 t 2) ) 1/2 ] ,
0
где x ≥ 0, kc, a, b - вещественные числа
Справочник по специальным функциям. Под ред. Абрамовица М. и Стиган И., Изд - во "Наука", 1979.
FUNCTION SF83R (X, RKC, A, B)
Параметры
|
X, RKC - A, B | заданные значения аргументов x, kc, a и b соответственно (тип: вещественный) |
Версии
| SF83D - | вычисление неполного эллиптического интеграла el2 (x, kc, a, b) второго рода общего вида в режиме удвоенной точности; при этом все параметры функции SF83D должны иметь тип DOUBLE PRECISION, а сама функция должна быть описана как DOUBLE PRECISION FUNCTION |
Вызываемые подпрограммы: нет
Замечания по использованию: нет
X = 3.14159265 / 3.0
RKC2 = 0.25
RKC = SQRT(RKC2)
A = 1.0
RES = SF83R (X, RKC, A, RKC2)
Результат: RES = 0.746321