[shiqiang]

t=linspace(0,2*pi,64);
r(1,64)=1;
polar(t,r);
hold on
for i=1:64;
t=t+pi/32;
polar(t,r)
end