plot(r%3D1-sin(theta))