Breaking down a numerical integration into two parts

2 views (last 30 days)
Sorry for this stupic question.
Suppose we have this integral
Of couse this can be solved analytically. But I still want to solve it numerically as:
dx = 0.01;
x = 1:dx:4;
fx = x.^2 + 5*x;
I = trapz(x,fx)
I = 58.5001
Now, we break down this integral into two part. Analytically, we write it as:
I am confused about the lower limit of the second integral when solving it numerically. Should it be 2 or . If we use 2, will not it mean that we are counting one point two times, once in and once in ?

Accepted Answer

Torsten
Torsten on 6 Mar 2024
I1 = dx/2 * f(x0) + dx * sum_{i=1}^{n-1} f(xi) + dx/2 * f(xn)
I2 = dx/2 * f(xn) + dx * sum_{i=n+1}^{n+m-1} f(xi) + dx/2 * f(x_{n+m})
->
I1 + I2 = dx/2 * f(x0) + dx * sum_{i=1}^{n+m-1} f(xi) + dx/2 * f(x_{n+m})
Thus since the endpoints in the trapezoidal rule are only counted half, all is fine.

More Answers (0)

Categories

Find more on Numerical Integration and Differential Equations in Help Center and File Exchange

Products


Release

R2023b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!