Текст подпрограммы и версий ( Фортран ) qsk3r.zip qsk3d.zip |
Тексты тестовых примеров ( Фортран ) tqsk3r.zip tqsk3d.zip |
Текст подпрограммы и версий ( Си ) qsk3r_c.zip qsk3d_c.zip |
Тексты тестовых примеров ( Си ) tqsk3r_c.zip tqsk3d_c.zip |
Текст подпрограммы и версий ( Паскаль ) qsk3r_p.zip qsk3e_p.zip |
Тексты тестовых примеров ( Паскаль ) tqsk3r_p.zip tqsk3e_p.zip |
Вычисление определенного интеграла от векторной функции по обобщенной квадратурной формуле Ньютона - Котесса шестого порядка точности.
Подпрограмма QSK3R вычисляет интеграл B ∫ F(x) dx от векторной функции F(x) = ( f1 (x), ... , fN (x) ) A с покомпонентной погрешностью B ЕРS ( 1 + | ∫ fi (x) dx | ) A
по обобщенной квадратурной формуле Ньютона - Котеса шестого порядка точности с шагом h = (B - A)/M, где M - число частичных отрезков разбиения.
Я.М.Жилейкин, М.В.Соколовский. Набор стандартных программ для вычисления интегралов от векторных функций. Сб. "Методы и алгоритмы в численном анализе", Изд-во МГУ, 1981.
SUBROUTINE QSK3R (RINT, A, B, F, N, EPS, K, IND, R, IERR)
Параметры
RINT - | вещественный вектоp длины N, содержащий вычисленные значения интегралов; |
A, B - | заданные нижний и верхний пределы интегрирования (тип: вещественный); |
F - | имя подпрограммы, вычисляющей подинтегральные функции; |
N - | размерность вектоpа - функции F (x) (тип: целый); |
EPS - | заданная меpа погрешности вычисления интеграла (тип: вещественный); |
K - | целая переменная, задающая начальное число частичных отрезков разбиения; |
IND - | целый вектоp длины N, каждая компонента которого pавна числу частичных отрезков разбиения, при котоpом достигается заданная точность; |
R - | вещественный рабочий вектоp длины 4N; |
IERR - | целая переменная, служащая для диагностических сообщений: |
IERR = 65 - | когда заданная точность не может быть достигнута при максимально возможном числе (1048576) частичных отрезков разбиения. |
Версии
QSK3D - | вычисление с удвоенной точностью определенного интеграла от векторной функции по обобщенной квадратурной формуле Ньютона - Котесса шестого порядка точности. |
Вызываемые подпрограммы
UTQS10 - | подпрограмма выдачи диагностических сообщений при работе подпрограммы QSK3R. |
UTQS12 - | подпрограмма выдачи диагностических сообщений при работе подпрограммы QSK3D. |
Замечания по использованию
1. |
Первый оператор подпрограммы F должен иметь вид: | |
2. |
При EPS ≤ 0 происходит только один просчет при заданном значении K. | |
3. |
Вычисление интеграла от отдельной компоненты прекращается как только на этой компоненте достигается заданная точность. | |
4. |
Если для каких - либо компонент векторной функции заданная погрешность интегрирования не может быть достигнута, то соответствующие компоненты вектоpа IND полагаются равными нулю. | |
5. |
Погрешность определяется по модулю разности двух просчетов по M и 2M частичным отрезкам разбиения. Если заданная точность не достигнута, то значение M удваивается. | |
6. | В подпрограмме QSK3D параметры RINT, A, B, F, EPS, R имеют тип DOUBLE PRECISION. |
SUBROUTINE FUN(X, Y, I) DIMENSION Y(5) Y = X**3 DO 1 J = 2, 5 1 Y(J) = Y(J-1)*X RETURN END DIMENSION RINT(5), R(20), IND(5) EXTERNAL FUN CALL QSK3R (RINT, 0., 1., FUN, 5, 1.E-6, 1, IND, R, IERR) Результаты: RINT IND 2.5000000000-01 1 1.9999999999-01 1 1.6666666666-01 1 1.4285714427-01 4 1.2500000496-01 4 IERR = 0