@spiceErrorCheck<EOL>def dpgrdr(body, x, y, z, re, f):
body = stypes.stringToCharP(body)<EOL>x = ctypes.c_double(x)<EOL>y = ctypes.c_double(y)<EOL>z = ctypes.c_double(z)<EOL>re = ctypes.c_double(re)<EOL>f = ctypes.c_double(f)<EOL>jacobi = stypes.emptyDoubleMatrix()<EOL>libspice.dpgrdr_c(body, x, y, z, re, f, jacobi)<EOL>return stypes.cMatrixToNumpy(jacobi)<EOL>
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