You can't really compare the two pictures as they are for completely different structures. The classical picture is for a classical structure, i.e. a bulk p-n junction (in Si).
0.9 is a very high charge on the electrodes, if you work out the equivalent doping concentration. You need to work backwards from a reasonable doping level and convert that to a number of electrons.
Actually you didn't attach the script for plotting, but rather for generating the data, but perhaps there is a wrong unit conversion or you are plotting grid indices instead of Z coordinates.
PS: Please don't ask the question "is my result correct". It doesn't have an answer.