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