Did you mean: plot (r =1 -sin (zeta ))?
plot (r =1-sin (theta))
plot(r=1 - sin(theta))
Want to compute something more complicated? Try a full Python/SymPy console at SymPy Live.