Abstract: | A linearized three‐level difference scheme on nonuniform meshes is derived by the method of the reduction of order for the Dirichlet boundary value problem of the nonlinear parabolic systems. It is proved that the difference scheme is uniquely solvable and second order convergent in L∞‐norm. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 638–652, 2003 |