N=160;
XY=zeros（1N）;
wave1=zeros（1N）;?wave1（1:N*3/4）=sin（（N*3/4:-1:1）/N*5*pi）;
wave2=zeros（1N）;
wave3=zeros（1N）;
for?loop1=N*3/4+1:4*N
????wave2=[wave2（2:N）?wave1（N）];
????wave1（2:N）=wave1（1:N-1）;
????wave1（1）=sin（loop1/N*5*pi）;
????E_Total=10*（wave2-wave1）;
????H_Total=10*（wave2+wave1）;
%?E&H
????figure（1）;?quiver3（XYXY