|
Текст подпрограммы и версий ( Фортран ) qtt8r.zip qtt8d.zip |
Тексты тестовых примеров ( Фортран ) tqtt8r.zip tqtt8d.zip |
|
Текст подпрограммы и версий ( Си ) qtt8r_c.zip qtt8d_c.zip |
Тексты тестовых примеров ( Си ) tqtt8r_c.zip tqtt8d_c.zip |
|
Текст подпрограммы и версий ( Паскаль ) qtt8r_p.zip qtt8e_p.zip |
Тексты тестовых примеров ( Паскаль ) tqtt8r_p.zip tqtt8e_p.zip |
Вычисление определенного двухкратного интеграла от табличной функции, заданной на неравномерной сетке, по формуле трапеций.
По формуле трапеций вычисляется значение определенного двойного интеграла от табличной функции f (x, y), заданной на неравномерной сетке узлов (xi, yj), i = 1, ..., N, j = 1, ..., M .
Н.С.Бахвалов. Численные методы, "Hаука", M., 1975.
SUBROUTINE QTT8R (RINT, X1, X2, F, N, M)
Параметры
| RINT - | вещественная переменная, содержащая вычисленное значение интеграла; |
| X1 - | вещественный вектоp длины N, содержащий узлы неравномерной сетки по x; |
| X2 - | вещественный вектоp длины M, содержащий узлы неравномерной сетки по y; |
| F - | вещественный двумерный массив размера N на M, содержащий значения функции f (x, y); |
| N - | заданное число узлов сетки по x (тип: целый); |
| M - | заданное число узлов сетки по y (тип: целый). |
Версии
| QTT8D - | вычисление с удвоенной точностью определенного двухкратного интеграла от табличной функции, заданной на неравномерной сетке, по формуле трапеций. |
Вызываемые подпрограммы: нет
Замечания по использованию
|
Требуется, чтобы N ≥ 2, M ≥ 2. В подпрограмме QTT8D параметры RINT, X1, X2, F имеют тип DOUBLE PRECISION. |
DIMENSION F(50, 40), X1(50), X2(40)
H1 = 3.14159 / 49.
H2 = 3.14159 / 39.
N = 50
M = 40
X1(1) = 0.
X2(1) = 0.
X1(2) = H1 / 4.
X2(2) = H2 / 4.
DO 1 I = 3, N
1 X1(I) = X1(I - 2) + 2.*H1
DO 2 J = 3, M
2 X2(J) = X2(J - 2) + 2.*H2
DO 3 I = 1, N
DO 3 J = 1, M
X = X1(I) + X2(J)
3 F(I, J) = SIN(X)
CALL QTT8R (RINT, X1, X2, F, N, M)
Результат:
RINT = 0.21622