I have an absorption curve that looks like this: http://inspirehep.net/record/1230861/files/MgII_all_stacked.png The x-axis is wavelength and the y-axis is intensity. Both my x and y-axes are lists containing a huge number of data. Example: Wavelength = [1,2,3,4,......] --> x-axis Intensity = [6,4,8,3,......] --> y-axis

I am trying to overplot a gaussian curve on top of it. This is my Gaussian equation: exp(-x^(2)/2sigma^(2)). How do I overplot this Gaussian into my absorption curve? I've tried several things already, but it keeps giving me a straight line. So sorry I'm not very good in Python. Would really appreciate any help that you can provide.