function [x_theta,x_phi,r]=cart_2_sph(nodi)

% come nella funzione cart2sph di matlab, ma con la seguente differenza:
% phi (angolo sull'asse z) sta in [0,pi] (e non in [-pi/2,pi/2])
% theta (angolo sul piano x,y) sta in [0,2*pi] (e non in [-pi,pi])

[x_theta,x_phi,r]=cart2sph(nodi(:,1),nodi(:,2),nodi(:,3));

for i=1:1:size(x_theta)
    if x_theta(i) < 0
        x_theta(i)=2*pi+x_theta(i);
    end
end
for i=1:1:size(x_phi)
    if x_phi(i)>=0
        x_phi(i)=-(x_phi(i)-pi/2);
    end
    if x_phi(i)<0
        x_phi(i)=(-x_phi(i)+pi/2);
    end
end