clc; clear;close;
% x=0:0.01:2*pi;
% y=zeros(size(x));
% 
% y1=(x>pi/4 & x<pi/2).*1;%(x-pi/4);
% y2=(x>=pi/2&x<3*pi/4).*1;%(-x+3*pi/4);
% y=y1+y2;
% 
% 
% id=plot(x,y);
% %axis([0 2*pi -1 1])
% yy=y;
% while(1)% ii=1:500
% yy(1:end-1)=y(2:end);
% yy(end)=y(1);
% y=yy;
% set(id,'YData',y)
% pause(0.01)
% end

[x,y]=meshgrid(0:0.01:1,-1:0.01:1);
r=4*pi*sqrt((x+1).^2+y.^2);
z1=besselj(0,r)+i*bessely(0,r);%exp(i*4*pi*x);
z=real(z1);
id=surface(x,y,z)
for ii=1:1000;
z=real(z1.*exp(-i*0.01*ii));
set(id,'ZData',z)
pause(0.01)

end