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