1

Тема: задача з інтегралами на сі

#include <stdio.h>
#include <math.h>

const float E = 1.e-3;

float fn1(float x) {
    return sqrt(1 + log(x));
}
float fn2(float x) {
    return log(1 + pow(x, 2)) / (1 + pow(x, 2));
}
float ft(int n, float a, float b, float fun(float)) {
    int i;
    float s1, h, s=0;
    do
    {
        s1 = s;
        h = (b - a) / n;
        s = (fun(a) + fun(b)) / 2;
        for(int i = 1; i <= n - 1; i++)
        s += fun(a + i * h);
        s * = h;
        n * = 2;
    }
    while(fabs(s - s1) > E);
    return s;
}
int main () {
    double y;
    y = ft(20, 2.2, 3.0, fn1) + ft(20, 0, 1.0, fn2);
    printf("%.5f \n", y);
    return 0;
}

Скласти програму з використанням функції обчислення інтегралів методом трапецій (точність обчислення е = 10-3)
http://cpp.dp.ua/uploads/posts/2015-12/1450287731_9_4.png