La integración numérica, mediante los métodos trapezoidal y de Simpson, exigen poco esfuerzo de programación en C++. Aquí los usamos para encontrar el área, entre 0 y 6, para la parábola y = x² (aunque no necesitamos la forma explícita de la función; sólo los valores de y como función de x). Los valores de la función se establecen para valores de x igualmente espaciados (h fijo) y para un número par de sub-intervalos h (número impar de puntos). Estos valores se disponen en el archivo datos.in (carpeta datos); tal como se señala a continuación:

7
1

0
1
4
9
16
25
36

y se llaman desde el programa principal. Los dos primeros valores corresponden al número total de datos de f(x) y h; respectivamente. El programa (*.zip) se puede bajar de aquí:

http://blogel.espaciolinux.info/wp-content/archivos/zeka/INTEGRACION.zip

Después de desempaquetar, en cónsola nos movemos hacia la carpeta INTEGRACION y allí compilamos con:

g++ integracion.c++ -o integracion

La ejecución, con ./integracion, produciría una salida similar a esta:

El área calculada por el método de Simpson (función en simpson.h) es más exacta que la determinada por la regla trapezoidal (función en trapezoide.h).