Convert your polar coordinates to Cartesian using pol2cart, and then you can plot everything in Cartesian coordinates. Check the below pseudo code.
theta = 0:0.01:2*pi;
rho = sin(theta).*cos(theta);
% convert polar coordinates to cartesian
[x,y] = pol2cart(theta,rho) ;
N = 10 ;
cx = (5-1)*rand(1,N)+5 ;
cy = (5-1)*rand(1,N)+5 ;
figure
hold on
for i = 1:N
plot(cx(i),cy(i),'sb')
plot(cx(i)+x,cy(i)+y,'r')
end
