حل معادله به روش سیمپسون

# include "stdio.h"
# include "conio.h"
main()
{
float   a,b,h,i1,h2,n,i2,x,i3,i;
float   w1,w2,w3,k;
int j;
clrscr();
printf ("enter 1 number for A:");
scanf ("%f",&a);
printf ("enter 1 number for B:");
scanf ("%f",&b);
printf ("enter 1 number for eps:");
scanf ("%f",&h);
n=(b-a)/ h;
printf("n:= %f",n);
h2=h*2;
i1=(1/(1+a*a))+(1/(1+b*b));
i2=0;
i3=0;
x=a+h;
j=1;
while (j<=(n/2))
{
i2=i2+(1/(1+x*x));
x=x+h2;
j=j+1;
}
x=a;
j=1;
while (j<=((n/2)-1))
{
x=x+h2;
i3=i3+(1/(1+x*x));
j=j+1;
}
w1=h/3;
w2=4*i2;
w3=2*i3;
i=(w1*(i1+w2+w3));
printf ("\n               (f(1)  +   f(2)    +    f(3)   +    f(4)   +   f(5))");
printf("\n---------------------------------------------------------------------------");
printf ("\ns(h)=(h/3)*(f(%f",a,")");
for (j=1;j<=4;j++)
{
printf ("+");
k=j*h;
printf("f(%f",k);
printf(")");
}
printf(")=%f",i);
getch();
}