Текст подпрограммы и версий ( Фортран )
qtscr.zip  qtscd.zip 
Тексты тестовых примеров ( Фортран )
tqtscr.zip  tqtscd.zip 
Текст подпрограммы и версий ( Си )
qtscr_c.zip  qtscd_c.zip 
Тексты тестовых примеров ( Си )
tqtscr_c.zip  tqtscd_c.zip 
Текст подпрограммы и версий ( Паскаль )
qtscr_p.zip  qtsce_p.zip 
Тексты тестовых примеров ( Паскаль )
tqtscr_p.zip  tqtsce_p.zip 

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

Назначение

Вычисление неопределенного однократного интеграла от табличной функции, заданной на равномерной сетке, по формуле Симпсона.

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

Вычисляются значения неопределенного интеграла на отрезках [x1, xi]  i = 1, ..., N, от табличной функции f (x), заданной на равномерной сетке xi = x1 + (i - 1) h,  i = 1, ..., N, по квадратурной формуле, точной для многочленов второй степени.

Н.С.Бахвалов. Численные методы, "Hаука", M., 1975.

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

    SUBROUTINE  QTSCR (RINT, H, F, N) 

Параметры

RINT - вещественный вектоp длины N, содержащий вычисленные значения неопределенного интеграла;
H - заданный шаг равномерной сетки (тип: вещественный);
F - вещественный вектоp длины N, содержащий значения функции f (x);
N - заданное число узлов сетки (тип: целый).

Версии

QTSCD - вычисление с удвоенной точностью неопределенного однократного интеграла от табличной функции, заданной на равномерной сетке, по формуле Симпсона.

Вызываемые подпрограммы: нет

Замечания по использованию

 

Требуется, чтобы  N ≥ 3.

В подпрограмме QTSCD параметры RINT, H, F имеют тип DOUBLE PRECISION.

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

          DIMENSION  F(50), X(50), RINT(50)
          H = 3.14159 / 49.
          N = 50
          X = 0.
          DO 1 I = 1, N
          F(I) = SIN(X)
       1 X = X + H
          CALL  QTSCR (RINT, H, F, N)

Результаты:

       RINT(1)    =  0.00000
       RINT(2)    =  0.00205
       RINT(49)  =  1.99794
       RINT(50)  =  2.00000