Текст подпрограммы и версий ( Фортран ) rcr1r.zip , rcr1d.zip |
Тексты тестовых примеров ( Фортран ) trcr1r.zip , trcr1d.zip |
Текст подпрограммы и версий ( Си ) rcr1r_c.zip , rcr1d_c.zip |
Тексты тестовых примеров ( Си ) trcr1r_c.zip , trcr1d_c.zip |
Текст подпрограммы и версий ( Паскаль ) rcr1r_p.zip , rcr1e_p.zip |
Тексты тестовых примеров ( Паскаль ) trcr1r_p.zip , trcr1e_p.zip |
Суммирование знакопеременных рядов при помощи преобразования Эйлера.
Подпрограмма RCR1R накапливает частичную сумму SUM знакопеременного ряда по заданному текущему члену TERM ряда и его номеру JTERM. При первом обращении к подпрограмме на вход подается первый член ряда (JTERM = 1), при втором обращении на вход подается второй член ряда (JTERM = 2) со знаком, противоположным знаку первого члена, и т.д. При суммировании используется преобразование Эйлера.
Справочник по специальным функциям. Под ред. Абрамовица М. и Стигана И. Изд - во "Наука", 1979.
SUBROUTINE RCR1R (SUM, TERM, JTERM, RAB)
Параметры
SUM - | вещественная переменная, содержащая на выходе из подпрограммы текущее значение частичной суммы знакопеременного ряда; |
TERM - | значение члена ряда с номером JTERM (тип: вещественный); |
JTERM - | значение номера текущего члена ряда (тип: целый); |
RAB - | вещественный рабочий вектор, длина которого равна максимальному количеству членов ряда, предназначенных для суммирования. |
Версии
RCR1D - | суммирование знакопеременных рядов при помощи преобразования Эйлера в режиме удвоенной точности; при этом параметры SUM, TERM и RAB должны иметь тип DOUBLE PRECISION. |
Вызываемые подпрограммы: нет
Замечания по использованию
Данная подпрограмма может быть применена не только для сходящихся знакопеременных рядов. Ее можно применять для суммирования асимптотических рядов. |
В данном примере суммируются 10 членов ряда:
ln (1 + x) = x - x2/2 + x3/3 - ... , | x | ≤ 1 , x ≠ - 1 DIMENSION RAB(10) X = 0.5 DO 1 K = 1, 10 TERM = (- 1)**(K - 1)*X**K/K JTERM = K CALL RCR1R (SUM, TERM, JTERM, RAB) 1 CONTINUE Результат: SUM = 0.4054769