y=u.*sin(u*cos(g).^2.*(a + p*cot(m) - 1.0*u*cos(g).*cos(v)) + cos(g).*(p*v + u*sin(g)).*(p*sin(v) ...