MHB Asymptotic error formula for the trapezoidal rule

  • Thread starter Thread starter kalish1
  • Start date Start date
  • Tags Tags
    Error Formula
AI Thread Summary
The discussion focuses on using the asymptotic error formula for the trapezoidal rule to estimate the number of subdivisions needed to evaluate the integral of e^{-x^2} from 0 to 2 with a specified accuracy of ε=10^{-10}. The initial attempt calculates the error and suggests n ≈ 6767, but a subsequent correction indicates that n should actually be around 15628 to meet the accuracy requirement. The trapezoidal scheme is then set up with this number of sub-intervals to approximate the integral. The final expression for the integral involves a summation of function values at calculated points, demonstrating the application of the trapezoidal rule. The discussion concludes with a clear method for achieving the desired accuracy in the integral approximation.
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.)
 
Suppose ,instead of the usual x,y coordinate system with an I basis vector along the x -axis and a corresponding j basis vector along the y-axis we instead have a different pair of basis vectors ,call them e and f along their respective axes. I have seen that this is an important subject in maths My question is what physical applications does such a model apply to? I am asking here because I have devoted quite a lot of time in the past to understanding convectors and the dual...
Insights auto threads is broken atm, so I'm manually creating these for new Insight articles. In Dirac’s Principles of Quantum Mechanics published in 1930 he introduced a “convenient notation” he referred to as a “delta function” which he treated as a continuum analog to the discrete Kronecker delta. The Kronecker delta is simply the indexed components of the identity operator in matrix algebra Source: https://www.physicsforums.com/insights/what-exactly-is-diracs-delta-function/ by...
Back
Top