Abstract
In this study, the Taylor-Galerkin method is modified to take into consideration the third order term in the Taylor series of the fundamental variable. In the Taylor-Galerkin method, after expressing the governing equation of motion in conservation form, the temporal discretization is done first and then spatial discretization follows in contrast to the conventional approaches. A predictor-corrector type algorithm has been developed previously by the same author. A new computationally efficient direct algorithm is proposed in this study. A study on convergency and accuracy of the solution is carried out. Numerical examples show that this new algorithm exhibits the same order of accuracy with less computational effort.