//Trapezoidal Method for the evaluation of Definite Integrals
#include<iostream>
#include<cmath>
using namespace std;
double f(double x) //write the function whose definite integral is to be calcuated here
{
double a=1/(1+x*x);
return a;
}
int main()
{
int n,i; //n is for subintervals and i is for loop
double a,b,h,sum=0,integral;
cout<<"Enter the limits of integration,\nInitial limit,a="; //get the limits of integration
cin>>a;
cout<<"Final limit, b=";
cin>>b;
cout<<"Enter the no. of subintervals, n="; //get the no. of subintervals
cin>>n;
double x[n+1],y[n+1];
h=(b-a)/n; //get the width of the subintervals
for (i=0;i<=n;i++)
{ //loop to evaluate x0,...xn and y0,...yn
x[i]=a+i*h; //and store them in arrays
y[i]=f(x[i]);
}
for (i=1;i<n;i++) //loop to evaluate h*(y1+...+yn-1)
{
sum=sum+h*y[i];
}
integral=h/2.0*(y[0]+y[n])+sum; //h/2*[y0+yn+2(y1+y2+y3+...yn-1)]
cout<<"The definite integral is "<<integral<<endl;
return 0;
}
Output:
Video Explanation of the code:
I’m a physicist specializing in computational material science with a PhD in Physics from Friedrich-Schiller University Jena, Germany. I write efficient codes for simulating light-matter interactions at atomic scales. I like to develop Physics, DFT, and Machine Learning related apps and software from time to time. Can code in most of the popular languages. I like to share my knowledge in Physics and applications using this Blog and a YouTube channel.
[wpedon id="7041" align="center"]
Hello, I want to find area under curve, but with indefinite integration, can you share to me the code using C++ ??
Thank you.
Hi Rizky!
Area under the curve always implies definite integration.
In the above code we achieve this by using a Numerical Method called Trapezoidal Method.
Similarly there are other Numerical methods to calculate the definite integrals, like SImpson Rules, Gauss Quadrature, etc.
All these methods are Numerical.
If you want to calculate the area under the curve or some definite integral in the Symbolic(Analytical) way, then it is very hard to using C++ and not very useful.
You can perform symbolic integration using Mathematica.
how to find centroid using trapezoidal rule and how to transfer it into the c++ code ? I got problem to solve this problem. Maybe you or anyone could help me. Thank you.
It will not work
do you have a C++ program for Romberg integration with trapezoidal rule??