Тема: задача з інтегралами на сі
#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)