Integración numérica: Simpson 1/3 (un tercio), código C | Numerical integration: Simpson 1/3, code C

#include <stdio.h>
#include <math.h>
#include <conio.h>
#include <stdlib.h>
#define F(x)(x*x)-4
#define FD(x)(2*x)
#define FP(x)(x*x)+x-4
#define FE(x,y)(x*y)+1

void main()
{
clrscr();
double li,ls,suma=0,h,x1,x2,AT,s1,s2;
int i=1,n;
gotoxy(10,6); printf("            INTEGRACION NUMERICA SIMPSON 1/3 ");
gotoxy(10,7);printf("------------------------------------------------------------");
gotoxy(10,8);printf(" Este codigo fue hecho por Jair Beltran USCO Garzon ");
gotoxy(10,9);printf("Se hallara el area de la funcion f(x)=x^2-4");
gotoxy(10,10);printf ("Digite el punto inferior ");
scanf ("%lf",&li);
gotoxy(10,11);printf ("Digite el punto superior ");
scanf ("%lf",&ls);
gotoxy(10,12);printf ("Digite el numero de sub-areas ");
scanf ("%d",&n);
h=(ls-li)/(2.0*n);
while (i<=n)
{ x1=li+((2*i-1)*h);
x2=li+(2*i*h);
s1=F(x1);s2=F(x2);
suma=suma+2.0*s1+s2;
i=i+1;
}
s1=F(li);s2=F(ls);
AT=(h/3.0)*(s1-s2+2.0*suma);
gotoxy(10,14);printf ("El area es: %f \n",AT);
gotoxy(10,15);printf("Presione una tecla para continuar");
getch();
}