clc; clear all; close all;
f=0:0.001:2*pi;
rad=5
x=rad*cos(f);
y=rad*sin(f);
plot(x,y,'o-')
x1=[0 0];
y1=[0 4];
%hold on
axis equal
rad1=sqrt(x1(2)^2+y1(2)^2);
for i=0:10
alfa=2*pi/60*(-i)+pi/2;
x1(2)=rad1*cos(alfa);
y1(2)=rad1*sin(alfa);
plot(x,y,'o-');
hold on
plot(x1,y1)
axis equal
hold off
pause(1)
end