1; %EM field definition Ep=0.01; %probe Epc=conj(Ep); Ed=.1; %drive Edc=conj(Ed); 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];