<SYSTEM_TASK:>
Solve the following problem using Python, implementing the functions described below, one line at a time
<END_TASK>
<USER_TASK:>
Description:
 def rquad(a, b, c):
 """ 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 """

    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)