#include<stdio.h>
#include<math.h>
double f(double x){
return 1/(5+(3*x));
}
void main(){
int n,i;
double a,b,h,x,sum=0,integral;
printf("\nEnter the no. of sub-intervals: ");
scanf("%d",&n);
printf("\nEnter the initial limit: ");
scanf("%lf",&a);
printf("\nEnter the final limit: ");
scanf("%lf",&b);
h=fabs(b-a)/n;
for(i=1;i<n;i++){
x=a+i*h;
if(i%3==0)
{
sum=sum+2*f(x);
}
else
{
sum=sum+3*f(x);
}
}
integral=(3*(h/8))*(f(a)+f(b)+sum);
printf("\nThe integral is: %lf\n",integral);
}
Comments
Post a Comment