Content deleted Content added
Line 21:
where <math>\delta_p^2</math> is the central difference operator for the ''p''-coordinate.
After performing a stability analysis, it can be shown that this method will be stable for any r because this method is implicit.
But problem is the final matrix to solve is not tridiagonal, but is banded. Depending on the ordering (the numbering of the points), the bandwidth of matrix can be very large, and this causes the solution procedure to be more difficult.
| |||