|
Текст подпрограммы и версий qsl1r_c.zip qsl1d_c.zip |
Тексты тестовых примеров tqsl1r_c.zip tqsl1d_c.zip |
Вычисление определенного интеграла по формулам Лобатто.
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