Текст подпрограммы и версий
qtsvr_c.zip  qtsvd_c.zip 
Тексты тестовых примеров
tqtsvr_c.zip  tqtsvd_c.zip 

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

Назначение

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

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

Вычисляются значения неопр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