#include
#include

float trapez(float a, float b, float (*funkcja)(float x), int n);

float fun(float x)
{
return 1./(1.+x*x);
}
void main()
{
int n,m,k;
float a, b, canal, blt ,cnt;
FILE *plik;
plik=fopen("wynik.dat","w");

printf("podaj k:\n");
scanf("%d",&k);
printf("podaj a, b:\n");
scanf("%f%f",&a,&b);

for(m=1;m<=k;m++)
{
n=pow(2,m);
cnt=trapez(a,b,fun,n);
canal=atan(b)-atan(a);
blt=fabs(cnt-canal);
fprintf(plik,"%d %f\n",n,blt);

}


}

float trapez( float a, float b, float(* funkcja) (float x) , int n )
{
float value;
/* Calkowanie metoda trapezow */
float h, x, sum;
int i;

h = (b - a) / (float) n; // krok calkowania

sum = 0.0;
for( i=1; i x = a + i * h;
sum += funkcja(x);
}

value = ( funkcja(a) + 2.0*sum + funkcja(b) ) * h / 2.0;

return value;
}