From: Andrej V. <an...@us...> - 2009-01-07 09:56:20
|
Update of /cvsroot/maxima/maxima/share/contrib/solve_rec In directory 23jxhf1.ch3.sourceforge.com:/tmp/cvs-serv26765 Modified Files: rtest_simplify_sum.mac simplify_sum.mac simplify_sum_test.mac Log Message: - to_hgfred: factor the expression before computing shift quotient - added new tests Index: rtest_simplify_sum.mac =================================================================== RCS file: /cvsroot/maxima/maxima/share/contrib/solve_rec/rtest_simplify_sum.mac,v retrieving revision 1.14 retrieving revision 1.15 diff -u -d -r1.14 -r1.15 --- rtest_simplify_sum.mac 6 Jan 2009 10:27:39 -0000 1.14 +++ rtest_simplify_sum.mac 7 Jan 2009 09:56:15 -0000 1.15 @@ -448,3 +448,21 @@ -(17*sqrt(2)*log(3-2*sqrt(2))-24*log(2)-68)/8 ); 0; + +test_sum( + sum(((x-1)^(2*n+1)*(x+1)^(-2*n-1))/(2*n+1),n,0,inf), + [x<1, x>0], + false, + [rootscontract, lambda([u], scanmap(factor, u)), ratsimp], + log(x)/2 +); +0; + +test_sum( + sum(-1/(2*k)+1/(2*2*k+1)+1/(2*2*k+3), k, 1, inf), + [], + [ratsimp], + (9*log(2)-8)/6 +); +0; + \ No newline at end of file Index: simplify_sum.mac =================================================================== RCS file: /cvsroot/maxima/maxima/share/contrib/solve_rec/simplify_sum.mac,v retrieving revision 1.22 retrieving revision 1.23 diff -u -d -r1.22 -r1.23 --- simplify_sum.mac 6 Jan 2009 10:27:39 -0000 1.22 +++ simplify_sum.mac 7 Jan 2009 09:56:15 -0000 1.23 @@ -565,7 +565,7 @@ while subst(var=lo, expr)=0 do lo : lo+1, expr:subst(var=var+lo, expr), - quo : ratsimp(shiftQuo(makefact(expr)*var!, var)), + quo : ratsimp(shiftQuo(factor(makefact(expr)*var!), var)), if not(?ratp(quo, var)) then return(false), ss_print_message(2, "Shift quotient", quo), Index: simplify_sum_test.mac =================================================================== RCS file: /cvsroot/maxima/maxima/share/contrib/solve_rec/simplify_sum_test.mac,v retrieving revision 1.9 retrieving revision 1.10 diff -u -d -r1.9 -r1.10 --- simplify_sum_test.mac 6 Jan 2009 10:27:39 -0000 1.9 +++ simplify_sum_test.mac 7 Jan 2009 09:56:15 -0000 1.10 @@ -375,3 +375,18 @@ false, [] ); + + +test_sum( + sum(((x-1)^(2*n+1)*(x+1)^(-2*n-1))/(2*n+1),n,0,inf), + [x<1, x>0], + false, + [rootscontract, lambda([u], scanmap(factor, u)), ratsimp] +); + +test_sum( + sum(-1/(2*k)+1/(2*2*k+1)+1/(2*2*k+3), k, 1, inf), + [], + false, + [ratsimp] +); |