Content deleted Content added
wikilinks |
replace r by \Delta t since r is not defined in the article |
||
Line 19:
\left(u_{ij}^{n+1}+u_{ij}^n\right)</math>
where <math>\delta_p</math> is the central difference operator for the ''p''-coordinate. After performing a [[Von Neumann stability analysis|stability analysis]], it can be shown that this method will be stable for any <math>\Delta t</math>.
A disadvantage of the Crank–Nicolson method is that the matrix in the above equation is [[band matrix|banded]] with a band width that is generally quite large. This makes the solution of the equation quite costly.
|