[octave]График функции

Надо построить график функции F: R^3->R^3 (в 3d, короче) в octave.

Функция задана параметрически. Параметры:

r = sqrt(x^2 + y^2 + z^2)
theta = atan2( sqrt(x^2+y^2), z )
phi = atan2(y,x)

Соответственно:

x = f (r,theta,phi)
y = g (r,theta, phi)
z = h (r, theta, phi)

пробовал так:
octave:1> x = [-3:0.1:3]
octave:2> z=y=x
octave:3> [xx,yy,zz] = meshgrid (x,y,z)
warning: broken pipe — some output may be lost
octave:4> r = sqrt (xx.^2 + yy.^2+zz.^2)
warning: broken pipe — some output may be lost
octave:5> theta = atan2(sqrt (xx.^2+yy.^2),zz)
error: invalid conversion of NDArray to Matrix
error: evaluating assignment expression near line 5, column 7

В общем sqrt (..) - это NDArray, а zz - матрица. Как их сделать одного типа?? И вообще в том ли направлении я двигаюсь?