acer

Physics, PDEtools:-undeclare...

Answered: acer 32440

17 hours ago

1 1

The hang (in Maple 2024 at least), seems to depend on having both

with(Physics):

and that,

PDEtools:-undeclare(...):

In fact it seems related just to the 2D prettyprinting of the constructed partial derivatives within the Vx Matrix result under default extended typesetting (and that loaded environment).

(You can check that claim by suppressing the relevant output by making the problematic line end with a full colon as statement terminator. You can also lprint it.)

So, one workaround is to not make that undeclare call.

I didn't check whether the latest Physics update add-on the Maple 2024 fixes it.

Here is a shorter example that illustrates the hang (which can be avoided with several kinds of workaround...),

>

diff(diff(omega(x,t),t),x)

>	Printing this hangs

>

Download derivativve_ex01.mw

collect, normal...

Answered: acer 32440

18 hours ago

0 6

nb. You could also get by here with just,

collect(F2,[theta1],normal)

The normal call smartens up those negative signs in the "coefficient".

Download janh_01.mw

Scale...

Answered: acer 32440

Yesterday at 6:30 PM

1 5

You could scale your image to some desired size (in pixels).

Img_Scale.mw

ps. If you create your Matrix M with the datatype=float[8] option then you could Embed/Scale that, without having to convert to `Image`.

sequence of two expressions...

Answered: acer 32440

Yesterday at 2:29 AM

2 1

Your result from solve is an expression sequence (ie. a sequence of two expressions).

You can't just pass that sequence straight to the simplify command, as the arguments after the first would get interpreted as incorrectly specifying options. That's why you got that error message.

Instead, you could put those two results from solve into a list, and then simplify that.

Download Q_Simplify_ac.mw

ScrollableMathTableOutput = false...

Answered: acer 32440

September 21 2025

1 0

Your worksheet contains the input,

ScrollableMathTableOutput = false

but that would do nothing as Maple input. It's not an assignment and, more importantly, there is no such programmatic option. It can't do anything relevant here.

There can be a meaningful such line, if you add it using an external text editor (preferably while the Maple GUI is not running), in the GUI's stored configuration/preferences file. That is not a file of Maple commands. It's a file of stored GUI options/preferences.

I don't know whether that would resolve your problem. But note that your worksheet did not, in and of itself, disable matrix scrolling. So you may well not have actually successfully tried that yet.

ps. Your worksheet attachment was last saved using Maple 2025.0. It may not fix this particular issue, but you might otherwise benefit from the Maple 2025.1 point-release update.

the changed behavior...

Answered: acer 32440

September 21 2025

1 0

The relevant thing that changed in Maple's behavior is that lowercase int began to try floating-point quadrature (aka numeric integration) in the case that it detects floats in the integrand. (Previously, it did that if only the numeric range contained floats.) This happens to have exposed a weakness in a particular numeric integration method.

So now your example (with the float in the integrand) is attempting numeric rather than symbolic integration, and it goes wrong using the _d01amc method.

Note however that you can disable that computational route, for the int command here, by supplying the numeric=false option. Then you'd get successful symbolic integration here, even for this example with a float in the integrand.

You could also get successful numeric integration here by forcing some other method (or raising Digits).

I believe that covers the relevant details. You can see the commonalities below, run with Maple 2022.2 versus Maple 2024.2. It is only the first variant below which changed, which behaves like the OP's example.

>

restart;

>	kernelopts(version);

>	infolevel[`evalf/int`]:=1:

>	int(1e6exp(-1e6t), t=0..infinity);

>	infolevel[`evalf/int`]:=1:

>	int(1e6exp(-1e6t), t=0.0..infinity);

evalf/int/control: attempting Levin method

Control: Entering NAGInt

trying d01amc (nag_1d_quad_inf)

d01amc: result=0.

>

restart;

>	infolevel[`evalf/int`]:=1:

>	evalf(Int(1e6exp(-1e6t), t=0..infinity));

evalf/int/control: attempting Levin method

Control: Entering NAGInt

trying d01amc (nag_1d_quad_inf)

d01amc: result=0.

>

restart;

>	infolevel[`evalf/int`]:=1:

>	int(1e6exp(-1e6t), t=0.0..infinity, numeric=false);

>

restart;

>	map(M->int(1e6exp(-1e6t), t=0.0..infinity, method=M), [_d01amc,gquad,_Dexp,_d01ajc]);

>	restart; Digits:=16:

>	map(M->int(1e6exp(-1e6t), t=0.0..infinity, method=M), [_d01amc,gquad,_Dexp,_d01ajc]);

Download int_deta_2022.2.mw

>

restart;

>	kernelopts(version);

>	infolevel[`evalf/int`]:=1:

>	int(1e6exp(-1e6t), t=0..infinity);

evalf/int/control: attempting Levin method

Control: Entering NAGInt

trying d01amc (nag_1d_quad_inf)

d01amc: result=0.

>	infolevel[`evalf/int`]:=1:

>	int(1e6exp(-1e6t), t=0.0..infinity);

evalf/int/control: attempting Levin method

Control: Entering NAGInt

trying d01amc (nag_1d_quad_inf)

d01amc: result=0.

>

restart;

>	infolevel[`evalf/int`]:=1:

>	evalf(Int(1e6exp(-1e6t), t=0..infinity));

evalf/int/control: attempting Levin method

Control: Entering NAGInt

trying d01amc (nag_1d_quad_inf)

d01amc: result=0.

>

restart;

>	infolevel[`evalf/int`]:=1:

>	int(1e6exp(-1e6t), t=0.0..infinity, numeric=false);

>

restart;

>	map(M->int(1e6exp(-1e6t), t=0.0..infinity, method=M), [_d01amc,gquad,_Dexp,_d01ajc]);

>	restart; Digits:=16:

>	map(M->int(1e6exp(-1e6t), t=0.0..infinity, method=M), [_d01amc,gquad,_Dexp,_d01ajc]);

Download int_deta_2024.2.mw

structured type...

Answered: acer 32440

September 20 2025

1 0

For example,

type(y(x), And(typefunc(symbol,symbol),
               satisfies(u->nops(u)=1)));

                true

If you actually intended name instead of (either instance of) symbol then you could change accordingly.

animations...

Answered: acer 32440

September 17 2025

0 0

The animate export of the solution modules returned by pdsolve,numeric is considerately more efficient than generating frames as separate plots and thus repeatedly pushing the computation further out.

The results of separate animations can also be ripped out of each structure and then recombined so that they show together.

The following seems reasonably quick, taking about 12 seconds to do the computational work on my machine, ie. generating the 4*2*2*35=560 curves.

I've used Maple 2018, as the OP did.

>

$OdeSys := z4*(X*(diff(Theta(X, tau), X, X))+diff(Theta(X, tau), X))/z3-(diff(Theta(X, tau), tau))-Nc*(Theta(X, tau)-Thetaa)^2*z5/z1-M^2*(Theta(X, tau)-Thetaa)^(p+1)*(1+n1*z2/z3)/(z3*(1-Thetaa)^p)-Nr*(Theta(X, tau)^4-Thetaa^4)/z3+QG*X*(1+lambda1*(Theta(X, tau)-Thetaa)+lambda2*(Theta(X, tau)-Thetaa)^2+lambda3*(Theta(X, tau)-Thetaa)^3)/z3; Cond := {Theta(0, tau) = 1+a*sin(alpha1), Theta(X, 0) = 0, (D[1](Theta))(1, tau) = 0}$

>

$OdeSys1 := z4*(X*(diff(Theta(X, tau), X, X))+diff(Theta(X, tau), X))/z3-(diff(Theta(X, tau), tau))-Nc*(Theta(X, tau)-Thetaa)^2*z5/z1-M^2*(Theta(X, tau)-Thetaa)^(p+1)*(1+n2*z2/z3)/(z3*(1-Thetaa)^p)-Nr*(Theta(X, tau)^4-Thetaa^4)/z3+QG*X*(1+lambda1*(Theta(X, tau)-Thetaa)+lambda2*(Theta(X, tau)-Thetaa)^2+lambda3*(Theta(X, tau)-Thetaa)^3)/z3; Cond1 := {Theta(0, tau) = 1+a*sin(alpha1), Theta(X, 0) = 0, (D[1](Theta))(1, tau) = 0}$

>

$OdeSysa := z4*(X*(diff(Theta(X, tau), X, X))+diff(Theta(X, tau), X))/z3-(diff(Theta(X, tau), tau))-Nc*(Theta(X, tau)-Thetaa)^2*z5/z1-M^2*(Theta(X, tau)-Thetaa)^(p+1)*(1+n1*z2/z3)/(z3*(1-Thetaa)^p)-Nr*(Theta(X, tau)^4-Thetaa^4)/z3+QG*X*(1+lambda1*(Theta(X, tau)-Thetaa)+lambda2*(Theta(X, tau)-Thetaa)^2+lambda3*(Theta(X, tau)-Thetaa)^3)/z3; Conda := {Theta(0, tau) = 1+a*cos(alpha1), Theta(X, 0) = 0, (D[1](Theta))(1, tau) = 0}$

>

$OdeSysa1 := z4*(X*(diff(Theta(X, tau), X, X))+diff(Theta(X, tau), X))/z3-(diff(Theta(X, tau), tau))-Nc*(Theta(X, tau)-Thetaa)^2*z5/z1-M^2*(Theta(X, tau)-Thetaa)^(p+1)*(1+n2*z2/z3)/(z3*(1-Thetaa)^p)-Nr*(Theta(X, tau)^4-Thetaa^4)/z3+QG*X*(1+lambda1*(Theta(X, tau)-Thetaa)+lambda2*(Theta(X, tau)-Thetaa)^2+lambda3*(Theta(X, tau)-Thetaa)^3)/z3; Conda1 := {Theta(0, tau) = 1+a*cos(alpha1), Theta(X, 0) = 0, (D[1](Theta))(1, tau) = 0}$

>

NN := 35:
str := time[real]():
for j from 1 to nops(NrVals) do
Ans := pdsolve((eval([OdeSys,Cond],Nr=NrVals[j]))[],numeric,spacestep=1/200,timestep=1/100):
Ans1 := pdsolve((eval([OdeSys1,Cond1],Nr=NrVals[j]))[],numeric,spacestep=1/200,timestep=1/100):
Ansa := pdsolve((eval([OdeSysa,Conda],Nr=NrVals[j]))[],numeric,spacestep=1/200,timestep=1/100):
Ansa1 := pdsolve((eval([OdeSysa1,Conda1], Nr=NrVals[j]))[],numeric,spacestep=1/200,timestep=1/100):

eta[j] := Ans:-animate((int(z3*z5*Nc*(Theta(X, tau)-Thetaa)^2/(z1*z4)+NrVals[j]*(Theta(X, tau)^4-Thetaa^4)/z4+M^2*(Theta(X, tau)-Thetaa)^(p+1)*(1+n1*z2/z3)/(z4*(1-Thetaa)^p), X = 0 .. 1,numeric,epsilon=1e-5,method=_Dexp))/(z3*z5*Nc*(1-Thetaa)^2/(z1*z4)+NrVals[j]*(-Thetaa^4+1)/z4+M^2*(1-Thetaa)*(1+n1*z2/z3)/z4), tau = 0 .. 2, X = 0.01..1.0, frames=NN, linestyle = "solid", 'axes' = 'boxed', labels = [" τ ", " η "], color = colour[j], title="X = %4.3f"):
eta1[j] := Ans1:-animate((int(z3*z5*Nc*(Theta(X, tau)-Thetaa)^2/(z1*z4)+NrVals[j]*(Theta(X, tau)^4-Thetaa^4)/z4+M^2*(Theta(X, tau)-Thetaa)^(p+1)*(1+n2*z2/z3)/(z4*(1-Thetaa)^p), X = 0 .. 1,numeric,epsilon=1e-5,method=_Dexp))/(z3*z5*Nc*(1-Thetaa)^2/(z1*z4)+NrVals[j]*(-Thetaa^4+1)/z4+M^2*(1-Thetaa)*(1+n2*z2/z3)/z4), tau = 0 .. 2, X = 0.01..1.0, frames=NN, linestyle = "dash", 'axes' = 'boxed', labels = [" τ ", " η "], color = colour[j], title="X = %4.3f");
etaa[j] := Ansa:-animate((int(z3*z5*Nc*(Theta(X, tau)-Thetaa)^2/(z1*z4)+NrVals[j]*(Theta(X, tau)^4- Thetaa^4)/z4+M^2*(Theta(X, tau)-Thetaa)^(p+1)*(1+n1*z2/z3)/(z4*(1-Thetaa)^p), X = 0 .. 1,numeric,epsilon=1e-5,method=_Dexp))/(z3*z5*Nc*(1-Thetaa)^2/(z1*z4)+NrVals[j]*(-Thetaa^4+1)/z4+M^2*(1-Thetaa)*(1+n1*z2/z3)/z4), tau = 0 .. 2, X = 0.01..1.0, frames=NN, linestyle = "solid", 'axes' = 'boxed', labels = [" τ ", " η "], color = colour1[j], title="X = %4.3f");
etaa1[j] := Ansa1:-animate((int(z3*z5*Nc*(Theta(X, tau)-Thetaa)^2/(z1*z4)+NrVals[j]*(Theta(X, tau)^4-Thetaa^4)/z4+M^2*(Theta(X, tau)-Thetaa)^(p+1)*(1+n2*z2/z3)/(z4*(1-Thetaa)^p), X = 0 .. 1,numeric,epsilon=1e-5,method=_Dexp))/(z3*z5*Nc*(1-Thetaa)^2/(z1*z4)+NrVals[j]*(-Thetaa^4+1)/z4+M^2*(1-Thetaa)*(1+n2*z2/z3)/z4), tau = 0 .. 2, X = 0.01..1.0, frames=NN, linestyle = "dash", 'axes' = 'boxed', labels = [" τ ", " η "], color = colour1[j], title="X = %4.3f");

Plots1[j] := Ans1:-animate(Theta(X, tau), tau = 0..2, linestyle = "dash", labels = ["Y", Theta(Y, tau)], color = colour[j], 'axes' = 'boxed', frames=NN, title="tau = %4.3f");
Plotsa[j] := Ansa:-animate(Theta(X, tau), tau = 0..2, linestyle = "solid", labels = ["Y", Theta(Y, tau)], color = colour1[j], 'axes' = 'boxed', frames=NN, title="tau = %4.3f");
Plots[j] := Ans:-animate(Theta(X, tau), tau = 0..2, linestyle = "solid", labels = ["Y", Theta(Y, tau)], color = colour[j], 'axes' = 'boxed', frames=NN, title="tau = %4.3f");
Plotsa1[j] := Ansa1:-animate(Theta(X, tau), tau = 0..2, linestyle = "dash", labels = ["Y", Theta(Y, tau)], color = colour1[j], 'axes' = 'boxed', frames=NN, title="tau = %4.3f");

end do:
(time[real]()-str)*`seconds`;

>	plots:-display(seq(plots:-display(seq([seq(PLOT(op([1,KK,..],RR[j]),op(2..,RR[j])), RR=[Plots,Plotsa,Plots1,Plotsa1])][], j=1..nops(NrVals))), KK=1..NN),insequence,size=[700,700]);

>	plots:-display(seq(plots:-display(seq([seq(PLOT(op([1,KK,..],RR[j]),op(2..,RR[j])), RR=[eta,etaa,eta1,etaa1])][], j=1..nops(NrVals))), KK=1..NN),insequence,size=[700,700]);

Download paper2_new_efficiency_plots_2025_acc.mw

ranges...

Answered: acer 32440

September 09 2025

1 0

You can even get a nice closed form for the expression that represents this curve, in terms of either variable.

>

restart;

>	f3:=-1/(((-qh-1)Tch+qc+1)deltac+(qh+1)Tch)(6alpha^2(alpha-1)(alpha+1)(((-qh-1)Tch+qc+1)^2deltac^2+((-7/3(qh^2)-10/3qh-1)Tch^2+(2(qh+1))(qc+1)Tch+(1/3)qc^2-2qc(1/3)-1)deltac+4qhTch^2(qh+1)(1/3)))=0:

>	params:=[qh=0.64,deltac=0.5,deltah=0.5,alpha=0.8,gamma=0.1,M=1.6205];

>	Eq1:=eval(collect(f3,[Tch,qc]), params);

>	plots:-implicitplot(Eq1, qc=0..1, Tch=0..1);

>	S1:=solve({(lhs-rhs)(eval(f3,convert(params,rational)))=0, qc>0, Tch>0}, [qc,Tch]);

[[qc < 3/5, 0 < qc, Tch = -15*qc-15+(10/41)*(3526*qc^2+7462*qc+3936)^(1/2)]]

>	S2:=solve({(lhs-rhs)(eval(f3,convert(params,rational)))=0, qc>0, Tch>0}, [Tch,qc]);

[[Tch < -15+(40/41)*246^(1/2), 0 < Tch, qc = -(123/125)*Tch-1/5+(2/125)*(3526*Tch^2-6150*Tch+2500)^(1/2)]]

>	indets([S1,S2],`<`); evalf(%);

${0 < Tch, 0 < qc, Tch < -15+(40/41)*246^(1/2), qc < 3/5}$

Download 1_6_ac.mw

The range qc=0.34..0.5, Tch=0.01..0.05 that you were using doesn't contain what you're after. It just misses the curve where the expression is zero...

convert example...

Answered: acer 32440

September 08 2025

0 0

You could augment the convert command.

Download conv_sinc_ex.mw

And you could also make such a convert(...,sinc) call for a compound expression involving a mix of sin&cos&tan, naturally.

ps. I changed your Post into a Question.

pps. You could also handle csc/sec/cot , in a few ways, eg. by converting to sincos as intermediate form.

ppps. It's possible to simply replace name sin/cos/tan by your operators, under eval. I deliberately used evalindets instead, to target only function calls of those names, and not any (other) instances of the bare names. Examples where is might matter could be rare, and possibly look contrived...

Explore...

Answered: acer 32440

September 08 2025

1 4

An Explore example, using your attachment g1.mw .

The signifant first change here is to remove the fixed view option from the plot3d call (since the surface's scale changes so much with the parameters) and to supply the adaptview=false option to Explore so that it allow each frame to get its own view independently. See magenta text edits.

g1_ac.mw

>

G := proc( alpha__1,alpha__2,alpha__3,alpha__4,alpha__5 )
   global last;
   last := [[:-alpha[1] = alpha__1,:-alpha[2] = alpha__2,:-alpha[3] = alpha__3,
             :-alpha[4] = alpha__4,:-alpha[5] = alpha__5],
            eval(M, [:-alpha[1] = alpha__1,:-alpha[2] = alpha__2,:-alpha[3] = alpha__3,
                     :-alpha[4] = alpha__4,:-alpha[5] = alpha__5])];
   plot3d(Im(eval(last[2])), a = -10 .. 10, l = -10 .. 10,
          #view = -10 .. 10,
          grid = [100, 100],
          colorscheme="turbo",
          style = surface, adaptmesh = false, size = [500, 500]); end proc:

>

last := 'last'; Explore(G(`α__1`, `α__2`, `α__3`, `α__4`, `α__5`), `α__1` = -5.00000001 .. 5.000000001, `α__2` = -5.00000001 .. 5.000000001, `α__3` = -5.00000001 .. 5.000000001, `α__4` = -5.00000001 .. 5.000000001, `α__5` = -5.00000001 .. 5.000000001, initialvalues = [`α__1` = 1, `α__2` = 1, `α__3` = 1, `α__4` = 1, `α__5` = 5], placement = right, adaptview = false)

example...

Answered: acer 32440

September 07 2025

2 1

The following tests that the first indices (yours have just one) of the elements are four distinct values.

Download map_op_ex.mw

You haven't stated explicitly that all the lists in L will always have elements that are indexed by at least one value, but I have presumed it so.

ps. I didn't see dharr's Answer before posting my own. The fact that they are so similar is really because it's a natural way to approach it.
[edit] dharr's Answer seems to have gone missing. But it was fine, IMO, (presuming all elements have exactly one index, which is also reasonable in the absence of stated details),
select(x->nops({map(op,x)[]})=4,L);

[edit] Some explanation of the steps, so that next time you might build your own,

>	Turn the inner lists of indexed items into lists of just their indices. map( LL -> map2(op,1,LL), L );

>	But now also make them sets of the indices, to get rid of duplicates. map( LL -> {map2(op,1,LL)[]}, L );

>	But now take the nops of those sets, and set equal to 4. map( LL -> nops({map2(op,1,LL)[]})=4, L );

>	Now use those equations as boolean valued tests, in order to select with this action. select(LL->nops({map2(op,1,LL)[]})=4,L);

>	Similarly, but without using select. map( LL -> ifelse(nops({map2(op,1,LL)[]})=4, LL, NULL), L );

some ideas...

Answered: acer 32440

August 29 2025

1 4

I'm guessing that you'd prefer to have the Product call be shown with the assigned name, rather than with the Vector or list to which that name is assigned.

Here is something that gets that, for Vector or list, in Maple 2021.

>

restart;

>	kernelopts(version);

>	with(Statistics):

>	Dist := proc(N::{posint, name}, alpha::evaln) local K; uses Statistics; K := numelems(eval(alpha)); Distribution( ':-PDF' = (z -> Product(('alpha'[k]+z[k])^N, k=1..K)) ) end proc:

>	P := Vector(3, [13,2,4]);

Vector(3, {(1) = 13, (2) = 2, (3) = 4})

>	X := RandomVariable(Dist(N, P)):

>	res1 := PDF(X, x);

Product((P[k]+x[k])^N, k = 1 .. 3)

>	eval(res1,1);

Product((P[k]+x[k])^N, k = 1 .. 3)

>	value(eval(res1,1));

>	L := [3,17,22]; X2 := RandomVariable(Dist(N, L)):

>	res2 := PDF(X2, y);

Product((L[k]+y[k])^N, k = 1 .. 3)

>	value(eval(res2,1)); value(res2);

>	eval(res2,1); res2;

Product((L[k]+y[k])^N, k = 1 .. 3)

Product(([3, 17, 22][k]+y[k])^N, k = 1 .. 3)

>	V := Vector([1/2,1/3,1/4]): subs('P'='V', eval(res1,1)); value(eval(%,1));

Product((V[k]+x[k])^N, k = 1 .. 3)

Download Product_error_ac.mw

and a variant on that: Product_error_ac2.mw

Here is something that behaves ok for a list (but not Vector), in Maple 2015: Product_error_ac_2015_list.mw

one way...

Answered: acer 32440

August 26 2025

4 0

With a little effort one can guard against collision of the textplots.

Download Animation_Trigo_ac.mw

examples...

Answered: acer 32440

August 26 2025

0 0

Your curves don't change much for such small changes in E1, it seems.

AGM_method_single_line_plot_ac.mw

There's quite a bit of duplicated effort in your approaches (which I didn't try to improve, though I showed a couple or alternative syntax blobs).

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

32440 Reputation

29 Badges

Social Networks and Content at Maplesoft.com

MaplePrimes Activity

These are answers submitted by acer

Physics, PDEtools:-undeclare...

collect, normal...

Scale...

sequence of two expressions...

ScrollableMathTableOutput = false...

the changed behavior...

structured type...

animations...

ranges...

convert example...

Explore...

example...

some ideas...

one way...

examples...

Save this setting as your default sorting preference?

Ask a Question

Create a Post

Generating PDF…

Save this setting as your default sorting preference?
Note: You can change your preference any time in your account settings
Don't show this again