

tau = linspace(0, 2*pi, 101);
s = linspace(0, 4*pi, 201);

 %M = moviein(41);
     plot(s -sin(s), 1-cos(s), 'g')
     hold on
for j = 0:120
    t = j*pi/30;
    
     h1 = plot(t+cos(tau),1+sin(tau))
     h2 = plot(t-sin(t), 1-cos(t), '*r')
     pause
     axis([-1, 4*pi+1,-.5,6])
     delete(h1)
     delete(h2)
      %M(:,j+1) = getframe
 end
 
