a = 0.1 ;
b = 0.2 ;
N = 10 ;
x = (b-a).*rand(N,1) + a;
y = (b-a).*rand(N,1) + a;
mx = mean(x) ; my = mean(y) ; % main point/ target
plot(x,y,'.r') ;
hold on
plot(mx,my,'*b') ;
% path
tx = linspace(0,2*pi,100) ;
ty = sin(tx) ;
for i = 1:100
xi = x+tx(i) ; yi = y+ty(i) ;
mxi = mx+tx(i) ; myi = my+ty(i) ;
plot(xi,yi,'.r') ;
hold on
plot(mxi,myi,'*b') ;
plot(tx,ty,'b') ;
axis([min(tx) max(tx) min(ty) max(ty)])
drawnow
pause(0.1)
hold off
end
