Не работает метод трапеций с функцией - C (СИ)

#include <conio.h>
#include <stdio.h>
#include <math.h>
 
float func (float x)
{
float SF;
SF=pow(sin(x),2)/x;
return SF;
}
int main()
{
int n,i,k=0;
float x0,x,y,a,b,eps,h,fa,fb,S,S1,x1,x2,SF,fx1,fx2;
printf("vvedie 'a'\n'a' = "); scanf("%f",&a);
printf("vvedie 'b'\n'b' = "); scanf("%f",&b);
printf("vvedie 'eps'\n'eps' = "); scanf("%f",&eps);
printf("vvedie 'n'\n'n' = "); scanf("%f",&n);
fa=func(a);
fb=func(b);
S=0.5*(fa+fb)*(b-a);
h=(b-a)/n;
do
{
x1=a; S1=0; k++;
for (i=0;i<n;i++)
{
fx1=func(x1);
x2=x1+h;
fx2=func(x2);
S1+=0.5*(fx1+fx2)*h;
x1+=h;
}
n*=2; h=(b-a)/n;
}
while (fabs(S1-S)>eps);
printf("S = %f\n",S1);
printf("n = %d\n",k);
getch();
return 0;
}