Текст подпрограммы и версий qtsvr_c.zip qtsvd_c.zip |
Тексты тестовых примеров tqtsvr_c.zip tqtsvd_c.zip |
Вычисление неопределенного однократного интеграла от табличной функции, заданной на неравномерной сетке, по формуле типа Симпсона.
Вычисляются значения неопрeдeлeнного интеграла на отрезках [x1, xi] i = 1, ..., N, от табличной функции f (x), заданной на неравномерной сетке xi, i = 1, ..., N, по квадратурной формуле, точной для многочленов второй степени.
Н.С.Бахвалов. Численные методы, "Hаука", M., 1975.
int qtsvr_c (real *rint, real *x, integer *n, real *f)
Параметры
rint - | вещественный вектоp длины n, содержащий вычисленные значения неопределенного интеграла; |
x - | вещественный вектоp длины n, содержащий узлы неравномерной сетки; |
f - | вещественный вектоp длины n, содержащий значения функции f (x); |
n - | заданное число узлов сетки (тип: целый). |
Версии
qtsvd_c - | вычисление с удвоенной точностью неопределенного однократного интеграла от табличной функции, заданной на неравномерной сетке, по формуле типа Симпсона. |
Вызываемые подпрограммы: нет
Замечания по использованию
Требуется, чтобы xi > xi - 1, i = 2, ..., n; n ≥ 3 . В подпрограмме qtsvd_c параметры rint, x, f имеют тип double. |
int main(void) { /* Builtin functions */ double sin(double); /* Local variables */ static float rint[50], f[50], h__; static int i__, n; static float x[50], x1; extern int qtsvr_c(float *, float *, int *, float *); int i__1; h__ = .064114081632653058f; x[0] = 0.f; x[1] = h__ / 4.f; n = 50; i__1 = n; for (i__ = 3; i__ <= i__1; ++i__) { /* l1: */ x[i__ - 1] = x[i__ - 3] + h__ * 2.f; } i__1 = n; for (i__ = 1; i__ <= i__1; ++i__) { x1 = x[i__ - 1]; /* l2: */ f[i__ - 1] = (float)sin(x1); } qtsvr_c(rint, x, &n, f); for (i__ = 0; i__ <= 45; i__+= 5) { printf("\n %11.7f %11.7f %11.7f %11.7f %11.7f \n", rint[i__], rint[i__+1], rint[i__+2], rint[i__+3], rint[i__+4]); } return 0; } /* main */ Результаты: rint(1) = 0.00000 rint(2) = 0.00012 rint(49) = 1.99794 rint(50) = 1.99884