c=3e8;lam=2e-6; k=2*pi/lam; n=[1 1+0.25i]; t=3e-15*linspace(0,30,1000); x=1e-6*linspace(0,5,200); y=1e-6*linspace(0,5,200); w=c*abs(k); [xx,yy]= meshgrid(x,y); z=zeros(size(xx)); for m=1:length(t) z=exp(1i*k*(n(1)*xx+n(2)*yy)-1i*w*t(m)); surf(xx,yy,real(z)); % hold on % surf(xx,yy,imag(z)); % hold off % axis([0 5e-6 0 5e-6 -1 1]); shading interp; colormap gray; pause(0.01); end;