plot (r =1 -sin (zeta ))