write a program to integrate a function using Simpson 3/8 rule in c++.

write a program to integrate a function using Simpson 3/8 rule in c++.

CODE:

#include<iostream.h>

#include<conio.h>

float f(float x){

return(1/(1+x));

}

void main(){

clrscr();

float a,b,h,sum1=0,sum2=0,sum=0;

int n,i;

cout<<"enter lower limit : ";

cin>>a;

cout<<"enter upper limit : ";

cin>>b;

cout<<"enter number of intervals : ";

cin>>n;

h=(b-a)/n;

sum=f(a)+f(b);

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

if((i%3)!=0){

sum1+=f(a+i*h);

}}

sum1=3*sum1;

for(i=2;i<=n-3;i++){

if((i%3)==0){

sum2+=f(a+i*h);

}

}

sum2=2*sum2;

sum=sum+sum1+sum2;

sum=(sum*3*h)/8;

cout<<"x\ty"<<endl;

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

cout<<i<<"/6\t"<<f(a+i*h)<<endl;

}

cout<<"\n value of integral using simpson's 3\8 rule: "<<sum;

getch();

}

OUTPUT: