It seems to me that it would be easier for you to work in spherical coordinates. In the thread, Deriving formulae for spheres, I posted several equations for deriving properties of spheres. The equation for deriving a vector from the center of a sphere located at the origin to the surface of the sphere is
[math]\text{V}\left(\rho,\,\varphi,\,\theta\right)=\left\langle\,\rho\,\text{sin}\,\varphi\,\text{cos}\,\theta,\,\,\rho\,\text{sin}\,\varphi\,\text{sin}\,\theta,\,\,\rho\,\text{cos}\,\varphi\,\right\rangle[/math]
where [math]\rho[/math] is the radius [math]\left\{\rho\in\mathbb{R}\,|\,0\le\rho<\infty\right\}[/math], [math]\varphi[/math] is the inclination [math]\left\{\varphi\in\mathbb{R}\,|\,0\le\varphi\le \pi\right\}[/math], and [math]\theta[/math] is the azimuth [math]\left\{\theta\in\mathbb{R}\,|\,0\le\theta\le 2\pi\right\}[/math].
You can use the above equation to derive the vectors you need.