Maple 2024 Questions and Posts

Error, (in RootOf) _Z occurs but is not the depend...

Asked by: nm 11353

July 07 2024

2 3

I added radnormal(sol) to my solution to workaround bug in solve hanging.

But now new problem showed up. sometimes radnormal gives internal error when there are _Z's in solution.

radnormal(sol);
Error, (in RootOf) _Z occurs but is not the dependent variable

Attached worksheet. Sorry that the solution is very large and has lots of _Zs and RootOf, but this is the first one I can see so far in the log file of my program running, so I left it as is:

Should I check in my code that solution does not contain _Z before calling radnormal on it? Is this a bug or known limitation?

>

restart;

>	interface(version);

>	Physics:-Version();

>

sol:=1/6*(-a^3 - 3*RootOf(4*_Z^2 - 4*_Z*(8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(1/3) + 2*a*_Z + (8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(2/3) - a*(8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(1/3) + a^2)*a^2 + 6*RootOf(4*_Z^2 - 4*_Z*(8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(1/3) + 2*a*_Z + (8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(2/3) - a*(8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(1/3) + a^2)^2*a + 8*RootOf(4*_Z^2 - 4*_Z*(8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(1/3) + 2*a*_Z + (8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(2/3) - a*(8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(1/3) + a^2)^3 + 3*sqrt(3)*sqrt(-RootOf(4*_Z^2 - 4*_Z*(8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(1/3) + 2*a*_Z + (8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(2/3) - a*(8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(1/3) + a^2)*(RootOf(4*_Z^2 - 4*_Z*(8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(1/3) + 2*a*_Z + (8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(2/3) - a*(8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(1/3) + a^2)*a^4 + 4*RootOf(4*_Z^2 - 4*_Z*(8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(1/3) + 2*a*_Z + (8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(2/3) - a*(8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(1/3) + a^2)^2*a^3 + 4*RootOf(4*_Z^2 - 4*_Z*(8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(1/3) + 2*a*_Z + (8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(2/3) - a*(8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(1/3) + a^2)^3*a^2 + 4*a^3 + 12*RootOf(4*_Z^2 - 4*_Z*(8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(1/3) + 2*a*_Z + (8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(2/3) - a*(8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(1/3) + a^2)*a^2 - 24*RootOf(4*_Z^2 - 4*_Z*(8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(1/3) + 2*a*_Z + (8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(2/3) - a*(8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(1/3) + a^2)^2*a - 32*RootOf(4*_Z^2 - 4*_Z*(8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(1/3) + 2*a*_Z + (8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(2/3) - a*(8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(1/3) + a^2)^3 - 108*RootOf(4*_Z^2 - 4*_Z*(8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(1/3) + 2*a*_Z + (8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(2/3) - a*(8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(1/3) + a^2))) + 54*RootOf(4*_Z^2 - 4*_Z*(8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(1/3) + 2*a*_Z + (8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(2/3) - a*(8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(1/3) + a^2))^(1/3) + 1/6*(4*RootOf(4*_Z^2 - 4*_Z*(8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(1/3) + 2*a*_Z + (8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(2/3) - a*(8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(1/3) + a^2)^2 + 2*a*RootOf(4*_Z^2 - 4*_Z*(8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(1/3) + 2*a*_Z + (8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(2/3) - a*(8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(1/3) + a^2) + a^2)/(-a^3 - 3*RootOf(4*_Z^2 - 4*_Z*(8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(1/3) + 2*a*_Z + (8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(2/3) - a*(8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(1/3) + a^2)*a^2 + 6*RootOf(4*_Z^2 - 4*_Z*(8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(1/3) + 2*a*_Z + (8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(2/3) - a*(8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(1/3) + a^2)^2*a + 8*RootOf(4*_Z^2 - 4*_Z*(8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(1/3) + 2*a*_Z + (8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(2/3) - a*(8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(1/3) + a^2)^3 + 3*sqrt(3)*sqrt(-RootOf(4*_Z^2 - 4*_Z*(8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(1/3) + 2*a*_Z + (8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(2/3) - a*(8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(1/3) + a^2)*(RootOf(4*_Z^2 - 4*_Z*(8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(1/3) + 2*a*_Z + (8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(2/3) - a*(8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(1/3) + a^2)*a^4 + 4*RootOf(4*_Z^2 - 4*_Z*(8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(1/3) + 2*a*_Z + (8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(2/3) - a*(8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(1/3) + a^2)^2*a^3 + 4*RootOf(4*_Z^2 - 4*_Z*(8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(1/3) + 2*a*_Z + (8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(2/3) - a*(8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(1/3) + a^2)^3*a^2 + 4*a^3 + 12*RootOf(4*_Z^2 - 4*_Z*(8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(1/3) + 2*a*_Z + (8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(2/3) - a*(8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(1/3) + a^2)*a^2 - 24*RootOf(4*_Z^2 - 4*_Z*(8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(1/3) + 2*a*_Z + (8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(2/3) - a*(8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(1/3) + a^2)^2*a - 32*RootOf(4*_Z^2 - 4*_Z*(8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(1/3) + 2*a*_Z + (8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(2/3) - a*(8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(1/3) + a^2)^3 - 108*RootOf(4*_Z^2 - 4*_Z*(8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(1/3) + 2*a*_Z + (8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(2/3) - a*(8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(1/3) + a^2))) + 54*RootOf(4*_Z^2 - 4*_Z*(8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(1/3) + 2*a*_Z + (8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(2/3) - a*(8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(1/3) + a^2))^(1/3) - 1/6*a + 1/3*RootOf(4*_Z^2 - 4*_Z*(8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(1/3) + 2*a*_Z + (8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(2/3) - a*(8*_Z^3 + 6*_Z^2*a - 3*_Z*a^2 - a^3 + 3*sqrt(3)*sqrt(-4*_Z^4*a^2 - 4*_Z^3*a^3 - _Z^2*a^4 + 32*_Z^4 + 24*_Z^3*a - 12*_Z^2*a^2 - 4*_Z*a^3 + 108*_Z^2) + 54*_Z)^(1/3) + a^2):

>	radnormal(sol);

Error, (in RootOf) _Z occurs but is not the dependent variable

>

Download bug_Z.mw

2024.1 freezing or partial lockout...

Asked by: Ronan 1341

July 06 2024

1 12

I have had this a few times this week since updating to 2024.1 on Windows 10.

I get sudden freezes in a worksheet. The !!! button greys out. The ! button is ok, so the worksheet can be run by using ctrl A and click !

Has anyone else experienced this?

Export Array/Table/Dataframe of data to JPG...

Asked by: Scot Gould 1039

July 04 2024

0 0

I have a Dataframe of data, although I assume this question applies to any type of rTable-like structure.

What is a simple/elegant way to export the image of the data to a JPG file? I would be happy to see it in the format when I ask it to print the Dataframe, or when I use DocumentTools:-Tabulate.

Why does the interface command performs differentl...

Asked by: C_R 3402

June 28 2024

1 3

Maple formats output depending on typesetting options "extened" and "standard" for the GUI (or interface). An example taken from

restart;
ts_standard:=proc(k::anything)
     interface(typesetting=standard):
     print(k);
     interface(typesetting=extended): 
     NULL;
end proc:
k:=3/8*ln(55/52)+sin(x)+3/4*exp(x);
                    3   /55\            3       
               k := - ln|--| + sin(x) + - exp(x)
                    8   \52/            4       

ts_standard(k);

Why is the input two times returned and why one time as a list?
Somehow the first interface statement is responsible for that.

I am only interested in the reformated input inside the list.
Is it possible to fix the code?

Other observation with typesetting=standard:

restart;
interface(typesetting=standard);
expr:=cos(x)^2;
((x->x)=combine[trig])(expr);
                      

((x->x)=combine[trig])(expr);#subsequent call
                          2   1            1
                    cos(x)  = - cos(2 x) + -
                              2            2

(with Maple 2024.1 this only accurs in Math-1D)
but

restart;
expr:=cos(x)^2;
interface(typesetting=standard);
((x->x)=combine[trig])(expr);
                                      2
                        expr := cos(x) 

                            extended

                          2   1            1
                    cos(x)  = - cos(2 x) + -
                              2            2

Maple 2024 crash calling solve. Reproducible each ...

Asked by: nm 11353

June 18 2024

1 5

I am getting Maple server crash each time running this solve command.

Could others reproduce it? I am using windows 10. Maple 2024. Why does it happen?

Will report it to Maplesoft in case it is not known. Worksheet below.

Download maple_crash_calling_solve_june_18_2024.mw

This bug seems to have been introduced in Maple 2023 since it crashes there also.

But not in Maple 2022. No crash there. Same PC.

Download maple_NO_crash_calling_solve_june_18_maple_2022.mw

to simplify or not simplify before calling integra...

Asked by: nm 11353

June 15 2024

1 2

Is this a valid behvior by int?

int(A,x,method=_RETURNVERBOSE) hangs.

But int(simplify(A),x,method=_RETURNVERBOSE) returns in few seconds with "default" result same as int(A,x)

Should this have happen? I try to avoid calling simplify unless neccessary because it can add csgn's and signums and so on to the result.

But the question is: Should one really need to simplify the integrand to get the result in this example? Does this mean one should call simplify on the integrand to avoid the hang that can show up?

This only happens when using method=_RETURNVERBOSE

Just trying to find out if this is normal behavior and can be expected sometimes.

Download why_int_hang_unless_simplify_june_15_2024.mw

Where is the code that simplifies this ...

Asked by: janhardo 685

June 08 2024

1 43

Can't figure out what code makes this simplification.
If this simplification works, it will be a part of a larger simplication procedure ( if it not conflicts hopefully)
vereenvouding_hoe_-vraag_MPF.mw

PDEtools:-Solve vs. solve. RootOf form difference...

Asked by: nm 11353

June 08 2024

1 4

I was trying to find out why my solution was not validating for this ode. It turned out because I was using solve instead of PDEtools:-Solve. It took me sometime to find this.

This made huge difference on odetest to verify the solution.

This is very simple ode. We just need to integrate once. But first we have to solve for y'(x).

And here comes the difference. When I used solve to solve for y'(x), odetest did not verify the solution.

When using PDEtools:-Solve, it did.

The difference is how each returned the solution for y'(x). Both have RootOf but written differently and this made the difference.

1) Why solutions are written differently?

2) Is this to be expected? I have thought Solve uses same engine as solve below the cover.

3) is it possible to make solve give the same form as Solve or change to that form?

I am now changing more of my code to use PDEtools:-Solve because of this.

>	interface(version);

>	Physics:-Version();

$`The "Physics Updates" version in the MapleCloud is 1757. The version installed in this computer is 1756 created 2024, June 5, 19:39 hours Pacific Time, found in the directory C:\Users\Owner\maple\toolbox\2024\Physics Updates\lib\`$

Using solve

>

restart;

>	ode:=x-ln(diff(y(x),x))-sin(diff(y(x),x))=0; RHS:=solve(ode,diff(y(x),x));

>	mysol:= y(x) = Int(RHS,x)+c__1;

>	odetest(mysol,ode);

using PDEtools:-Solve (now it verifies) with no extra effort

>

restart;

>	ode:=x-ln(diff(y(x),x))-sin(diff(y(x),x))=0; RHS:=PDEtools:-Solve(ode,diff(y(x),x)): RHS:=rhs(%);

>	mysol:= y(x) = Int(RHS,x)+c__1;

>	odetest(mysol,ode);

Download PDEtools_Solve_vs_solve_june_8_2024.mw

Update

Here is a counter example. Where now it is the other way around.

Using solve makes odetest happy, but when using PDEtools:-Solve odetest does not verify the solution. Same exact ODE.

>	interface(version);

>	Physics:-Version()

Example, using solve works

>	ode:=exp(diff(y(x), x) - y(x)) - diff(y(x), x)^2 + 1 = 0; RHS:=solve(ode,diff(y(x),x)); RHS:=eval(RHS,y(x)=y); mysol:=Intat(eval(1/RHS,y=_a),_a=y(x))=x+c__1; odetest(mysol,ode);

exp(diff(y(x), x)-y(x))-(diff(y(x), x))^2+1 = 0

Warning, solutions may have been lost

Intat(1/RootOf(-exp(_Z-_a)+_Z^2-1), _a = y(x)) = x+c__1

Example, using PDEtools:-Solve fails

>	ode:=exp(diff(y(x), x) - y(x)) - diff(y(x), x)^2 + 1 = 0; RHS:=rhs(PDEtools:-Solve(ode,diff(y(x),x))); RHS:=eval(RHS,y(x)=y); mysol:=Intat(eval(1/RHS,y=_a),_a=y(x))=x+c__1; odetest(mysol,ode);

exp(diff(y(x), x)-y(x))-(diff(y(x), x))^2+1 = 0

Intat(1/RootOf(_Z^2*exp(_a)-exp(_Z)-exp(_a)), _a = y(x)) = x+c__1

>

Download PDEtools_Solve_vs_solve_june_9_2024.mw

So now I have no idea which to use. Sometimes solve works and sometimes Solve works. I guess I have to now solve the ode both ways each time and see which works.

should int() result match exactly what is given as...

Asked by: nm 11353

June 01 2024

2 5

Simple question. I hope it has simple answer. I have always thought that what int() returns should match exactly what "default" result shows when using int() with the option _RETURNVERBOSE

I mean exact match. But this below shows that int() result underwent some simplification as it is not the same as default.

restart;

integrand:=sin(x)/(sin(x) + 1);
maple_result_1 :=int(integrand,x);
maple_result_2 := int(integrand,x,'method'=':-_RETURNVERBOSE')[1]

Ofcourse maple_result_2 can be made the same as maple_result_1

simplify(rhs(maple_result_2)) assuming 0<x and x<Pi;

But this is beside the point. Why is "default" is not excatly the same as int() result? It seems that int() does something more after obtaining the :"default" result as shown.

Should default not match exactly result from int() ?

Maple 2024.

Numeric Formatting Causes Maple to Hang...

Asked by: jganding 135

May 31 2024

0 4

I put together the attached worksheet to help me determine the cheapest way to buy "refreshments" for a party by comparing price and volume of different bottle size options. The spreadsheet works fine as is. However, when I right click on the output of line (14) and format pct_difference as percent with 2 decimal places and execute the worksheet, Maple hangs on that line and progresses no further. This doesn't happen in Maple 2018 but the problem does show up in Maple 2024. Suggestions please?

cost_comparison_-_liquid_(v01MP).mw

Why does this command only works in Maple and not ...

Asked by: FDS 190

May 24 2024

1 2

I am stuck this command works seemlessly in Maple:

ThermophysicalData:-CoolProp:-Property(D, T = 20*Unit('degC'), P = 760*Unit('mmHg'), water)

but it does not work in Maple Flow. Does anyone knows why? Thank you so much for your help in the matter.

How come does "residualQQplot" give me a strange ...

Asked by: Mathrookie94 20

May 23 2024

1 12

Hello :)

I have a math problem, where I first need to use Linear regression to find the equation based on a set of data. I did that, no problem.

However, in the next part of the problem I need to check if the residuals are under "normal distribution". Usually, I check if a dataset is normally distributed via "QQ-plot", and there will be no problems. But this time, because I need to check the residuals, I need to use the "residualQQplot(data,LinReg)" command to make it happen. But when I read the mean-value, mu, it says "-0," and nothing else? I know it should be "-3,2752*10^-15.

The standard deviation is correct.

How do I fix this, so the residualQQplot shows me the right result?

I have attached the worksheet here. worksheet_-_linear_reg_and_residuals_for_normal_distribution.mw

Thank you!

Maple gives solutions may have been lost unless tw...

Asked by: nm 11353

May 22 2024

1 5

I can't understand this behavior. Any idea why it happens?

Solve is able to solve equation f(y)=x+A for y, but can't solve f(y)=x for y.

This is unexpected for me. I do not see why it can solve it when RHS is x+A but not when RHS is just x.

I can trick it to solve f(y)=x for y by asking it to solve f(y)=x+A for y and then set A=0 in the solution. But one should not have to do this. Is this a bug or Am I missing something?

Download why_solve_when_adding_term_only_may_22_2024.mw

warning (`rtable_alias`)...

Asked by: charlie_fcl 130

May 18 2024

1 12

Hi,

I'm trying out the 2024 version of Maple and I'm getting the following warning message:

Warning, not a built-in function (`rtable_alias`)

which I didn´t get for the 2023 version. I have no clue where it is coming from since it happens even when I start a new worksheet:

I've also attached print outs of the same worksheets (from Maple help examples and from Maple Portal), one using Maple 2023 version and the other one using Maple 2024 version so youcould see the warning and some other problems.

I really appreciate if someone would have an idea of what is going on here. Thanks very much in advance.

interpolation_2023.pdf

interpolation_2024.pdf

optimization_2023.pdf

optimization_2024.pdf

File open not working on Windows system...

Asked by: C_R 3402

May 14 2024

1 3

This is new:

Maple 2024 frozen on opening recent files

Maple 2023 frozen on opening

Maple 2022 frozen on opening start page

Maple 2021 blinking

Maple 2020 opening start page

The above system state is constant for about 30 min. Maple sessions without start page are working. I can enter code but file opeing and saving does not work. The fact that Maple 2020 is also not working makes it unlikely that the Java environement is part of the problem.

I have several times restarted the system. The rest of the system is working.

Something happened to the system and I have no clue what is was and what I can do about it.

Any ideas or suggestions what I could try? Windows 10.

E-Mail Address:
Password:
Remember Me:	Automatically sign in on future visits

E-Mail Address:
Password:
Remember Me:	Automatically sign in on future visits

Ask a Question

Create a Post