code: def rquad(a, b, c):
      
      a = ctypes.c_double(a)
      b = ctypes.c_double(b)
      c = ctypes.c_double(c)
      root1 = stypes.emptyDoubleVector(2)
      root2 = stypes.emptyDoubleVector(2)
      libspice.rquad_c(a, b, c, root1, root2)
      return stypes.cVectorToPython(root1), stypes.cVectorToPython(root2)
docs: Find the roots of a quadratic equation.
  
      http://naif.jpl.nasa.gov/pub/naif/toolkit_docs/C/cspice/rquad_c.html
  
      :param a: Coefficient of quadratic term.
      :type a: float
      :param b: Coefficient of linear term.
      :type b: float
      :param c: Constant.
      :type c: float
      :return: Root built from positive and negative discriminant term.
      :rtype: tuple