u0=4*pi*1e-7;

e0=1e-9/(36*pi);

Z0=(u0/e0)^0.5;

f=1e8;

w=2*pi*f;

k=w*(u0*e0)^0.5;

phi_E=0;

phi_H=pi;

EE=20;

HH=EE/Z0;

z=0:0.1:20;

m0=zeros(size(z));

gifname='mag_motion.gif';

figure

for t=0:1:100    

  Ex=EE*cos(k*z-w*t*1e-9+phi_E);    

  Hx=HH*cos(k*z-w*t*1e-9+phi_H);        

  plot3(z,m0,Ex,'b','LineWidth',2);    

  hold on    

  plot3(z,Hx,m0,'r','LineWidth',2); 

  hold off    

  xlabel('way')    

  ylabel('B-Hx')    

  zlabel('E-Ex')    

  title(['the work','t=',num2str(t),'ns'],'fontsize',14)    

  set(gca,'fontsize',12)    

  drawnow    

  frame=getframe(1);    

  im=frame2im(frame);    

  [imind,cm]=rgb2ind(im,500);    

  if t==0;        

       imwrite(imind,cm,gifname,'gif');    

  else        

       imwrite(imind,cm,gifname,'gif','WriteMode','append','DelayTime',0.1);    

  end

end