x=-2*pi:0.05:2*pi;
for n=-2*pi:2*pi;
subplot (2,2,1), plot(n*cos(x).*(1+cos(x)), n*sin(x).*(1+cos(x)))
subplot(2,2,2), plot(n*(cos(x)).^3,n*(sin(x)).^3)
subplot(2,2,3), plot(n*(x-sin(x)),n*(1-cos(x)))
subplot(2,2,4), plot(2*n*cos(x)-n*cos(2*x),2*n*sin(x)+n*sin(2*x))
end