Вычисление определенного интеграла методами трапеции, Симпсона (парабол), средних прямоугольников в borland си

#include <stdio.h>
#include <stdlib.h>
#include <conio.h>
#include <math.h>
main()
{
const int n=100;
double a=0, b=pow(3,(1/2));
double atan(x),x,h;
double F[n],I1,I2,I3,chet,nechet;
int i;
h=(b-a)/n;
//Трапеции
for(i=0,x=a; i<=n; i++,x=x+h)
 F[i] = x*atan(x);
I1 = 0;
for(i=1;i<n;i++)
 I1 =I1+F[i];
I1 =I1+ ((F[0] + F[n])/2);
I1 =I1*h;
//Симпсона
nechet = 0;
for(i=1; i<=n; i=i+2)
nechet  = nechet+F[i]*4;
I2 = nechet;
chet = 0;
for(i=2; i<=n-1; i=i+2)
 chet = chet+F[i];
I2 =I2+ (chet*2 + F[0] + F[n-1]);
I2 =I2* h/3;
//Средних прямоугольников
I3=0;
for(i=0,i<=n-1; i++;);
   F[i]=((F[i]+F[i+1])/2);
I3+=F[i];
 for (x=a; x<=b; x+=h);
I3+=(x*atan(x));
I3=I3*h;
printf("\n integral po formule trapeciy raven = %f \n",I1);
printf("\n integral po formule Simpsona raven = %f \n",I2);
printf("\n integral po formule srednih prymougolnikov raven = %f \n",I3);
getch();
return 0;
}