code: def dpgrdr(body, x, y, z, re, f):
      
      body = stypes.stringToCharP(body)
      x = ctypes.c_double(x)
      y = ctypes.c_double(y)
      z = ctypes.c_double(z)
      re = ctypes.c_double(re)
      f = ctypes.c_double(f)
      jacobi = stypes.emptyDoubleMatrix()
      libspice.dpgrdr_c(body, x, y, z, re, f, jacobi)
      return stypes.cMatrixToNumpy(jacobi)
docs: This routine computes the Jacobian matrix of the transformation
      from rectangular to planetographic coordinates.
  
      http://naif.jpl.nasa.gov/pub/naif/toolkit_docs/C/cspice/dpgrdr_c.html
  
      :param body: Body with which coordinate system is associated.
      :type body: str
      :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
      :param re: Equatorial radius of the reference spheroid.
      :type re: float
      :param f: Flattening coefficient.
      :type f: float
      :return: Matrix of partial derivatives.
      :rtype: 3x3-Element Array of floats