The numerical simulation of non-linear electron-positron oscillations is reported, showing the evolution of the electric field and the plasma number density for large amplitude disturbances. Sharp density gradients and changes in the oscillation frequency are demonstrated, and a new analytical framework is presented to illustrate these phenomena, particularly in the context of pulsar plasmas.