MHB Asymptotic error formula for the trapezoidal rule

  • Thread starter Thread starter kalish1
  • Start date Start date
  • Tags Tags
    Error Formula
Click For Summary
SUMMARY

The discussion centers on using the asymptotic error formula for the trapezoidal rule to estimate the number of subdivisions required to evaluate the integral $\int_{0}^{2}e^{-x^2}dx$ with an accuracy of $\epsilon=10^{-10}$. The user calculated the error estimate as $E_n^T(f)\approx \frac{h^2}{4e^4}$ and derived that $n \approx 6767$ is insufficient. A more accurate calculation led to $N \ge 15628$ for the trapezoidal scheme, ensuring the desired accuracy. The approximate integral is computed using the trapezoidal sum with these subdivisions.

PREREQUISITES
  • Understanding of the trapezoidal rule in numerical integration
  • Familiarity with asymptotic error analysis
  • Knowledge of calculus, specifically differentiation and integration
  • Basic programming skills for implementing numerical methods
NEXT STEPS
  • Study the derivation of the trapezoidal rule error formula
  • Learn about numerical integration techniques in Python using libraries like NumPy
  • Explore the implications of error bounds in numerical methods
  • Investigate alternative numerical integration methods such as Simpson's rule
USEFUL FOR

Mathematicians, engineers, and computer scientists involved in numerical analysis, particularly those focused on improving the accuracy of integral approximations.

kalish1
Messages
79
Reaction score
0
I need to use the asymptotic error formula for the trapezoidal rule to estimate the number n of subdivisions to evaluate $\int_{0}^{2}e^{-x^2}dx$ to the accuracy $\epsilon=10^{-10}$. I also need to find the approximate integral in this case. I would like to know if my attempt is correct. Thanks in advance for any help.

**My attempt:** $E_n^T(f)\approx -h^2/12[f'(b)-f'(a)]. f(x)=e^{-x^2}, f'(x)=-2xe^{-x^2}.$

So $E_n^T(f)\approx h^2/(4e^4)$ and since $h=1/n$, we have to find an n that satisfies the inequality $\frac{1}{4e^4n^2}\leq 10^{-10}$. We obtain $n \approx 6767$. The approximate integral is ??
 
Mathematics news on Phys.org
Hmm. I don't quite get what you get. I have
\begin{align*}
E_{N}(f)& \approx - \frac{(b-a)^{2}}{12 N^{2}}[f'(b)-f'(a)] \\
&= \frac{2^{2}}{12 N^{2}}[4e^{-4}-0] \\
&= \frac{4}{3 N^{2}e^{4}}.
\end{align*}
Hence, we need
$$ \frac{4}{3 N^{2}e^{4}} \le 10^{-10},$$
or
$$N \ge \frac{2 \cdot 10^{5}}{ \sqrt{3} \, e^{2}} \approx 15627.2,$$
so let $N=15628$. You must set up a trapezoidal scheme with this many sub-intervals, and evaluate. You're going to have:
$$A \approx \frac{2-0}{2 \cdot 15628} \sum_{j=0}^{15627}[f(x_{j})+f(x_{j+1})]
=\frac{1}{15628} \sum_{j=0}^{15627}[e^{-(2j/15628)^{2}}+e^{-(2(j+1)/15628)^{2}}]$$
$$=\frac{1}{15628} \sum_{j=0}^{15627}[e^{-(j/7814)^{2}}+e^{-((j+1)/7814)^{2}}].$$
 
Great! CCC. (Crystal clear clarification.)
 

Similar threads

Replies
6
Views
2K
  • · Replies 6 ·
Replies
6
Views
2K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 18 ·
Replies
18
Views
3K
  • · Replies 3 ·
Replies
3
Views
2K
Replies
6
Views
2K
  • · Replies 2 ·
Replies
2
Views
1K
  • · Replies 38 ·
2
Replies
38
Views
11K
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 4 ·
Replies
4
Views
1K