trapozoid_integration.cpp

# include <stdio.h>
# include <math.h>

double f(double x)
{
    return ( exp(x*x) );
}
int main()
{
    double l1=0;
    double l2=2;

    int i=0;
    double h = (l1+l2)/10;
    double x=l1;
    double y=0;
    double sum=0;
    while(x <= l2)
    {
     y = f(x);
        if (i==0 || i== 9)
        {
            sum +=y;
        }
        else
        {
            sum+=2*y;        
        }
        i+=1;
        x+=h;
    }
    sum =(double)h/2*sum;
    printf("%lf",sum);
}