**Source Code:**

//integration of given function using Simpson's 3/8 rule#include<stdio.h>float y(float x){return 1/(1+x*x); //function of which integration is to be calculated}int main(){float x0,xn,h,s;int i,n,j,flag;printf("Enter x0, xn, no. of subintervals: ");scanf("%f%f%d",&x0,&xn,&n);h = (xn-x0)/n;s = y(x0)+y(xn);for(i = 1; i<=n-1;i++){for(j=1;j<=n-1;j++){if(i==3*j){flag = 1;break;}elseflag = 0;}if(flag==0)s += 3*y(x0+i*h);elses += 2*y(x0+i*h);}printf("Value of integral is %6.4f\n",(3*h/8)*s);return 0;}

## No comments:

## Post a Comment