yyP = (-bb+np.sqrt(bb*bb - 4*x[2]*cc))/(2*x[2])