Before revision:  A novel approach to the error estimation is also proposed . Numerical examples are given to illustrate the effectiveness and practicability, with emphasis on the potential applications in plasma physics.
After revision:  A novel feature of this algorithm is the ability to estimate the error level automatically . Numerical examples are given to illustrate the effectiveness and practicability, with emphasis on the potential applications in plasma physics.