Integración numérica: Simpson 3/8 (tres octavos), código C | Numerical integration: Simpson 3/8, 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();
clrscr();
float s1,s2,s3,li,ls,n,dx,suma,x1,x2,x3,h,AT;
int i;
gotoxy(10,6); printf("            INTEGRACION NUMERICA SIMPSON 3/8 ");
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 ("%f",&li);
gotoxy(10,11);printf ("Digite el punto inferior ");
scanf ("%f",&ls);
gotoxy(10,12);printf ("Digite el numero de sub-areas ");
scanf ("%f",&n);
dx=(ls-li)/(3*n);
suma=0;
i=1;
while (i<=n)
{
x1=li+((3*i-2)*dx);
x2=li+((3*i-1)*dx);
x3=li+(3*i*dx);
s1=F(x1);s2=F(x2);s3=F(x3);
suma=suma+3*s1+3*s2+2*s3;
i=i+1;
}
s1=F(li);s2=F(ls);
AT=3*(dx/8)*(s1-s2+suma);
gotoxy(10,14);printf ("El area es: %f \n",AT);
gotoxy(10,15);printf("Presione una tecla para continuar");
getch();
}