1; %EM field definition Ed=0.1; %drive Edc=conj(Ed); Ep=0.06*Ed; %probe Epc=conj(Ep); Em=-Ep; % opposite sideband (resulting from EOM modulation of drive) Emc=conj(Em); %wd=w13; %wp=w12; %wm=wd-(wp-wd); %modulation_freq=[0, wp, wd, wm, -wp, -wd, -wm, wp-wd, wd-wp]; E_field =[0, Ep, Ed, Em, Epc, Edc, Emc, 0, 0 ]; %modulation_freq=[0, wp, -wp]; %E_field =[0, Ep, Epc];