I confess I haven't run the code, but I did read it, and I noticed this:
def original_electric_field(t, y):
#return plane_wave(t, y, angular_frequency)
return gaussian(t, y)
This describes a short voltage pulse, centred on frequency zero. However, what you actually want to model is a pulse of light, centred on a non-zero frequency. I think you need to multiply the Gaussian by the plane wave.
I haven't tried it, sorry. Maybe I have misunderstood.
I confess I haven't run the code, but I did read it, and I noticed this:
This describes a short voltage pulse, centred on frequency zero. However, what you actually want to model is a pulse of light, centred on a non-zero frequency. I think you need to multiply the Gaussian by the plane wave.
I haven't tried it, sorry. Maybe I have misunderstood.