C code for SIMPSON'S 1/3 RULE
| 1 | #include<stdio.h> | |||
| 2 | #include<conio.h> | |||
| 3 | #include<math.h> | |||
| 4 | float f(float x) | |||
| 5 | { | |||
| 6 | return(1/(1+x)); | |||
| 7 | } | |||
| 8 | void main() | |||
| 9 | { | |||
| 10 | int i,n; | |||
| 11 | float xn,x0,h,y[20],se,ans,s0,x[20]; | |||
| 12 | printf("\n Enter values of h,xn,x0: "); | |||
| 13 | scanf("%f%f%f",&x0,&h,&xn); | |||
| 14 | n=(xn-x0)/h; | |||
| 15 | if(n%2==1) | |||
| 16 | { | |||
| 17 | n=n+1; | |||
| 18 | } | |||
| 19 | h=(xn-x0)/n; | |||
| 20 | printf("\n Refined_value_of n and h are:%d %f\n",n,h); | |||
| 21 | printf("\n Y values: \n"); | |||
| 22 | for(i=0; i<=n; i++) | |||
| 23 | { | |||
| 24 | x[i]=x0+i*h; | |||
| 25 | y[i]=f(x[i]); | |||
| 26 | printf("\n %f\n",y[i]); | |||
| 27 | } | |||
| 28 | so=0; | |||
| 29 | se=0; | |||
| 30 | for(i=1; i<n; i++) | |||
| 31 | { | |||
| 32 | if(i%2==1) | |||
| 33 | { | |||
| 34 | so=so+y[i]; | |||
| 35 | } | |||
| 36 | else | |||
| 37 | { | |||
| 38 | se=se+y[i]; | |||
| 39 | } | |||
| 40 | } | |||
| 41 | ans=h/3*(y[0]+y[n]+4*so+2*se); | |||
| 42 | printf("\n Final_integration_is %f",ans); | |||
| 43 | getch(); | |||
| 44 | } | |||
