clear;
rc=1;rm=1/3;
t=0:360;
xo=(rc+rm)*cosd(t)-rm*cosd(t*(rc+rm)/rm);
yo=(rc+rm)*sind(t)-rm*sind(t*(rc+rm)/rm);
xi=(rc-rm)*cosd(t)+rm*cosd(t*(rc-rm)/rm);
yi=(rc-rm)*sind(t)-rm*sind(t*(rc-rm)/rm);
plot(xo,yo,'b',xi,yi,'r');
grid('on');
title('cycloid curve');
xlabel('x');ylabel('y');
axis('square');
axis([-2,2,-2,2]);
