maxima-bugs — Reports from the SF bug-tracking system

[Maxima-bugs] [maxima:bugs] #4601 5.48 regression in spherical_bessel_j integral

[Raymond T. <rt...@us...>]

- **status**: open --> closed
- **Comment**:

Fixed in commit [fc85ad]




---

**[bugs:#4601] 5.48 regression in spherical_bessel_j integral**

**Status:** closed
**Group:** None
**Labels:** integrate 
**Created:** Mon Aug 25, 2025 09:28 AM UTC by Antonio Rojas
**Last Updated:** Tue Aug 26, 2025 02:54 PM UTC
**Owner:** Raymond Toy


In 5.48, the improper integral
~~~
integrate(spherical_bessel_j(1,x)^2,x,0,inf);
~~~
returns 0 instead of the correct answer %pi/6, as it did in 5.47

Bisected to 6ed48036faca8c5fcf996db95647bf09e9d13913


---

Sent from sourceforge.net because max...@li... is subscribed to https://sourceforge.net/p/maxima/bugs/

To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options.  Or, if this is a mailing list, you can unsubscribe from the mailing list.

[Maxima-bugs] [maxima:bugs] #4601 5.48 regression in spherical_bessel_j integral

[Raymond T. <rt...@us...>]

- **labels**:  --> integrate
- **assigned_to**: Raymond Toy



---

**[bugs:#4601] 5.48 regression in spherical_bessel_j integral**

**Status:** open
**Group:** None
**Labels:** integrate 
**Created:** Mon Aug 25, 2025 09:28 AM UTC by Antonio Rojas
**Last Updated:** Tue Aug 26, 2025 02:53 PM UTC
**Owner:** Raymond Toy


In 5.48, the improper integral
~~~
integrate(spherical_bessel_j(1,x)^2,x,0,inf);
~~~
returns 0 instead of the correct answer %pi/6, as it did in 5.47

Bisected to 6ed48036faca8c5fcf996db95647bf09e9d13913


---

Sent from sourceforge.net because max...@li... is subscribed to https://sourceforge.net/p/maxima/bugs/

To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options.  Or, if this is a mailing list, you can unsubscribe from the mailing list.

[Maxima-bugs] [maxima:bugs] Re: #4601 5.48 regression in spherical_bessel_j integral

[Raymond T. <rt...@us...>]

I think I figured out the issue.  A stupid error on my part.  I've uploaded a fix, but want to do a bit more testing before landing.


---

**[bugs:#4601] 5.48 regression in spherical_bessel_j integral**

**Status:** open
**Group:** None
**Created:** Mon Aug 25, 2025 09:28 AM UTC by Antonio Rojas
**Last Updated:** Mon Aug 25, 2025 03:22 PM UTC
**Owner:** nobody


In 5.48, the improper integral
~~~
integrate(spherical_bessel_j(1,x)^2,x,0,inf);
~~~
returns 0 instead of the correct answer %pi/6, as it did in 5.47

Bisected to 6ed48036faca8c5fcf996db95647bf09e9d13913


---

Sent from sourceforge.net because max...@li... is subscribed to https://sourceforge.net/p/maxima/bugs/

To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options.  Or, if this is a mailing list, you can unsubscribe from the mailing list.

[Maxima-bugs] [maxima:bugs] #4601 5.48 regression in spherical_bessel_j integral

[Raymond T. <rt...@us...>]

Thanks for doing the bisection.  Looks like I messed it up.


---

**[bugs:#4601] 5.48 regression in spherical_bessel_j integral**

**Status:** open
**Group:** None
**Created:** Mon Aug 25, 2025 09:28 AM UTC by Antonio Rojas
**Last Updated:** Mon Aug 25, 2025 09:28 AM UTC
**Owner:** nobody


In 5.48, the improper integral
~~~
integrate(spherical_bessel_j(1,x)^2,x,0,inf);
~~~
returns 0 instead of the correct answer %pi/6, as it did in 5.47

Bisected to 6ed48036faca8c5fcf996db95647bf09e9d13913


---

Sent from sourceforge.net because max...@li... is subscribed to https://sourceforge.net/p/maxima/bugs/

To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options.  Or, if this is a mailing list, you can unsubscribe from the mailing list.

[Maxima-bugs] [maxima:bugs] #4602 taylorinfo ignores asymp

[Barton W. <wil...@us...>]

---

**[bugs:#4602] taylorinfo ignores asymp **

**Status:** open
**Group:** None
**Created:** Mon Aug 25, 2025 09:32 AM UTC by Barton Willis
**Last Updated:** Mon Aug 25, 2025 09:32 AM UTC
**Owner:** nobody


I think that `taylorinfo` should note that the expansion is in reciprocal powers
~~~
(%i4) taylor(exp(1/x),[x,0,3,'asymp])$

(%i5) taylorinfo(%);

(%o5) [[x,0,3]]
~~~

Guess:  in lines 3258-3259 in `hayat`, we should change `$asympt` to `$asymp`:

~~~
(setq acc (if (and (fourth qk) (consp (fourth qk)) (eq '$asympt (caar (fourth qk))))
				(list '$asympt) nil))
~~~



---

Sent from sourceforge.net because max...@li... is subscribed to https://sourceforge.net/p/maxima/bugs/

To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options.  Or, if this is a mailing list, you can unsubscribe from the mailing list.

[Maxima-bugs] [maxima:bugs] #4601 5.48 regression in spherical_bessel_j integral

[Antonio R. <aro...@us...>]

---

**[bugs:#4601] 5.48 regression in spherical_bessel_j integral**

**Status:** open
**Group:** None
**Created:** Mon Aug 25, 2025 09:28 AM UTC by Antonio Rojas
**Last Updated:** Mon Aug 25, 2025 09:28 AM UTC
**Owner:** nobody


In 5.48, the improper integral
~~~
integrate(spherical_bessel_j(1,x)^2,x,0,inf);
~~~
returns 0 instead of the correct answer %pi/6, as it did in 5.47

Bisected to 6ed48036faca8c5fcf996db95647bf09e9d13913


---

Sent from sourceforge.net because max...@li... is subscribed to https://sourceforge.net/p/maxima/bugs/

To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options.  Or, if this is a mailing list, you can unsubscribe from the mailing list.

[Maxima-bugs] [maxima:bugs] #4600 Vectorsimp has many undocumented variables controlling its behavior

[Raymond T. <rt...@us...>]

---

**[bugs:#4600] Vectorsimp has many undocumented variables controlling its behavior**

**Status:** open
**Group:** None
**Labels:** Documentation 
**Created:** Sun Aug 24, 2025 08:37 PM UTC by Raymond Toy
**Last Updated:** Sun Aug 24, 2025 08:37 PM UTC
**Owner:** nobody


The documentation for [`vectorsimp`](https://maxima.sourceforge.io/docs/manual/Matrices-and-Linear-Algebra.html#index-vectorsimp) lists many variables that control what `vectorsimp` does.  However, most don't seem to be documented anywhere.

Thees need to be documented.  They do appear to be used in the source code.


---

Sent from sourceforge.net because max...@li... is subscribed to https://sourceforge.net/p/maxima/bugs/

To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options.  Or, if this is a mailing list, you can unsubscribe from the mailing list.

[Maxima-bugs] [maxima:bugs] Re: #4599 intosum not fully recursive

[Stavros M. <mac...@us...>]

Sounds good.


---

**[bugs:#4599] intosum not fully recursive**

**Status:** open
**Group:** None
**Labels:** intosum 
**Created:** Fri Aug 22, 2025 01:25 PM UTC by Barton Willis
**Last Updated:** Fri Aug 22, 2025 04:12 PM UTC
**Owner:** nobody


Here we need to apply `intosum` twice.  The user documentation isn't entirely clear about this, but I think that `intosum` isn't fully recursive.

~~~
	(sumexpand  : true, cauchysum  : true, display2d : false)$
	powerseries(exp(x)/(x-5)^n,x,5);
(%o37) ('sum('sum((5^i10*(x-5)^(i9-i10))/(i10!*(i9-i10)!),i10,0,i9),i9,0,inf))

 /(x-5)^n
	intosum(%);
(%o38) 'sum(('sum((5^i10*(x-5)^(i9-i10))/(i10!*(i9-i10)!),i10,0,i9))/(x-5)^n,

            i9,0,inf)
	intosum(%);
(%o39) 'sum('sum((5^i10*(x-5)^(-n+i9-i10))/(i10!*(i9-i10)!),i10,0,i9),i9,0,inf)
~~~


---

Sent from sourceforge.net because max...@li... is subscribed to https://sourceforge.net/p/maxima/bugs/

To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options.  Or, if this is a mailing list, you can unsubscribe from the mailing list.

[Maxima-bugs] [maxima:bugs] #4599 intosum not fully recursive

[Barton W. <wil...@us...>]

OK, how about appending to the user documentation:

For iterated summations, users may need to apply `intosum` repeatedly to move a constant multiplicative factor inside the innermost summation. For example, for the double summation  
`a * sum(sum(f(i, j), j, j0, j1), i, i0, i1)`,  ` intosum` must be applied twice to move `a` inside both summations.




---

**[bugs:#4599] intosum not fully recursive**

**Status:** open
**Group:** None
**Labels:** intosum 
**Created:** Fri Aug 22, 2025 01:25 PM UTC by Barton Willis
**Last Updated:** Fri Aug 22, 2025 03:17 PM UTC
**Owner:** nobody


Here we need to apply `intosum` twice.  The user documentation isn't entirely clear about this, but I think that `intosum` isn't fully recursive.

~~~
	(sumexpand  : true, cauchysum  : true, display2d : false)$
	powerseries(exp(x)/(x-5)^n,x,5);
(%o37) ('sum('sum((5^i10*(x-5)^(i9-i10))/(i10!*(i9-i10)!),i10,0,i9),i9,0,inf))

 /(x-5)^n
	intosum(%);
(%o38) 'sum(('sum((5^i10*(x-5)^(i9-i10))/(i10!*(i9-i10)!),i10,0,i9))/(x-5)^n,

            i9,0,inf)
	intosum(%);
(%o39) 'sum('sum((5^i10*(x-5)^(-n+i9-i10))/(i10!*(i9-i10)!),i10,0,i9),i9,0,inf)
~~~


---

Sent from sourceforge.net because max...@li... is subscribed to https://sourceforge.net/p/maxima/bugs/

To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options.  Or, if this is a mailing list, you can unsubscribe from the mailing list.

[Maxima-bugs] [maxima:bugs] #4599 intosum not fully recursive

[Stavros M. <mac...@us...>]

Hmm. The current behavior gives the user more control -- they can *choose* to push the factor into as many levels of some as they want. If it were recursive, how could the user push it in just one level?


---

**[bugs:#4599] intosum not fully recursive**

**Status:** open
**Group:** None
**Labels:** intosum 
**Created:** Fri Aug 22, 2025 01:25 PM UTC by Barton Willis
**Last Updated:** Fri Aug 22, 2025 01:25 PM UTC
**Owner:** nobody


Here we need to apply `intosum` twice.  The user documentation isn't entirely clear about this, but I think that `intosum` isn't fully recursive.

~~~
	(sumexpand  : true, cauchysum  : true, display2d : false)$
	powerseries(exp(x)/(x-5)^n,x,5);
(%o37) ('sum('sum((5^i10*(x-5)^(i9-i10))/(i10!*(i9-i10)!),i10,0,i9),i9,0,inf))

 /(x-5)^n
	intosum(%);
(%o38) 'sum(('sum((5^i10*(x-5)^(i9-i10))/(i10!*(i9-i10)!),i10,0,i9))/(x-5)^n,

            i9,0,inf)
	intosum(%);
(%o39) 'sum('sum((5^i10*(x-5)^(-n+i9-i10))/(i10!*(i9-i10)!),i10,0,i9),i9,0,inf)
~~~


---

Sent from sourceforge.net because max...@li... is subscribed to https://sourceforge.net/p/maxima/bugs/

To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options.  Or, if this is a mailing list, you can unsubscribe from the mailing list.

[Maxima-bugs] [maxima:bugs] #4599 intosum not fully recursive

[Barton W. <wil...@us...>]

---

**[bugs:#4599] intosum not fully recursive**

**Status:** open
**Group:** None
**Labels:** intosum 
**Created:** Fri Aug 22, 2025 01:25 PM UTC by Barton Willis
**Last Updated:** Fri Aug 22, 2025 01:25 PM UTC
**Owner:** nobody


Here we need to apply `intosum` twice.  The user documentation isn't entirely clear about this, but I think that `intosum` isn't fully recursive.

~~~
	(sumexpand  : true, cauchysum  : true, display2d : false)$
	powerseries(exp(x)/(x-5)^n,x,5);
(%o37) ('sum('sum((5^i10*(x-5)^(i9-i10))/(i10!*(i9-i10)!),i10,0,i9),i9,0,inf))

 /(x-5)^n
	intosum(%);
(%o38) 'sum(('sum((5^i10*(x-5)^(i9-i10))/(i10!*(i9-i10)!),i10,0,i9))/(x-5)^n,

            i9,0,inf)
	intosum(%);
(%o39) 'sum('sum((5^i10*(x-5)^(-n+i9-i10))/(i10!*(i9-i10)!),i10,0,i9),i9,0,inf)
~~~


---

Sent from sourceforge.net because max...@li... is subscribed to https://sourceforge.net/p/maxima/bugs/

To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options.  Or, if this is a mailing list, you can unsubscribe from the mailing list.

[Maxima-bugs] [maxima:bugs] #4598 a definite integral with demoivre : false

[Barton W. <wil...@us...>]

---

**[bugs:#4598] a definite integral with demoivre : false **

**Status:** open
**Group:** None
**Labels:** defint 
**Created:** Sat Aug 16, 2025 11:22 PM UTC by Barton Willis
**Last Updated:** Sat Aug 16, 2025 11:22 PM UTC
**Owner:** nobody


OK
~~~
(%i3)	block([demoivre : false], integrate(x^8/(5+x^9)^2,x,0,inf));
(%o3)	1/45
~~~
Wrong
~~~
(%i4)	block([demoivre : true], integrate(x^8/(5+x^9)^2,x,0,inf));
(%o4)	1/135
~~~

I think this integral gets done using Cauchy's integral theorem.  


---

Sent from sourceforge.net because max...@li... is subscribed to https://sourceforge.net/p/maxima/bugs/

To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options.  Or, if this is a mailing list, you can unsubscribe from the mailing list.

[Maxima-bugs] [maxima:bugs] #4597 residue should check if taylor series is a power series

[Barton W. <wil...@us...>]

---

**[bugs:#4597] residue should check if taylor series is a power series**

**Status:** open
**Group:** None
**Created:** Sat Aug 16, 2025 10:09 PM UTC by Barton Willis
**Last Updated:** Sat Aug 16, 2025 10:09 PM UTC
**Owner:** nobody


Wrong--the residue of   `x -> log(1-x)/(1-x)` at 1 does not exist, but Maxima says it does:
~~~
(%i3)	residue(log(1-x)/(1-x),x,1);
(%o3)	-log(x-1)-log(-1)

(%i4)	build_info();
(%o4)	build_info(version="5.48.1",...)
~~~

The function `residue` assumes that `taylor` returns a power series, but for this expression it does not:
~~~
(%i6)	residue(log(1-x)/(1-x),x,1);
1" Call   "taylor[-(log(1-x)/(x-1)),x,1,1]

1" Return "taylor(-log(-1)-log(x-1)+...)/(x-1)+...

(%o6)	-log(x-1)-log(-1)
~~~



---

Sent from sourceforge.net because max...@li... is subscribed to https://sourceforge.net/p/maxima/bugs/

To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options.  Or, if this is a mailing list, you can unsubscribe from the mailing list.

[Maxima-bugs] [maxima:bugs] #4596 taylor with wrong number of arguments gives Lisp error

[Barton W. <wil...@us...>]

---

**[bugs:#4596] taylor with wrong number of arguments gives Lisp error**

**Status:** open
**Group:** None
**Labels:** taylor 
**Created:** Fri Aug 15, 2025 04:38 PM UTC by Barton Willis
**Last Updated:** Fri Aug 15, 2025 04:38 PM UTC
**Owner:** nobody


OK
~~~
(%i3)	taylor(42,x,%pi, 3);
(%o3)/T/	42+...
~~~

Not OK
~~~
(%i4)	taylor(42,x,%pi);
Maxima encountered a Lisp error:

 The value

   $TAYLOR

 is not of type

   LIST

Automatically continuing.

To enable the Lisp debugger set *debugger-hook* to nil.
~~~
This is Maxima 5.48.1


---

Sent from sourceforge.net because max...@li... is subscribed to https://sourceforge.net/p/maxima/bugs/

To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options.  Or, if this is a mailing list, you can unsubscribe from the mailing list.

[Maxima-bugs] [maxima:bugs] #4595 residue bug

[Barton W. <wil...@us...>]

---

**[bugs:#4595] residue bug**

**Status:** open
**Group:** None
**Labels:** residue 
**Created:** Thu Aug 14, 2025 03:18 PM UTC by Barton Willis
**Last Updated:** Thu Aug 14, 2025 03:18 PM UTC
**Owner:** nobody


Wrong
~~~
(%i4) residue(plog(x)/(x-%i*sin(%pi/9)-cos(%pi/9)),x,%e^((%i*%pi)/9));
(%o4) 0
~~~

OK, but should be more simple 
~~~
(%i5) residue(plog(x)/(x-%i*sin(%pi/9)-cos(%pi/9)),x,rectform(%e^((%i*%pi)/9)));
(%o5) (9*log(sin(%pi/9)^2+cos(%pi/9)^2)+2*%i*%pi)/18
~~~

Changing `plog` to `log` does not fix this problem.  

The testsuite sends this problem to `residue,` but I don't specifically know the test case that triggers this bug.

One possible fix is for the residue function to locally set  `taylor_simplifier` to `rectform`. 


---

Sent from sourceforge.net because max...@li... is subscribed to https://sourceforge.net/p/maxima/bugs/

To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options.  Or, if this is a mailing list, you can unsubscribe from the mailing list.

[Maxima-bugs] [maxima:bugs] #4594 Maxima-5.48.1 Breaks Draw3d With Surface_Hide=True

[Tomio A. <tom...@us...>]

I think that `set hidden3d nooffset` does not work expectedly on gnuplot 6.0.3.
So my workaround is to use **gnuplot 5.4.10** instead.



---

**[bugs:#4594] Maxima-5.48.1 Breaks Draw3d With Surface_Hide=True**

**Status:** open
**Group:** None
**Created:** Wed Aug 13, 2025 02:53 PM UTC by F Russell
**Last Updated:** Wed Aug 13, 2025 02:53 PM UTC
**Owner:** nobody
**Attachments:**

- [maxout8523.gnuplot](https://sourceforge.net/p/maxima/bugs/4594/attachment/maxout8523.gnuplot) (1.0 kB; application/x-gnuplot)


After upgrading to maxima-5.48.1 all of my draw3d graphics that contain the "surface_hide=true" option now broken.  Only the axes are shown.  The surface is not drawn.

The example command is:

load(draw)$
draw3d(xrange=[-2,2], yrange=[-2,2], surface_hide=true, explicit(x^2-y^2, x, -1, 1, y, -1, 1))$

I have traced the problem to the "set pm3d at s depthorder explicit" line in the maxoutXXX.gnuplot file.  This file is passed to gnuplot and the plot fails without error but no surface is shown.

If I comment https://sourceforge.net/p/maxima/bugs/4297/out the above mentioned line the surface plot will succeed.

This was supposed to have been fixed with this report:

https://sourceforge.net/p/maxima/bugs/4297/

However, the problem persists with draw3d.

I am using maxima-5.48.1 and gnuplot-6.0.3.



---

Sent from sourceforge.net because max...@li... is subscribed to https://sourceforge.net/p/maxima/bugs/

To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options.  Or, if this is a mailing list, you can unsubscribe from the mailing list.

[Maxima-bugs] [maxima:bugs] #4594 Maxima-5.48.1 Breaks Draw3d With Surface_Hide=True

[F R. <fru...@us...>]

---

**[bugs:#4594] Maxima-5.48.1 Breaks Draw3d With Surface_Hide=True**

**Status:** open
**Group:** None
**Created:** Wed Aug 13, 2025 02:53 PM UTC by F Russell
**Last Updated:** Wed Aug 13, 2025 02:53 PM UTC
**Owner:** nobody
**Attachments:**

- [maxout8523.gnuplot](https://sourceforge.net/p/maxima/bugs/4594/attachment/maxout8523.gnuplot) (1.0 kB; application/x-gnuplot)


After upgrading to maxima-5.48.1 all of my draw3d graphics that contain the "surface_hide=true" option now broken.  Only the axes are shown.  The surface is not drawn.

The example command is:

load(draw)$
draw3d(xrange=[-2,2], yrange=[-2,2], surface_hide=true, explicit(x^2-y^2, x, -1, 1, y, -1, 1))$

I have traced the problem to the "set pm3d at s depthorder explicit" line in the maxoutXXX.gnuplot file.  This file is passed to gnuplot and the plot fails without error but no surface is shown.

If I comment https://sourceforge.net/p/maxima/bugs/4297/out the above mentioned line the surface plot will succeed.

This was supposed to have been fixed with this report:

https://sourceforge.net/p/maxima/bugs/4297/

However, the problem persists with draw3d.

I am using maxima-5.48.1 and gnuplot-6.0.3.



---

Sent from sourceforge.net because max...@li... is subscribed to https://sourceforge.net/p/maxima/bugs/

To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options.  Or, if this is a mailing list, you can unsubscribe from the mailing list.

[Maxima-bugs] [maxima:bugs] #4593 tellrat called on non-monic polynomial

[Barton W. <wil...@us...>]

---

**[bugs:#4593] tellrat called on non-monic polynomial**

**Status:** open
**Group:** None
**Labels:** tellrat 
**Created:** Mon Aug 11, 2025 11:43 AM UTC by Barton Willis
**Last Updated:** Mon Aug 11, 2025 11:43 AM UTC
**Owner:** nobody


The user documentation for `tellrat` says that each argument to `tellrat` must be a "polynomial with integer coefficients,"  but it does not say that the polynomial must be monic.  But it seems that the arguments must be monic:

OK
~~~
(%i1)	tellrat(z^2 + 2);
(%o1)	[z^2+2]
~~~
Error
~~~
(%i2)	tellrat(2*x^2+1);
tellrat: minimal polynomial must be monic.
~~~
Either the user documentation needs to be changed, or the source code extended to allow non monic arguments.


---

Sent from sourceforge.net because max...@li... is subscribed to https://sourceforge.net/p/maxima/bugs/

To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options.  Or, if this is a mailing list, you can unsubscribe from the mailing list.

[Maxima-bugs] [maxima:bugs] Re: #4591 solve called on degree six gives 0^0 error

[Stavros M. <mac...@us...>]

I upgraded to 5.48.1 (using brew) and the result is identical.


---

**[bugs:#4591] solve called on degree six gives 0^0 error**

**Status:** open
**Group:** None
**Labels:** solve solvefactors 
**Created:** Tue Aug 05, 2025 10:45 AM UTC by Barton Willis
**Last Updated:** Wed Aug 06, 2025 04:40 PM UTC
**Owner:** nobody


OK:
~~~
(%i1) xxx : (x^2 + 1+%i)^3$
(%i2) solve(xxx,x);
(%o2) [x = - sqrt(- %i - 1), x = sqrt(- %i - 1)]
~~~

Not OK--this was done using the default value for solvefactors (true)
~~~
(%i3) solve(expand(xxx),x);                  
expt: undefined: 0^0
~~~




---

Sent from sourceforge.net because max...@li... is subscribed to https://sourceforge.net/p/maxima/bugs/

To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options.  Or, if this is a mailing list, you can unsubscribe from the mailing list.

[Maxima-bugs] [maxima:bugs] Re: #4591 solve called on degree six gives 0^0 error

[Stavros M. <mac...@us...>]

I'll update to 5.48.1 and see if that fixes it.
~~~~
Maxima 5.48.0 https://maxima.sourceforge.io
using Lisp SBCL 2.5.7
Distributed under the GNU Public License. See the file COPYING.
Dedicated to the memory of William Schelter.
The function bug_report() provides bug reporting information.
(%i1) ( display2d:false, leftjust:true, linel:70 )$

(%i2) xxx: x^3+3*%i*x^2+3*x^2+6*%i*x+2*%i-2$

(%i3) solve(xxx,x);
expt: undefined: 0^0
 -- an error. To debug this try: debugmode(true);

(%i4) algebraic:true;

(%o4) true
(%i5) solve(xxx,x);

(%o5) [x = (-1/2-(sqrt(3)*%i)/2)*((6*(3*%i+3)*%i-3*(2*%i-2))/6
                                 +(-1*(3*%i+3)^3)/27)
                                 ^(1/3)
         -((2*%i+(-1*(3*%i+3)^2)/9)*((sqrt(3)*%i)/2+(-1)/2))
          /((6*(3*%i+3)*%i-3*(2*%i-2))/6+(-1*(3*%i+3)^3)/27)^(1/3)
         +(-1*(3*%i+3))/3,
       x = ((sqrt(3)*%i)/2+(-1)/2)*((6*(3*%i+3)*%i-3*(2*%i-2))/6
                                   +(-1*(3*%i+3)^3)/27)
                                   ^(1/3)
         -((2*%i+(-1*(3*%i+3)^2)/9)*(-1/2-(sqrt(3)*%i)/2))
          /((6*(3*%i+3)*%i-3*(2*%i-2))/6+(-1*(3*%i+3)^3)/27)^(1/3)
         +(-1*(3*%i+3))/3,
       x = ((6*(3*%i+3)*%i-3*(2*%i-2))/6+(-1*(3*%i+3)^3)/27)^(1/3)
         -(2*%i+(-1*(3*%i+3)^2)/9)/((6*(3*%i+3)*%i-3*(2*%i-2))/6
                                   +(-1*(3*%i+3)^3)/27)
                                   ^(1/3)+(-1*(3*%i+3))/3]
(%i6) ratsimp(%);
`quotient' by `zero'
 -- an error. To debug this try: debugmode(true);
~~~~



---

**[bugs:#4591] solve called on degree six gives 0^0 error**

**Status:** open
**Group:** None
**Labels:** solve solvefactors 
**Created:** Tue Aug 05, 2025 10:45 AM UTC by Barton Willis
**Last Updated:** Tue Aug 05, 2025 10:01 PM UTC
**Owner:** nobody


OK:
~~~
(%i1) xxx : (x^2 + 1+%i)^3$
(%i2) solve(xxx,x);
(%o2) [x = - sqrt(- %i - 1), x = sqrt(- %i - 1)]
~~~

Not OK--this was done using the default value for solvefactors (true)
~~~
(%i3) solve(expand(xxx),x);                  
expt: undefined: 0^0
~~~




---

Sent from sourceforge.net because max...@li... is subscribed to https://sourceforge.net/p/maxima/bugs/

To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options.  Or, if this is a mailing list, you can unsubscribe from the mailing list.

[Maxima-bugs] [maxima:bugs] Re: #4591 solve called on degree six gives 0^0 error

[Raymond T. <rt...@us...>]

What version are you using?  With current (mostly) HEAD, I get:
```
(%i4) algebraic:true;
(%o4)                                true
(%i5) solve(expand(xxx),x);
(%o5)             [x = - sqrt(- %i - 1), x = sqrt(- %i - 1)]
```
I don't see how your answer of `-%i-1` is even a root of `xxx`.


---

**[bugs:#4591] solve called on degree six gives 0^0 error**

**Status:** open
**Group:** None
**Labels:** solve solvefactors 
**Created:** Tue Aug 05, 2025 10:45 AM UTC by Barton Willis
**Last Updated:** Tue Aug 05, 2025 06:22 PM UTC
**Owner:** nobody


OK:
~~~
(%i1) xxx : (x^2 + 1+%i)^3$
(%i2) solve(xxx,x);
(%o2) [x = - sqrt(- %i - 1), x = sqrt(- %i - 1)]
~~~

Not OK--this was done using the default value for solvefactors (true)
~~~
(%i3) solve(expand(xxx),x);                  
expt: undefined: 0^0
~~~




---

Sent from sourceforge.net because max...@li... is subscribed to https://sourceforge.net/p/maxima/bugs/

To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options.  Or, if this is a mailing list, you can unsubscribe from the mailing list.

[Maxima-bugs] [maxima:bugs] #4592 integrate causes sign: argument cannot be imaginary error

[Barton W. <wil...@us...>]

---

**[bugs:#4592] integrate causes sign: argument cannot be imaginary error**

**Status:** open
**Group:** None
**Labels:** integrate sign 
**Created:** Tue Aug 05, 2025 07:01 PM UTC by Barton Willis
**Last Updated:** Tue Aug 05, 2025 07:01 PM UTC
**Owner:** nobody


OK:
~~~
(%i12) integrate(1/(x + %i+1),x);
(%o12)                          log(x + %i + 1)
~~~
Not OK

~~~
(%i13) integrate(1/(x^4 + %i+1),x);
sign: argument cannot be imaginary; found %i
~~~



---

Sent from sourceforge.net because max...@li... is subscribed to https://sourceforge.net/p/maxima/bugs/

To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options.  Or, if this is a mailing list, you can unsubscribe from the mailing list.

[Maxima-bugs] [maxima:bugs] #4591 solve called on degree six gives 0^0 error

[Stavros M. <mac...@us...>]

~~~~
xxx : (x^2 + 1+%i)^3$
solve(expand(xxx),x)
~~~~
For all values of ``gcd``, ``solve(expand(xxx),x)`` gives an error with ``algebraic:false`` and the correct answer with ``algebraic:true``. Although the answer is simplified (i.e. ``expand(...,0,0)`` leaves it unchanged), it is ridiculously large. But ``expand`` and ``rectform`` reduce it to the simple answer.
~~~~
algebraic:true$ solve(expand(xxx),x)[1]
   =>
     - 1   sqrt(3) %i
x = (─── - ──────────)
      2        2
                                                 3
  6 (3 %i + 3) %i - 3 (2 %i - 2)   - 1 (3 %i + 3)  1/3
 (────────────────────────────── + ───────────────)
                6                        27
                              2
                - 1 (3 %i + 3)    sqrt(3) %i   - 1
        (2 %i + ───────────────) (────────── + ───)
                       9              2         2
 - ─────────────────────────────────────────────────────
                                                   3
    6 (3 %i + 3) %i - 3 (2 %i - 2)   - 1 (3 %i + 3)  1/3
   (────────────────────────────── + ───────────────)
                  6                        27
   - 1 (3 %i + 3)
 + ────────────── 
         3
----         
expand(%) => 
    x = - %i - 1
~~~~


---

**[bugs:#4591] solve called on degree six gives 0^0 error**

**Status:** open
**Group:** None
**Labels:** solve solvefactors 
**Created:** Tue Aug 05, 2025 10:45 AM UTC by Barton Willis
**Last Updated:** Tue Aug 05, 2025 11:53 AM UTC
**Owner:** nobody


OK:
~~~
(%i1) xxx : (x^2 + 1+%i)^3$
(%i2) solve(xxx,x);
(%o2) [x = - sqrt(- %i - 1), x = sqrt(- %i - 1)]
~~~

Not OK--this was done using the default value for solvefactors (true)
~~~
(%i3) solve(expand(xxx),x);                  
expt: undefined: 0^0
~~~




---

Sent from sourceforge.net because max...@li... is subscribed to https://sourceforge.net/p/maxima/bugs/

To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options.  Or, if this is a mailing list, you can unsubscribe from the mailing list.

[Maxima-bugs] [maxima:bugs] #4591 solve called on degree six gives 0^0 error

[Stavros M. <mac...@us...>]

same error for xxx:(x+%i+1)^3


---

**[bugs:#4591] solve called on degree six gives 0^0 error**

**Status:** open
**Group:** None
**Labels:** solve solvefactors 
**Created:** Tue Aug 05, 2025 10:45 AM UTC by Barton Willis
**Last Updated:** Tue Aug 05, 2025 10:45 AM UTC
**Owner:** nobody


OK:
~~~
(%i1) xxx : (x^2 + 1+%i)^3$
(%i2) solve(xxx,x);
(%o2) [x = - sqrt(- %i - 1), x = sqrt(- %i - 1)]
~~~

Not OK--this was done using the default value for solvefactors (true)
~~~
(%i3) solve(expand(xxx),x);                  
expt: undefined: 0^0
~~~




---

Sent from sourceforge.net because max...@li... is subscribed to https://sourceforge.net/p/maxima/bugs/

To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options.  Or, if this is a mailing list, you can unsubscribe from the mailing list.

[Maxima-bugs] [maxima:bugs] #4591 solve called on degree six gives 0^0 error

[Barton W. <wil...@us...>]

---

**[bugs:#4591] solve called on degree six gives 0^0 error**

**Status:** open
**Group:** None
**Labels:** solve solvefactors 
**Created:** Tue Aug 05, 2025 10:45 AM UTC by Barton Willis
**Last Updated:** Tue Aug 05, 2025 10:45 AM UTC
**Owner:** nobody


OK:
~~~
(%i1) xxx : (x^2 + 1+%i)^3$
(%i2) solve(xxx,x);
(%o2) [x = - sqrt(- %i - 1), x = sqrt(- %i - 1)]
~~~

Not OK--this was done using the default value for solvefactors (true)
~~~
(%i3) solve(expand(xxx),x);                  
expt: undefined: 0^0
~~~




---

Sent from sourceforge.net because max...@li... is subscribed to https://sourceforge.net/p/maxima/bugs/

To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options.  Or, if this is a mailing list, you can unsubscribe from the mailing list.