| 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 |
|
} |
|
|
