  In this work, we suggest an easy-to-code higher-order finite volume
semi-discrete scheme to analyze the nonlinear behavior of the electron-plasma
oscillations by solving electron fluid equations numerically. The present
method employs a fourth-order accurate centrally weighted essentially
nonoscillatory reconstruction (CWENO4) polynomial for estimating the numerical
flux at the grid-cell interfaces, and a fourth-order Runge-Kutta method for the
time integration. The numerical implementation is validated by reproducing
earlier results for both non-dissipative and dissipative cold plasmas. The
stability of the present scheme is illustrated by evolving the nonlinear
electron plasma oscillations in a cold non-dissipative plasma for hundred
plasma periods, which also display a negligible numerical dissipation. The
fourth-order accuracy of the existing approach is also confirmed by evaluating
the convergence of errors for the nonlinear electron plasma oscillations in a
cold non-dissipative plasma.
