x1 = (-bc+np.sqrt(bc**2-4*a[0]*cc))/(2*a[0])