The Kurganov and Levy scheme [SIAM J.Sci.Comput.22, p1467,2000] which is a semi-discrete numerical solution of hyperbolic systems is applied to the system of one dimensional electrostatic fluid equations and modified to model perturbations and shockwaves in an electrostatic plasma. For an initial density perturbations in the system, we illustrate how the method captures the formation and evolut...