An O(k<sup>2</sup>+kh<sup>2</sup>+h<sup>2</sup>) Accurate Two-level Implicit Cubic Spline Method for One Space Dimensional Quasi-linear Parabolic Equations

In this piece of work, using three spatial grid points, we discuss a new two-level implicit cubic spline method of O(k2 + kh2 + h4) for the solution of quasi-linear parabolic equation , 0< x <1, t > 0 subject to appropriate initial and Dirichlet boundary conditions, where h > 0, k > 0 are grid sizes in space and time-directions, respectively. The cubic spline approximation produces at each time level a spline function which may be used to obtain the solution at any point in the range of the space variable. The proposed cubic spline method is applicable to parabolic equations having singularity. The stability analysis for diffusion- convection equation shows the unconditionally stable character of the cubic spline method. The numerical tests are performed and comparative results are provided to illustrate the usefulness of the proposed method.