how to plot measurements data

I have a data set containing photometric measurement from galaxy, you can think of this like intensity. now, let's say I know that the data is coming from normal distribution and I know it's mean. how can I plot the curve that will fit to the distribution based on the data. or to be precise, I know that plotting the measurements suppose to look like this

when every color is different data, but when I'm trying to plot I get this is what I get

I upload the data file, Ignore I'd and z coulombs. the code I try is

x = linspace(4858-50/2,4858+50,numel(gp));
plot(x,gp)

[Birdman]

Can you share your data?

[Jan]

And the code you use for plotting?

[avner yakov]

I update the question