Текст подпрограммы и версий ( Фортран )
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

Назначение

Суммирование знакопеременных рядов при помощи преобразования Эйлера.

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

Подпрограмма 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