Take the 2-minute tour ×
Mathematics Stack Exchange is a question and answer site for people studying math at any level and professionals in related fields. It's 100% free, no registration required.

In my Numerical Methods for PDEs textbook by Ari Uscher L-stability and stiff decay are introduced by considering a generalized test equation:

$y' = \lambda (y - g(t)), 0 < t < b$

where g(t) is a bounded but otherwise arbitrary function. As $Re(\lambda) \rightarrow -\infty$ the exact solution satisfies $y(t) \rightarrow g(t), 0< t < b$. My first question (1) is that this is not obvious to me and no explanation is given. How can one conclude that?

It goes on to define stiff decar for a method as $|y_n - g(t_n)| \rightarrow 0$ as $\Delta t Re(\lambda) \rightarrow 0$. My second question (2) is how this is found in practice. If I want to determine if an ODE is stiff, how do I find/choose a g(t) to satisfy the above requirement? Or does any arbitrary bounded g satisfy it if my ODE is stiff as the test equation supposedly shows (which I don't understand and is listed as question 1).

share|improve this question

Your Answer

 
discard

By posting your answer, you agree to the privacy policy and terms of service.

Browse other questions tagged or ask your own question.