Solution (source code)

= Solution
{c}

The <ellipsoid> is the <level set> $F(x,y,z)=2x^2+y^2+z^2=1$. Its <gradient> is $(4x,2y,2z)$, which never vanishes on the <ellipsoid>. Half this <gradient> gives the outward <normal vector>
$$
\boxed{\mathbf n(\theta,\phi)=(\sqrt2\cos\theta,\sin\theta\cos\phi,\sin\theta\sin\phi).}
$$
This <normal vector> remains valid at $\theta=0,\pi$, where the given <parametrization> is singular.