function [theta, phi] = conversion(s) phi = atan2(s(1,:), s(2,:)); theta = acos(s(3,:)); end