Uno de los metodos para el calculo numerico de la integral de una funcion f(x) en el intervalo a < b, es el llamado metodo de los rectangulos :
[a] : se contruye una funcion de retorna los valores de f(x).
: se divide el intervalo a-b en partes iguales, por ejemplo dx = (b-a)/10.
[c] : nuestra intencion es calcular la suma de las areas de los rectangulos
f(x)*dx, para obtener el valor del area total bajo la curva de f(x) entre a y b.
[d] : hacemos float sum = 0 y calculamos :
for(float x = a;x < b;x += dx)
sum += f(x);
sum *= dx; nos da el valor parcial del area buscada;
[c] : hacemos dx \= 10 y repetimos el proceso hasta que
[1] : llegamos a un valor predefinido de dx;
[2] : o la diferencia entre dos valores sucesivos de sum es menor que la aproximacion pedida.
Espero haber sido claro, es mas complicado escribir la explicacion que mandar el programa, pero en este foro no permiten hacerlo.
Por cualquier duda, estoy a tu disposicion.