docs: def drdgeo(lon, lat, alt, re, f):
      """
      This routine computes the Jacobian of the transformation from
      geodetic to rectangular coordinates.
  
      http://naif.jpl.nasa.gov/pub/naif/toolkit_docs/C/cspice/drdgeo_c.html
  
      :param lon: Geodetic longitude of point (radians).
      :type lon: float
      :param lat: Geodetic latitude of point (radians).
      :type lat: float
      :param alt: Altitude of point above the 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
      """
      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.drdgeo_c(lon, lat, alt, re, f, jacobi)
      return stypes.cMatrixToNumpy(jacobi)