% rotation(a,n) displays n points of the orbit of 0 under rotation by a

function [void] = rotation(a,n)

cla
t=0:0.01:1;
%plot(cos(2*pi*t),sin(2*pi*t));
hold on

for i=0:n,
  plot(cos(2*pi*i*a),sin(2*pi*i*a),'.','markersize',12);
end

axis off

text(0,0,['\alpha = ' num2str(a,8)],'fontsize',20,'horizontalalignment','center');