Fixed bug in the Q matrix calculation and the univariate factorization function.

[ginac.git] / check / exam_factor.cpp

diff --git a/check/exam_factor.cpp b/check/exam_factor.cpp
index ed064588e11bbabb42ff14829a946bfa42f89ec5..bf1751a572c02e7bf1cb13e0cad5d12514bbb78c 100644 (file)
--- a/check/exam_factor.cpp
+++ b/check/exam_factor.cpp
@@ -58,6 +58,9 @@ static unsigned exam_factor1()
        e = ex("(1+x)^3", syms);
        result += check_factor(e);
 
        e = ex("(1+x)^3", syms);
        result += check_factor(e);
 

+       e = ex("(x+1)*(x+4)", syms);
+       result += check_factor(e);
+

        e = ex("x^6-3*x^5+x^4-3*x^3-x^2-3*x+1", syms);
        result += check_factor(e);
 
        e = ex("x^6-3*x^5+x^4-3*x^3-x^2-3*x+1", syms);
        result += check_factor(e);
 


@@ -79,7 +82,11 @@ static unsigned exam_factor1()
        e = ex("x^8-40*x^6+352*x^4-960*x^2+576", syms);
        result += check_factor(e);
 
        e = ex("x^8-40*x^6+352*x^4-960*x^2+576", syms);
        result += check_factor(e);
 

-       e = ex("x*(2+x^2)*(1+x+x^3+x^2+x^6+x^5+x^4)*(1+x^3)^2*(-1+x)", syms);
+       e = ex("x*(2+x^2)*(1+x+x^3+x^2+x^6+x^5+x^4)*(1+x)^2*(1-x+x^2)^2*(-1+x)", syms);
+
+       result += check_factor(e);
+
+       e = ex("(x+4+x^2-x^3+43*x^4)*(x+1-x^2-3*x^3+4*x^4)", syms);

        result += check_factor(e);
 
        return result;
        result += check_factor(e);
 
        return result;