sorry, I didn’t understand.

Is this code true:

closeall;

clear;

w0 = 2*pi*c/1482e-9; # center frequency

Imin = 10; Imax = 10;

Ey = getdata(‘monitor’,‘Ey’);

Esy = getdata(‘monitor_7’,‘Ey’);

t = getdata(‘monitor’,‘t’);

I = linspace(Imin, Imax, length(t));

zero_pad = 2^16;

w = fftw(t,1,zero_pad);

filter1 = 2*exp(-(w-w0)^2/(150e12)^2);

Ex = getdata(‘monitor’,‘Ex’);

Esx = getdata(‘monitor_7’,‘Ex’);

t = getdata(‘monitor’,‘t’);

I = linspace(Imin, Imax, length(t));

zero_pad = 2^16;

w = fftw(t,1,zero_pad);

filter1 = 2*exp(-(w-w0)^2/(150e12)^2);

# remove high frequency components y

Ey_through_w = 2*((1:length(w))<=(length(w)/2+0.1))*fft(pinch(Ey),1,zero_pad);

Ey_through_t = invfft(pinch(Ey_through_w)*filter1);*

Ey_through_t = Ey_through_t(1:length(t));

Esy_through_w = 2((1:length(w))<=(length(w)/2+0.1))*fft(pinch(Esy),1,zero_pad);

Esy_through_t = invfft(pinch(Esy_through_w)*filter1);

Esy_through_t = Esy_through_t(1:length(t));

# remove high frequency components x

Ex_through_w = 2*((1:length(w))<=(length(w)/2+0.1))*fft(pinch(Ey),1,zero_pad);

Ex_through_t = invfft(pinch(Ey_through_w)*filter1);*

Ex_through_t = Ey_through_t(1:length(t));

Esx_through_w = 2((1:length(w))<=(length(w)/2+0.1))*fft(pinch(Esy),1,zero_pad);

Esx_through_t = invfft(pinch(Esy_through_w)*filter1);

Esx_through_t = Esy_through_t(1:length(t));

# plot results

plot(I, sqrt(abs(Ex)^2+abs(Ey)^2)/sqrt(abs(Esx)^2+abs(Esy)^2),‘I (MW/cm^2)’,’(E_out/E_source)^2’);

legend(“Increasing Source Intensity”);