1; %EM field definition Ed=10.12; %drive Edc=conj(Ed); Ep=0.8*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); %light_positive_freq = [wp, wd, wp-wd]; E_field = [Ep, Ed, 0 ];