polar plot r=1+sin(theta)