xt1=@(t) 16*(sin(t)).^3;
yt1=@(t) 13*cos(t)-5*cos(2*t)-cos(4*t);
fplot(xt1,yt1,[-6,8])
hold on
xt2=@(t) t-sin(t);
yt2=@(t) 1-cos(t);
fplot(xt2,yt2,[0,7])
hold off