code: def spherical(coordinates):
      
      c = coordinates
      r = N.linalg.norm(c,axis=0)
      theta = N.arccos(c[2]/r)
      phi = N.arctan2(c[1],c[0])
      return N.column_stack((r,theta,phi))
docs: No error is propagated