Integración numérica: trapecios, código C | Numerical integration: trapezoids, 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();
float a,b,n,dx,suma,xi,AT;
int i;
gotoxy(10,6); printf("            INTEGRACION NUMERICA TRAPEZOIDAL ");
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",&a);
gotoxy(10,11);printf ("Digite el punto superior ");
scanf ("%f",&b);
gotoxy(10,12);printf ("Digite el numero de sub-areas ");
scanf ("%f",&n);
dx=(b-a)/n;
suma=0;
i=1;
while (i<=n-1)
{
xi=a+i*dx;
suma=suma+F(xi);
i=i+1;
}
AT=(b-a)*((F(a)+2*suma+F(b))/(2*n));
gotoxy(10,14);printf ("El area es: %f \n",AT);
gotoxy(10,15);printf("Presione una tecla para continuar");
getch();
}