[shiqiang]

theta=linspace(0,2*pi,64);
x=sin(theta);
y=cos(theta);
plot(x,y);
axis equal;
hold on;
for i=1:64;
a=[0;x(i)];
b=[0;y(i)];
line(a,b);
end