Integración numérica: rectángulos, código C++ | Numerical integration: rectangles, code C++

#include <iostream>
#include <cstdio>
#include <cmath>
using namespace std;
float fun (float x);
int main() {
float s,a,b,n,dx,i,xi,AT;
printf ("Este codigo fue hecho por Jair Beltran USCO Garzon");
printf ("Digite el punto inferior ");
scanf ("%f",&a);
printf ("Digite el punto superior ");
scanf ("%f",&b);
printf ("Digite el numero de sub-areas ");
scanf ("%f",&n);
dx=(b-a)/n;
s=0;
i=1;
while (i<=n)
{
xi=a+(2*i-1)*(dx/2);
s=s+fun(xi);
i=i+1;
}
AT=dx*s;
printf ("El area es:%f \n",AT);
return 0;
}
float fun(float x)
{float f;
f=(x*x)-4;
return f;
}