/* fichero integra.c */
/* cálculo de la integral definida de una función */
#include<stdio.h>
void main(void)
{
        int n, k;
        double area=0.0, h;
        double f(double);
        printf("\nIntroduzca el número de trapecios de la aproximación: ");
        scanf("%d", &n);
        h = 3./n;
        for (k=0; k<n; k++)
                area += ((f(k*h)+f((k+1)*h))*h)/2;
        printf("\nEl área calculada es %9.3lf", area);
}
double f(double x){
        double y;
        y = x*x*x+x*x-5*x+3;
        return y;
}