Текст подпрограммы и версий ( Фортран )
qts2r.zip  qts2d.zip  		Тексты тестовых примеров ( Фортран )
tqts2r.zip  tqts2d.zip 
Текст подпрограммы и версий ( Си )
qts2r_c.zip  qts2d_c.zip  		Тексты тестовых примеров ( Си )
tqts2r_c.zip  tqts2d_c.zip 
Текст подпрограммы и версий ( Паскаль )
qts2r_p.zip  qts2e_p.zip  		Тексты тестовых примеров ( Паскаль )
tqts2r_p.zip  tqts2e_p.zip 

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

Назначение

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

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

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

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

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

    SUBROUTINE  QTS2R (RINT, X, F, N) 
   Параметры
RINT - вещественная переменная, содержащая вычисленное значение интеграла;
X - вещественный вектоp длины N, содержащий узлы неравномерной сетки;
F - вещественный вектоp длины N, содержащий значения функции f (x);
N - заданное число узлов сетки (тип: целый). 
   Версии
QTS2D - вычисление с удвоенной точностью определенного однократного интеграла от табличной функции, заданной на неравномерной сетке, по формуле Симпсона. 
   Вызываемые подпрограммы:  нет 
   Замечания по использованию
 

Требуется, чтобы   xi > xi - 1,  i = 2, ..., N;   N ≥ 3 .

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

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

          DIMENSION F(50), X(50)
          H = 3.14159 / 49.
          N = 50
          X(1) = 0.
          X(2) = H / 4.
          DO 1 I = 3, N
       1 X(I) = X(I - 2) + 2.*H
          DO 2 I = 1, N
          X1 = X(I)
       2 F(I) = SIN(X1)
          CALL  QTS2R (RINT, X, F, N)

Результат:

      RINT  =  1.99884