Numerical instability
 | This article or section is in a state of significant expansion or restructuring. You are welcome to assist in its construction by editing it as well. If this article or section has not been edited in several days, please remove this template. If you are the editor who added this template and you are actively editing, please be sure to replace this template with {{in use}} during the active editing session. Click on the link for template parameters to use.
This redirect was last edited by Alapidus (talk | contribs) 19 years ago. (Update timer) |
Numerical Instability deals with an unexpected phenomenon, whose provenence is very obscure, which arises from the numerical solution of differential equations. The usual notion of instability includes concepts such as fickleness, unreliability, and oversensitivity. If a computation is numerically unstable it gives inaccurate results. A major source of numerical instability results from a cumulative magnification effect of small errors on the numerical results.
There are open questions at the frontiers of numerical analysis which try to deal with issues of numerical instability, such as nonlinear instability. This article is written to be accessible to those who need only understand the simplest definition of a differential equation.
A typical computation will be analyzed algebraically to reveal how the phenomenon arises. A special case of the solution of the equation
y'(x)=f(x,y(x)) with y(0)=1 on the closed interval [0,1]
will be examined. This equation is studied very closely in the subject of Ordinary Differential Equations, but we'll solve it simply using plausible notions found in elementary calculus and then evaluate the numerical stability of the method.
For the simplified (but useful) approach this grid is used:
GN={x | xn = n*h n=0,1,...,N with h=1/N}. N is a positive integer.
Other simplifications can be made as we go along, in order to elucidate the main point. For example, consider approximate values of y at xn. These approximations of y(xn) will be denoted by yn.