MHB Bound of Euler method- nonuniform partition

evinda
Gold Member
MHB
Messages
3,741
Reaction score
0
Hello! (Wave)

Consider a nonuniform partition $a=t_0< t_1< \dots < t_{\nu}=b$ and assume that if $h_n=t^{n+1}-t^n, 0 \leq n \leq N-1 $ is the changeable step, then $\min_{n} h_n > \lambda \max_{n} h_n, \lambda>0$ independent of $n$.
Show a bound of the error of Euler method analogous to

$$||y^n-y(t^n)|| \leq h c_1 \frac{e^{L(t^n-a)}}{2L} \max_{t \in [a,t^n]} ||y^{(2)} (t)||$$

where $h=\max_{n} h_n$.

According to my lecture notes:We will use the following lemma:

Let $\delta$ be a positive number and $K, d_0, d_1, \dots$ non-negative numbers such that

$$d_{i+1} \leq (1+ \delta) d_i+K, i=0,1, \dots$$

Then it holds

$$d_n \leq e^{n \delta} d_0+ K \frac{e^{n \delta}-1}{\delta}, n=0,1,2, \dots$$

$$\left\{\begin{matrix}
y^{n+1}=y^n+h_n f(t^n,y^n) &, n=0,1, \dots, N-1 \\
y^0=y(0) &
\end{matrix}\right.$$

Expanding $y$ with Taylor, we have:

$$y(t^{n+1})=y(t^n)+ h_n y'(t^n)+ \frac{h^2}{2} y''(\xi_n), \xi_n \in (t^n, t^{n+1})$$

From the differential equation we have:

$$y(t^{n+1})=y(t^n)+ h_n f(t^n, y(t^n))+ \frac{h^2}{2} y''(\xi_n), \xi_n \in (t^n, t^{n+1})$$

$$y(t^{n+1})-y^{n+1}=y(t^n)-y^n+h_n [f(t^n, y(t^n))-f(t^n,y^n)]+ \frac{h h_n}{2} y''( \xi_n)$$

$$\epsilon^{n+1}= \epsilon^n+ h_n [f(t^n, y(t^n))-f(t^n,y^n)]+ \frac{h h_n}{2} y''(\xi_n)$$

and so on...

Why is it as follows?

$$y(t^{n+1})=y(t^n)+ h_n y'(t^n)+ \frac{h^2}{2} y''(\xi_n), \xi_n \in (t^n, t^{n+1})$$Shouldn't we have $h_n$ at each point where we have the step, i.e. shouldn't it be as follows?

$$y(t^{n+1})=y(t^n)+ h_n y'(t^n)+ \frac{h_n^2}{2} y''(\xi_n), \xi_n \in (t^n, t^{n+1})$$

Or am I wrong? (Worried)
 
Mathematics news on Phys.org
Hey! (Mmm)

evinda said:
Shouldn't we have $h_n$ at each point where we have the step, i.e. shouldn't it be as follows?

$$y(t^{n+1})=y(t^n)+ h_n y'(t^n)+ \frac{h_n^2}{2} y''(\xi_n), \xi_n \in (t^n, t^{n+1})$$

Or am I wrong? (Worried)

You're right. Looks like a bit of sloppiness. (Nod)
 
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...
Fermat's Last Theorem has long been one of the most famous mathematical problems, and is now one of the most famous theorems. It simply states that the equation $$ a^n+b^n=c^n $$ has no solutions with positive integers if ##n>2.## It was named after Pierre de Fermat (1607-1665). The problem itself stems from the book Arithmetica by Diophantus of Alexandria. It gained popularity because Fermat noted in his copy "Cubum autem in duos cubos, aut quadratoquadratum in duos quadratoquadratos, et...
I'm interested to know whether the equation $$1 = 2 - \frac{1}{2 - \frac{1}{2 - \cdots}}$$ is true or not. It can be shown easily that if the continued fraction converges, it cannot converge to anything else than 1. It seems that if the continued fraction converges, the convergence is very slow. The apparent slowness of the convergence makes it difficult to estimate the presence of true convergence numerically. At the moment I don't know whether this converges or not.

Similar threads

Replies
2
Views
1K
Replies
1
Views
1K
Replies
38
Views
11K
Replies
1
Views
2K
Replies
3
Views
1K
Replies
3
Views
2K
Back
Top