Текст подпрограммы и версий
qsl1r_c.zip  qsl1d_c.zip 
Тексты тестовых примеров
tqsl1r_c.zip  tqsl1d_c.zip 

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

Назначение

Вычисление определенного интеграла по формулам Лобатто.

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

qsl1r_c вычисляет интеграл
            B
           ∫ f (x) dx 
          A 

с заданной абсолютной погрешностью E по формулам Лобатто по семи узлам (формула 11 - й степени точности) с автоматическим выбором шага.

О.В.Брушлинская, Л.Г.Васильева. Набор стандартных программ приближенного вычисления однократных интегралов с автоматическим выбором шага. Сб. "Численный анализ на ФОРТРАНе", вып. 8, Изд-во МГУ, 1974.

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

    int qsl1r_c (real *rint, real *a, real *b, R_fp f, real *e)

Параметры

rint - переменная, содержащая вычисленное значение интеграла (тип: вещественный);
a, b - заданные нижний и верхний пределы интегрирования (тип: вещественный);
f - имя вещественной подпрограммы - функции вычисления подинтегральной функции  f (x);
e - заданная абсолютная погрешность вычисления интеграла (тип: вещественный).

Версии

qsl1d_c - вычисление с удвоенной точностью определенного интеграла по формулам Лобатто.

Вызываемые подпрограммы: нет

Замечания по использованию

 

Подпрограмма не гарантирует заданную погрешность, хотя в большинстве случаев интегралы будут вычислены с тpебуемой точностью.

В подпрограмме qsl1d_c параметры rint, a, b, f, e имеют тип double.

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

int main(void)
{
    /* Local variables */
    static float rint, a, b;
    extern int qsl1r_c(float *, float *, float *, R_fp, float *);
    static float e;
    extern float f_c();

    a = 0.f;
    b = 1.f;
    e = 1e-8f;
    qsl1r_c(&rint, &a, &b, (R_fp)f_c, &e);

    printf("\n %16.7e \n",rint);
    return 0;
} /* main */

float f_c(float *x)
{
    /* System generated locals */
    float ret_val;

    /* Builtin functions */
    double sqrt(double);

    if (*x != 0.f) {
        goto l1;
    } else {
        goto l2;
    }
l2:
    ret_val = 0.f;
    goto l3;
l1:
    ret_val = 1.f / (float)sqrt(*x);
l3:
    return ret_val;
} /* f_c */


Результат:   rint  =  1.9999963224