def dlatdr(x, y, z):
    """
    This routine computes the Jacobian of the transformation from
    rectangular to latitudinal coordinates.

    http://naif.jpl.nasa.gov/pub/naif/toolkit_docs/C/cspice/dlatdr_c.html

    :param x: X-coordinate of point.
    :type x: float
    :param y: Y-coordinate of point.
    :type y: float
    :param z: Z-coord
    :type z: float
    :return: Matrix of partial derivatives.
    :rtype: 3x3-Element Array of floats
    """
    x = ctypes.c_double(x)
    y = ctypes.c_double(y)
    z = ctypes.c_double(z)
    jacobi = stypes.emptyDoubleMatrix()
    libspice.dlatdr_c(x, y, z, jacobi)
    return stypes.cMatrixToNumpy(jacobi)
This routine computes the Jacobian of the transformation from
    rectangular to latitudinal coordinates.

    http://naif.jpl.nasa.gov/pub/naif/toolkit_docs/C/cspice/dlatdr_c.html

    :param x: X-coordinate of point.
    :type x: float
    :param y: Y-coordinate of point.
    :type y: float
    :param z: Z-coord
    :type z: float
    :return: Matrix of partial derivatives.
    :rtype: 3x3-Element Array of floats