Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

I called dsolve with timelimit on this ode.  All of Maple instantly locks up.

I do not mean the worksheet, but everything. I have each worksheet using its own engine,. Not able to open new worksheet, nothing clicks. The whole front end locks up. Can't click on anything. 

Killing every mserver.exe does not even resolve this. I had to terminate all of Maple from the task manager. The strange thing, is that looking at task manager I see mserver.exe doing nothing. zero CPU. Only the front end process is at high CPU (the one with the Java icon).  It looks like the Java frontend is locked up for some reason.

 

Do others see the same thing? does this happen on the mac also? Please make sure to save all your work before trying to run this on your PC.

Maple 2024.1 on windows 10. This is first order Riccati ode.

interface(version);

`Standard Worksheet Interface, Maple 2024.1, Windows 10, June 25 2024 Build ID 1835466`

Physics:-Version();

`The "Physics Updates" version in the MapleCloud is 1774 and is the same as the version installed in this computer, created 2024, July 16, 17:18 hours Pacific Time.`

ode:=(a2*x^2+b2*x+c2)*diff(y(x),x)=y(x)^2+(a1*x+b1)*y(x)+a0*x^2+b0*x+c0;

(a2*x^2+b2*x+c2)*(diff(y(x), x)) = y(x)^2+(a1*x+b1)*y(x)+a0*x^2+b0*x+c0

DEtools:-odeadvisor(ode);

[_rational, _Riccati]

#WARNING. This will freeze all of Maple.
timelimit(30,dsolve(ode));


Download maple_lock_up_july_17_2024.mw

These are GUI options on my end.

 

 

 

 

UPDATE 1

in Maple 2023.2 also on windows 10,  there is no lock up at all. same code. The solution to the ode is very large. I tried using both typesetting level  to EXTENDED  and  typesetting level MAPLE STANDARD and no hang.

I wonder if this has something to do with why Maple 2024.1 locks up? But I have not changed anything in my end. I made no change in options or anything else. same PC, same graphics card. 

UPDATE 2

Found the BUG!!  

In Maple 2024.1, if I change the display option typesetting level  to EXTENDED, the front end do not hang.

If I change the display typesetting level MAPLE STANDARD, the ftont end hangs.

But why??


 

interface(version);

`Standard Worksheet Interface, Maple 2024.1, Windows 10, June 25 2024 Build ID 1835466`

Physics:-Version();

`The "Physics Updates" version in the MapleCloud is 1774 and is the same as the version installed in this computer, created 2024, July 16, 17:18 hours Pacific Time.`

ode:=(a2*x^2+b2*x+c2)*diff(y(x),x)=y(x)^2+(a1*x+b1)*y(x)+a0*x^2+b0*x+c0;

(a2*x^2+b2*x+c2)*(diff(y(x), x)) = y(x)^2+(a1*x+b1)*y(x)+a0*x^2+b0*x+c0

DEtools:-odeadvisor(ode);

[_rational, _Riccati]

interface(typesetting=extended):

timelimit(30,dsolve(ode));

`[Length of output exceeds limit of 1000000]`

interface(typesetting=standard):

#WARNING, this will hang MAPLE now. WHY??
timelimit(30,dsolve(ode));


 

Download maple_lock_up_july_17_2024_2.mw

ps. Reported to Maplesoft support. July 17, 2024. Hopefully this will be fixed in Maple 2024.2

 

Good day, all.

Please I want to solve the following delay differential equation:

ODE := diff(y(t), t$2) = (2*(1-y(t-1)^2))*(diff(y(t), t))-y(t)

ics := y(0) = 1, (D(y))(0) = 0

using the following codes but there is an error. Please kindly help to modify the codes.

restart;
Digits:=30:

f:=proc(n)
	2*(1-(y[n-1])^2)*delta[n]+y[n]:
end proc:

g:=proc(n)
	-4*y[n-1]*delta[n-1]+2*(1-(y[n-1])^2)*f(n)-delta[n]:
end proc:


e1:=y[n+2] = -y[n]+2*y[n+1]+(1/120)*h^2*(-3*h*g(n+2)+3*g(n)*h+16*f(n+2)+16*f(n)+88*f(n+1)):
e2:=h*delta[n] = -y[n]+y[n+1]-(1/1680)*h^2*(-128*h*g(n+1)-11*h*g(n+2)+59*g(n)*h+40*f(n+2)+520*f(n)+280*f(n+1)):
e3:=h*delta[n+1] = -y[n]+y[n+1]+(1/1680)*h^2*(-152*h*g(n+1)-10*h*g(n+2)+32*g(n)*h+37*f(n+2)+187*f(n)+616*f(n+1)):
e4:=h*delta[n+2] = -y[n]+y[n+1]+(1/1680)*h^2*(128*h*g(n+1)-101*h*g(n+2)+53*g(n)*h+744*f(n+2)+264*f(n)+1512*f(n+1)):

inx:=0:
ind:=0:
iny:=1:
h:=1/2:
n:=1:
omega:=10:
u:=omega*h:
N:=solve(h*p = 10, p):

err := Vector(round(N)):
exy_lst := Vector(round(N)):
numerical_y1:=Vector(round(N)):

c:=1:
for j from 0 to 2 do
	t[j]:=inx+j*h:
end do:

vars:=y[n+1],y[n+2],delta[n+1],delta[n+2]:

step := [seq](eval(x, x=c*h), c=1..N):
printf("%6s%45s%45s\n", 
	"h","Num.y","Num.z");
#eval(<vars>, solve({e||(1..4)},{vars}));


st := time():
for k from 1 to N/2 do

	par1:=x[0]=t[0],x[1]=t[1],x[2]=t[2]:
	par2:=y[n]=iny,delta[n]=ind:
	res:=eval(<vars>, fsolve(eval({e||(1..4)},[par1,par2]), {vars}));

	for i from 1 to 2 do
		printf("%6.5f%45.30f%45.30f\n", 
		h*c,res[i],res[i+2]):
		
		numerical_y1[c] := res[i]:
		
		c:=c+1:
	end do:
	iny:=res[2]:
	ind:=res[4]:
	inx:=t[2]:
	for j from 0 to 2 do
		t[j]:=inx + j*h:
	end do:
end do:
v:=time() - st;
v/4;
printf("Maximum error is %.13g\n", max(err));
NFE=evalf((N/4*3)+1);
#get array of numerical and exact solutions for y1
numerical_array_y1 := [seq(numerical_y1[i], i = 1 .. N)]:
#exact_array_y1 := [seq(exy[i], i = 1 .. N)]:

#get array of time steps
time_t := [seq](step[i], i = 1 .. N):

#display graphs for y1
with(plots):
numerical_plot_y1 := plot(time_t, numerical_array_y1, style = [point], symbol = [asterisk],
				color = [blue,blue],symbolsize = 20, legend = ["TFIBF"]);

 Thank you, and best regards.

Dear all,

I'd like to explore graphically a polynomial surface depending on two parameters a and b.

The problem is that, as soon as I start playing with the sliders, Maple freezes and I have to 'force quit'.

Can you please tell me if you have the same problem with this example?

Thanks.

restart:

with(plots):

K := 1 - y*x - (1 - x)*(b*x^3 + a*x^2 + x + 1)*(1 - y)*(b*y^3 + a*y^2 + y + 1)

1-y*x-(1-x)*(b*x^3+a*x^2+x+1)*(1-y)*(b*y^3+a*y^2+y+1)

(1)

Explore(plot3d(K, x=0..1, y=0..1,  font= [Times, bold, 20], labels= ['x', 'y','z'],   labelfont= [Times, bold, 40],  title = "K(x,y)"), b = 0..1., a=0..1.);

 

 

Download Explore.mw

Hello Dear Professional users,

I have a question regarding the "fsolve" command and also the "assign" command in Maple.

In my previous codes, I just used one time from "fsolve" and then "assign" command.

Today, I want to use "fsolve" and "assign" in a for loop. But, I can not get the results correctly.

For example, previously I reach a system of algebraic equations and the solve my problem easily:

N:=8:
y:=sum(a[n]*t^n,n=0..N):
y:=unapply(y,t):
*****some calculations with a[n] coefficients *****
A:={a set of nonlinear algebraic equations in terms of a[n]}%The number of equations is N+1 (same as the number of a[n] for n=0..N)
sol:=fsolve(A):
assign(sol)
plot(y(t),t=0..1)

****************************************************************************************************
****************************************************************************************************
****************************************************************************************************

Today, my problem is:

N:=8:
M:=4:
for i from 1 to M do
y[i]:=sum(a[i,n]*t^n,n=0..N):
y[i]:=unapply(y[i],t):
end do
for i from 1 to M do
*****some calculations with a[n] coefficients *****
@A[i]:={a set of nonlinear algebraic equations in terms of a[i,n]}%The number of equations is N+1 (same as the number of a[n] for n=0..N)
@sol[i]:=fsolve(A[i]):
@assign(sol[i])
@plot(y[i](t),t=0..1)
end do
****************************************************************************************************
****************************************************************************************************
****************************************************************************************************

What I should write instead of "A[i]", "sol[i]", and "assign(sol[i])" the lines started with
@A[i],....
@sol[i],...
@assign(sol[i]),....

Thanks for your attention in advance

With kind regards,
Emran

I am trying to install Syrup in my home computer (I have it installed in my work computer).  I followed the instructions in the Readme file: 

From Standard Maple:
    Open the file Syrup-Installer.mla.
    To do so, use File -> Open, choose file type 
    "Maple Library Archive (.mla)", select the file, and
    click "Open".

Everything seemed to work and the help page opened up:  But, it is not the syrup help page.  furthermore, when I type ?Syrup, it doesn't open it either.

 

When I try to run a worksheet that uses Syrup (that works on my work computer), I get these errors:

I"m going to reboot now and try again.  

Jorge

I was working with a Dataframe when I wanted to change the datatype of multiple columns at the same time as this is quite a large dataframe. I found in the helpfile that I can change datatype by the following command: 

SubsDatatype(Data, plts, float) which then change the datatype of "plts" into float. I had hoped that using multiple columns in the command would work in this way: SubsDatatype(Data, [plts, act], float)  but apparently not. Is there a way to do this or do I have to do it column by column?

Additionally I have another question about dataframes. I would like to replace "0" in the dataframe by a "blank" as you can do in excel. How do you do this in a dataframe?

Thanks in advance for any help given!

there is any way for define conformable fractional derivative in partial differential equation

restart;
with(PDEtools);
pde := a*diff(psi(x, t), x $ 2) + (b*abs(psi(x, t))^(-2*n) + c*abs(psi(x, t))^(-n) + d*abs(psi(x, t))^n + f*abs(psi(x, t))^(2*n))*psi(x, t) = 0;
pde + i*diff(u(x, t), [t $ beta]) = 0;

how define a  fractional derivative in sense of conformable derivative

I notice that the command gcd(a,b), if a and b are large degree polynomials, takes too much time and often crashes Windows (not only Maple).

As the euclidean algorithm is very efficient even for large numbers,why not for polynomials?

And how could I calculate the gcd between polynomials with a large degree?

Thanks Michele

i did a solution of this ODE equation but the solution of paper is different from mine also in other some equation i have same problem i can't get exactly and pretty solution

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);

`Standard Worksheet Interface, Maple 2024.1, Windows 10, June 25 2024 Build ID 1835466`

Physics:-Version();

`The "Physics Updates" version in the MapleCloud is 1767 and is the same as the version installed in this computer, created 2024, June 28, 12:19 hours Pacific Time.`

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

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?

Dear maple user,for defining the piecewise function please rectify this

h:z-> piecewise(do+Lo<z<do+4*Lo+0.3,    1-cos(2*pi*(z-L), other wise 1)

I have four symmetries as

1] \frac{\partial}{\partial t)

2] \frac{\partial}{\partial x)

3] \frac{\partial}{\partial y)

4] 2t \frac{\partial}{\partial t)+x\frac{\partial}{\partial x)+y\frac{\partial}{\partial y)-2u\frac{\partial}{\partial u)

Kindly help me out to find 1D optimal system with structural constants.

I will be greatful.

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.

I have created two graphs*, each is a combination of multiple plots. From this, I use "plots:-display" to create an array of the two graphs.  I wish to export this array of graphs to a jpg file. Is this possible?  (Yes, I can export each graph, but I would like it combined in the array image that Maple displays.)

* Why the asterisk: I make two plots. When I combine them through the plots:-display procedure, what type of output is produced? I've called it a graph. The help implies the output is a plot array. However, that is different from displaying an array of graphs through plots: display (Array([plot1, plot2])). Correct?

First 17 18 19 20 21 22 23 Last Page 19 of 2215