1; %EM field definition Ep=1.0; Epc=conj(Ep); Ed=10; Edc=conj(Ed); wd=w13; wp=w12; modulation_freq=[0, wp, wd, -wp, -wd, wp-wd, wd-wp]; E_field =[0, Ep, Ed, Epc, Edc, 0, 0 ]; %modulation_freq=[0, wp, -wp]; %E_field =[0, Ep, Epc];