a = 2;
t = linspace(0, 2*pi, 1000);

x = a*cos(3*t).*cos(t);
y = a*cos(3*t).*sin(t);

plot(x,y)
title('Three-Petal Rose')
xlabel('x')
ylabel('y')
axis equal