Moved plotting and advancing the simulation into simulate() function and changed initial charge distribution