After committing, I did change a few things that seemed to be bugs. With my
version
(not yet on CVS), I get
(%i1) load("altsimp.lisp")$
(%i2) nusum(n*n!,n,0,n);
Dependent equations eliminated: (1)
(%o2) (n+1)!-1
(%i3) nusum(n^4*4^n/binomial(2*n,n),n,0,n);
(%o3) (2*(n+1)*(63*n^4+112*n^3+18*n^2-22*n+3)*4^n)/(693*binomial
(2*n,n))-2/(33*7)
(%i4) unsum(%,n);
(%o4) (n^4*4^n)/binomial(2*n,n)
(%i5) unsum(product(i^2,i,1,n),n);
(%o5) (product(i^2,i,1,n-1))*(n-1)*(n+1)
I'll do some more tests; if everything checks out, I'll make these changes.
I don't want altsimp to be a distraction from fixing bugs & extending
Maxima
is useful ways. Altsimp is pretty much a 100% re-write. Over the years, a
great
deal of effort has gone into simp.lisp and simplus. A from scratch re-write
might
be a bad idea. A more conservative method for making Maxima infinity
correct might
be better.
Barton
-----maxima-bounces at math.utexas.edu wrote: -----
>I?have?downloaded?altsimp.lisp?and?tried?to?execute?the?code.
>
>In?a?Maxima?session?I?loaded?the?code?with?load("altsimp.lisp").?But?with
>gcl
>2.6.8?and?CLISP?2.46?on?my?system?I?get?fatal?errors?and?the?Maxima
>session?is
>killed.
>
>I?have?the?same?problems?when?I?compile?with?compile_file("altsimp.lisp")
>and
>load?the?object-file.
>
>With?CLISP?I?get?a?"Program?stack?overflow"?in?rtest9.mac.
>
>When?I?use?a?batch-mode?I?get?the?following?error?with?GCL?2.6.8?in
>rexamples.mac?(at?this?point?the?testsuite?seems?to?kill?my?Maxima
>session):
>
>**********************?Problem?79?***************
>Input:
>(kill(all),?example(nusum))
>
>
>Result:
>(%i1)?nusum(n*n!,n,0,n)
>
>Dependent?equations?eliminated:??(1)
>(%o1)????????????????????????????(n?+?1)!?-?1
>(%i2)?nusum(n^4*4^n/binomial(2*n,n),n,0,n)
>???????????????????????????4????????3???????2??????????????n
>????????????2?(n?+?1)?(63?n??+?112?n??+?18?n??-?22?n?+?3)?4?????2
>(%o2)???????------------------------------------------------?-?----
>??????????????????????????693?binomial(2?n,?n)?????????????????33?7
>(%i3)?unsum(%,n)
>?????????????????????????????????????4??n
>????????????????????????????????????n??4
>(%o3)??????????????????????????----------------
>???????????????????????????????binomial(2?n,?n)
>(%i4)?unsum(product(i^2,i,1,n),n)
>error-catch
>
>This?differed?from?the?expected?result:
>Done
>
>Do?you?have?an?idea?what?I?should?do?