[1], [2], [3]

                
                        delta = b*b - 4*a*c;
a1 = (-b + raiz(delta)) / (2 * a);
a2 = (-b - raiz(delta)) / (2 * a);