How Does the Implicit Euler Method Affect the Amplification Factor for y' = λy?

  • Context: MHB 
  • Thread starter Thread starter ra_forever8
  • Start date Start date
  • Tags Tags
    advanced Odes Topic
Click For Summary
SUMMARY

The discussion focuses on the application of the implicit Euler method to the differential equation y' = λy, where λ is negative. The modified equation derived from the implicit Euler scheme is y_(n+1) = y_n / (1 + hλ), indicating that the amplification factor is influenced by the stability condition |1 - hλ| ≥ 1. When λ is negative, the solution exhibits decay, confirming stability. The analysis concludes that for λ = ±1, the stability of the solution can be further examined.

PREREQUISITES
  • Understanding of differential equations, specifically first-order linear equations.
  • Familiarity with numerical methods, particularly the implicit Euler method.
  • Knowledge of stability analysis in numerical solutions.
  • Basic proficiency in logarithmic and exponential functions.
NEXT STEPS
  • Study the implications of the implicit Euler method on stability for various values of λ.
  • Explore the derivation and application of the modified equation for different numerical methods.
  • Investigate the behavior of solutions for nonlinear differential equations using numerical methods.
  • Learn about the convergence criteria for numerical methods applied to differential equations.
USEFUL FOR

Mathematicians, engineers, and students in applied mathematics or numerical analysis who are interested in the stability of numerical methods for solving differential equations.

ra_forever8
Messages
106
Reaction score
0
Fint the modified equation when the implicit Euler method is applied to y'= f(y). If f(y)=λy, where λ is negative. what is the effect on the amplication factor?
=>

y ' = lambda * y
dy / dx = lambda * y
dy / y = lambda dx
ln y = lambda * x + C
y = [ C / lambda ]e^( lambda * x )

This is what I have done so far. Can anyone help me if Iam going right direction to answer the question?

 
Physics news on Phys.org
grandy said:
Fint the modified equation when the implicit Euler method is applied to y'= f(y). If f(y)=λy, where λ is negative. what is the effect on the amplication factor?
=>

y ' = lambda * y
dy / dx = lambda * y
dy / y = lambda dx
ln y = lambda * x + C
y = [ C / lambda ]e^( lambda * x )

This is what I have done so far. Can anyone help me if Iam going right direction to answer the question?



I agree with the solution up to $\displaystyle \begin{align*} \ln{|y|} = \lambda x + C \end{align*}$, from there you should get $\displaystyle \begin{align*} y = A e^{\lambda x } \end{align*}$, the constant factor does not depend on $\displaystyle \begin{align*} \lambda \end{align*}$...
 
y ' = λ * y
dy / dx = λ * y
dy / y = λ dx
ln y = λ* x + C
y = Ae^( λ* x ), the constant factor does not depend on λ.
i SOLVE THIS FOR THE ACTUAL SOLUTIONS
NOW,
The implicit Euler scheme is given by:
y_(n+1)= y_n +hf(y_n+1 , t_n+1 )
For f(y)=λ y, we have:
y_(n+1)= y_n +hf(y_n+1 , t_n+1 )= y_n + h λ y_n+1
Solving this for y_n+1 (in general, this is not possible), we arrive at:
y_n+1 = y_n / (1-h λ).....(eqn 1)

From (eqn 1), we can see that if |1−h λ|≥1, the solution is decaying (stable). Compare this to the actual solution of y(x)=Ae^(λ x).
If we have λ being negative, we would have:
y_n+1 = y_n / (1+h λ).....(eqn 2)

Compare this to the actual solution of y(x)=Ae^(−λ x). What conclusion can i draw?
Trying it for λ= ±1 , what happens to stability.
 

Similar threads

Replies
3
Views
2K
Replies
6
Views
2K
  • · Replies 2 ·
Replies
2
Views
5K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 1 ·
Replies
1
Views
8K
  • · Replies 4 ·
Replies
4
Views
2K
Replies
3
Views
3K
  • · Replies 8 ·
Replies
8
Views
7K
Replies
1
Views
4K