Error in integral for 5.19.0
- Subject: Error in integral for 5.19.0
- From: Robert Dodier
- Date: Sun, 9 Aug 2009 22:38:06 -0600
On 8/9/09, Richard Hennessy <rich.hennessy at verizon.net> wrote:
> (%i1) integrate(sqrt(2*m*(E[n]-U(x))),x,-x[0],x[0])=(n-1/2)*%pi*hbar;
>
> Maxima encountered a Lisp error:
>
> Error in PROGN [or a callee]: Bind stack overflow.
Looks like a bug all right --- please file a bug report if you haven't already.
Part of a stack trace is shown below.
The calls from SRATSIMP through SIGN-PREP keep repeating.
Why does it jump from SIGN-PREP to SRATSIMP again?
Incidentally (SRATSIMP '((MTIMES SIMP) -1 ((%IMAGPART)
(($X SIMP ARRAY) 0)))) by itself doesn't cause any trouble.
FWIW
Robert Dodier
PS.
53446: (FR-ARGS ((%IMAGPART) (($X SIMP ARRAY) 0)))
53447: (FR1 ((%IMAGPART) (($X SIMP ARRAY) 0)) NIL)
53448: (FULLRATSIMP #<unavailable argument>)
53449: (SRATSIMP ((%IMAGPART) (($X SIMP ARRAY) 0)))
53450: (SIGN-PREP #<unavailable argument>)
53451: (SIGN01 (($X SIMP ARRAY) 0))
53452: ($SIGN (($X SIMP ARRAY) 0))
53453: ($CSIGN (($X SIMP ARRAY) 0))
53454: (SIMP-IMAGPART
((%IMAGPART) (($X SIMP ARRAY) 0))
#<unavailable argument>
NIL)
53455: (FR-ARGS ((%IMAGPART) (($X SIMP ARRAY) 0)))
53456: (FR1 ((%IMAGPART) (($X SIMP ARRAY) 0)) NIL)
53457: (FULLRATSIMP #<unavailable argument>)
53458: (SRATSIMP ((%IMAGPART) (($X SIMP ARRAY) 0)))
53459: (SIGN-PREP #<unavailable argument>)
53460: (SIGN01 (($X SIMP ARRAY) 0))
53461: ($SIGN (($X SIMP ARRAY) 0))
53462: ($CSIGN (($X SIMP ARRAY) 0))
53463: (SIMP-IMAGPART
((%IMAGPART) (($X SIMP ARRAY) 0))
#<unavailable argument>
NIL)
53464: (FR-ARGS ((%IMAGPART) (($X SIMP ARRAY) 0)))
53465: (FR1 ((MTIMES SIMP) -1 ((%IMAGPART) (($X SIMP ARRAY) 0))) NIL)
53466: (FULLRATSIMP #<unavailable argument>)
53467: (SRATSIMP ((MTIMES SIMP) -1 ((%IMAGPART) (($X SIMP ARRAY) 0))))
53468: ($DEFINT
((MTIMES SIMP) ((MEXPT SIMP) 2 ((RAT SIMP) 1 2))
((MEXPT SIMP) ((MTIMES SIMP) $M ((MPLUS SIMP) (# $N) (# -1 #)))
((RAT SIMP) 1 2)))
$X
((MTIMES SIMP) -1 (($X SIMP ARRAY) 0))
(($X SIMP ARRAY) 0))