code: def dsphdr(x, y, z):
      
      x = ctypes.c_double(x)
      y = ctypes.c_double(y)
      z = ctypes.c_double(z)
      jacobi = stypes.emptyDoubleMatrix()
      libspice.dsphdr_c(x, y, z, jacobi)
      return stypes.cMatrixToNumpy(jacobi)
docs: This routine computes the Jacobian of the transformation from
      rectangular to spherical coordinates.
  
      http://naif.jpl.nasa.gov/pub/naif/toolkit_docs/C/cspice/dsphdr_c.html
  
  
      :param x: X-coordinate of point.
      :type x: float
      :param y: Y-coordinate of point.
      :type y: float
      :param z: Z-coordinate of point.
      :type z: float
      :return: Matrix of partial derivatives.
      :rtype: 3x3-Element Array of floats