code: def drdpgr(body, lon, lat, alt, re, f):
      
      body = stypes.stringToCharP(body)
      lon = ctypes.c_double(lon)
      lat = ctypes.c_double(lat)
      alt = ctypes.c_double(alt)
      re = ctypes.c_double(re)
      f = ctypes.c_double(f)
      jacobi = stypes.emptyDoubleMatrix()
      libspice.drdpgr_c(body, lon, lat, alt, re, f, jacobi)
      return stypes.cMatrixToNumpy(jacobi)
docs: This routine computes the Jacobian matrix of the transformation
      from planetographic to rectangular coordinates.
  
      http://naif.jpl.nasa.gov/pub/naif/toolkit_docs/C/cspice/drdpgr_c.html
  
      :param body: Body with which coordinate system is associated.
      :type body: str
      :param lon: Planetographic longitude of a point (radians).
      :type lon: float
      :param lat: Planetographic latitude of a point (radians).
      :type lat: float
      :param alt: Altitude of a point above reference spheroid.
      :type alt: 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