theta = 0:0.01:pi;
rho = ((cos(1/2 *pi*cos(theta))-cos(1/2 *pi))./(sin(theta))).^2;
PolarPlot(theta,rho)