グーグルで検索していてここにたどり着きました。
>>405-408
のプログラムを利用させてもらっています。
これに改良を加え座標を方程式で出してその値を極座標に変換できるように
しようと試みたのですがうまくいきません。
次のようにしてみました。

<< Calculus`VectorAnalysis`

Solve[{x + y + z == 1,
　 x - y + z == 3,
x + y - z == 5},
{x, y, z}]

CoordinatesFromCartesian[{x, y, z} /. %, Spherical] {1,
HoldForm[Degree]/Degree, HoldForm[Degree]/Degree}

返数は
CoordinatesFromCartesian[{{4, -1, -2}}, Spherical]
（以下同じような感じで二つ返ってきてしまいます。）
どうもただ式のx, y, zに方程式の解が代入されただけの感じのものです。

どうすれば成功するのでしょうか。アドバイスもらえないでしょうか。