1; %EM field definition Ed=0.1; %drive Edc=conj(Ed); Ep=1.0*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 ];