Simpson 1/3

 #include <iostream>

using namespace std;

int main()

{

    int i,j,k,n;

    float x[10],y[10],sum,h;

    cout<<"Enter no of values in table: ";

    cin>>n;

    cout<<"Enter X corresponding Y\n";

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

    {

        cin>>x[i]>>y[i];

    }

    h = x[1] - x[0];

    sum = y[0]+y[n-1];

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

    {

        if(i%2 == 0)

        {

            sum+=2*y[i];

        }

        else

        {

            sum+=4*y[i];

        }

    }

    cout<<"Integral vlaue = \n"<<(h/3)*sum;

    return 0;

}