code: def drdsph(r, colat, lon):
      
      r = ctypes.c_double(r)
      colat = ctypes.c_double(colat)
      lon = ctypes.c_double(lon)
      jacobi = stypes.emptyDoubleMatrix()
      libspice.drdsph_c(r, colat, lon, jacobi)
      return stypes.cMatrixToNumpy(jacobi)
docs: This routine computes the Jacobian of the transformation from
      spherical to rectangular coordinates.
  
      http://naif.jpl.nasa.gov/pub/naif/toolkit_docs/C/cspice/drdsph_c.html
  
      :param r: Distance of a point from the origin.
      :type r: float
      :param colat: Angle of the point from the positive z-axis.
      :type colat: float
      :param lon: Angle of the point from the xy plane.
      :type lon: float
      :return: Matrix of partial derivatives.
      :rtype: 3x3-Element Array of floats