C/Visual C - Métodos numéricos en C

#include <iostream>

#include <cmath>

// Definir la función de la que se quiere calcular el área

double f(double x) {

    return x * x; // Ejemplo: f(x) = x^2

}

// Implementar el método del trapecio

double trapecio(double (*func)(double), double a, double b, int n) {

    double h = (b - a) / n; // Ancho de cada subintervalo

    double area = 0.5 * (func(a) + func(b)); // Área inicial con los extremos

    for (int i = 1; i < n; i++) {

        area += func(a + i * h); // Sumar las áreas de los trapecios

    }

    area *= h; // Multiplicar por el ancho del subintervalo

    return area;

}

int main() {

    double a = 0; // Límite inferior

    double b = 1; // Límite superior

    int n = 10;   // Número de subintervalos

    double area = trapecio(f, a, b, n);

    std::cout << "El área bajo la curva es: " << area << std::endl;

    return 0;

}