anglev=3000*2*pi/60.0 %angular velocity rad/sec
dt=2*pi/anglev/600; %T/600
t=(0:dt:60); %sec %one minute total
v=25; %mm/sec
xo=0+v*t; %center of rotation
yo=xo*0;
R=10 ; %10mm away from xo, yo
angle=0+t*anglev; %angle relative to x axis
xc=xo+R*cos(angle); %center of cirlcle
yc=yo+R*sin(mod(angle,2*pi));
%draw the circle
r=5; %5mm circle
angc=(0:0.1:2*pi);
for i=1:length(t)
x=xc(i)+cos(angc)*r; %circle
y=yc(i)+sin(angc)*r;
x(length(angc)+1)=x(1); %close the circle
y(length(angc)+1)=y(1);
plot(x,y,'r-'); %draw circle
axis([xo(i)-20,xo(i)+20,yo(i)-20,yo(i)+20]);
hold on;
plot([xc(i),xo(i)],[yc(i),yo(i)],'b-') %rotation arm
plot(0,0,'k+'); %origin
xlabel('x(mm)');
ylabel('y(mm)');
title(['time=',num2str(t(i)),' sec']);
pause(0.1)
hold off;
end