Tuesday 20 September 2011

Simpson's 1.3 Rule



/* This Program finds area by Simpson's 1/3 rule*/

#include<stdio.h>
#include<math.h>

float area(float x)
{
float f;
f=log(1+x*x)/(1+x*x);
return(f);
}

void main()
{
int i,n;
float a,b,h,sum1,sum2,sum3,sum;
clrscr();
printf("Type the lower limit a and upper limit b\n");
scanf("%f  %f",&a,&b);
printf("Type number of interval n\n");
scanf("%d",&n);
sum1=area(a);
sum1+=area(b);
sum2=0;
sum3=0;
h=(b-a)/n;
for(i=1;i<n;i+=2)
sum2+=area(a+i*h);
for(i=2;i<n-1;i+=2)
sum3+=area(a+i*h);
sum=h*(sum1+4*sum2+2*sum3)/3;
printf("The required area = %8.4f\n",sum);
getch();
}

No comments:

Post a Comment