On 6/7/07, George Leeman <george.leeman at sbcglobal.net> wrote:
> eibe*%e^(%i*phi)*rho*s+%e^(2*%i*phi+2*%i*%pi*m+%i*al)*r*rho^2+2*%e^(%i*%pi*m
> )*r*rho+%e^(2*%i*phi+2*%i*%pi*m+%i*al)*r
>
> (%i2) solve( %, eibe );
>
> Maxima encountered a Lisp error: Error in PROGN [or a callee]: Bind stack
> overflow.
Ugh. That's pretty bad.
It appears that solve (or some callee) calls RATDISREP over and over.
1 Enter ?ratdisrep
[eibe*%e^(%i*phi)*rho*s+%e^(2*%i*phi+2*%i*%pi*m+%i*al)*r*rho^2
+2*%e^(%i*%pi*m)*r*rho
+%e^(2*%i*phi+2*%i*%pi*m+%i*al)*r]
2 Enter ?ratdisrep [%i]
2 Exit ?ratdisrep %i
2 Enter ?ratdisrep [phi]
2 Exit ?ratdisrep phi
2 Enter ?ratdisrep [%i]
2 Exit ?ratdisrep %i
2 Enter ?ratdisrep [2*phi+2*%pi*m+al]
2 Exit ?ratdisrep 2*phi+2*%pi*m+al
2 Enter ?ratdisrep [%i]
2 Exit ?ratdisrep %i
2 Enter ?ratdisrep [2*phi+2*%pi*m+al-%pi]
2 Exit ?ratdisrep 2*phi+2*%pi*m+al-%pi
2 Enter ?ratdisrep [%i]
2 Exit ?ratdisrep %i
2 Enter ?ratdisrep [2*phi+2*%pi*m+al-2*%pi]
2 Exit ?ratdisrep 2*phi+2*%pi*m+al-2*%pi
2 Enter ?ratdisrep [%i]
2 Exit ?ratdisrep %i
2 Enter ?ratdisrep [2*phi+2*%pi*m+al-3*%pi]
2 Exit ?ratdisrep 2*phi+2*%pi*m+al-3*%pi
2 Enter ?ratdisrep [%i]
2 Exit ?ratdisrep %i
2 Enter ?ratdisrep [2*phi+2*%pi*m+al-4*%pi]
2 Exit ?ratdisrep 2*phi+2*%pi*m+al-4*%pi
2 Enter ?ratdisrep [%i]
2 Exit ?ratdisrep %i
2 Enter ?ratdisrep [2*phi+2*%pi*m+al-5*%pi]
2 Exit ?ratdisrep 2*phi+2*%pi*m+al-5*%pi
Not sure what's up with subtracting %pi over and over.
The functions TMS and TIMESIN (not sure what they do)
might also be causing trouble here.
For the record, it looks like linsolve solves this problem successfully.
(I thought solve would call linsolve for linear problems, but that
doesn't appear to be the case.)
HTH
Robert Dodier