Simpson 3/8

 #include <iostream>

#include <math.h>

using namespace std;

float f(float x)

{

    return 1/(1+x*x);

}

int main()

{

    int i,n;

    float x0,xn,h,sum;

    cout<<"Enter no ofintervals: ";

    cin>>n;

    cout<<"Enter x0 & x1:";

    cin>>x0>>xn;

    h = (xn-x0)/n;

    cout<<h<<"\n";

    sum = f(x0) + f(xn);

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

    {

        if(i%3 == 0)

        {

            sum+=2*f(x0+i*h);

        }

        else

        {

            sum+=3*f(x0+i * h);

        }

    }

    cout<<"Integral value = "<<(3*(h/8))*sum;

    return 0;

}