Hello,

I have this simple problem which doesn't want to work.

restart;
d := g -> (D@@2)(g) - x^2*g;
((d@@2)(g) assuming x::constant);

The result of the last line is:

(D@@4)(g) - 2*(D@@2)(x)*x*g - 2*D(x)^2*g - 4*D(x)*x*D(g) - x^2*(D@@2)(g) - x^2*((D@@2)(g) - x^2*g)

so Maple doesn't set D(x) to 0. On the other hand if I just write

D(x) assuming x::constant

then Mapel returns 0.

Similarly

D(f^k) assuming k::constant

just returns D(f^k) and not k*D(f)*f^(k-1) as the example in HELP suggests.

If I select Insert | Header Footer ... and choose the Custom Header tab, I don't see a way to create a multi-line header. How do I create a custom multi-line header?

restart;
with(Optimization);
with(Student[Calculus1]);
z1 := 0.15;
                           z1 := 0.15

z2 := 0.85;
                           z2 := 0.85

NULL;
Q := Matrix(3, 3, [1, 0, 0, 1, 1, 0, 1, 3, 3]);
P_S := Matrix(1, 3, [ps0, ps1, ps2]);
Phi_S := Matrix(3, 1, [1, t^(1 + gs1), t^(2 + gs2)]);
dz1Phi_S := t^z1*Matrix(3, 3, [0, 0, 0, 0, Gamma(gs1 + 2)/Gamma(gs1 + 2 - z1), 0, 0, 0, Gamma(gs2 + 3)/Gamma(gs2 + 3 - z1)]);
dz2Phi_S := t^z2*Matrix(3, 3, [0, 0, 0, 0, Gamma(gs1 + 2)/Gamma(gs1 + 2 - z2), 0, 0, 0, Gamma(gs2 + 3)/Gamma(gs2 + 3 - z2)]);
NULL;
P_I := Matrix(1, 3, [p_i0, p_i1, p_i2]);
Phi_I := Matrix(3, 1, [1, t^(1 + gi1), t^(2 + gi2)]);
NULL;
dz1Phi_I := t^z1*Matrix(3, 3, [0, 0, 0, 0, Gamma(gi1 + 2)/Gamma(gi1 + 2 - z1), 0, 0, 0, Gamma(gi2 + 3)/Gamma(gi2 + 3 - z1)]);
dz2Phi_I := t^z2*Matrix(3, 3, [0, 0, 0, 0, Gamma(gi1 + 2)/Gamma(gi1 + 2 - z2), 0, 0, 0, Gamma(gi2 + 3)/Gamma(gi2 + 3 - z2)]);
P_H := Matrix(1, 3, [ph0, ph1, ph2]);
Phi_H := Matrix(3, 1, [1, t^(1 + gh1), t^(2 + gh2)]);
dz1Phi_H := t^z1*Matrix(3, 3, [0, 0, 0, 0, Gamma(gh1 + 2)/Gamma(gh1 + 2 - z1), 0, 0, 0, Gamma(gh2 + 3)/Gamma(gh2 + 3 - z1)]);
dz2Phi_H := t^z1*Matrix(3, 3, [0, 0, 0, 0, Gamma(gh1 + 2)/Gamma(gh1 + 2 - z2), 0, 0, 0, Gamma(gh2 + 3)/Gamma(gh2 + 3 - z2)]);
P_L := Matrix(1, 3, [pl0, pl1, pl2]);
Phi_L := Matrix(3, 1, [1, t^(1 + gl1), t^(2 + gl2)]);
dz1Phi_L := t^z1*Matrix(3, 3, [0, 0, 0, 0, Gamma(gl1 + 2)/Gamma(gl1 + 2 - z1), 0, 0, 0, Gamma(gl2 + 3)/Gamma(gl2 + 3 - z1)]);
dz2Phi_L := t^z1*Matrix(3, 3, [0, 0, 0, 0, Gamma(gl1 + 2)/Gamma(gl1 + 2 - z2), 0, 0, 0, Gamma(gl2 + 3)/Gamma(gl2 + 3 - z2)]);
S := (P_S . Q) . Phi_S;
NULL;
ds1 := simplify(((P_S . Q) . dz1Phi_S) . Phi_S, GAMMA);
ds2 := simplify(((P_S . Q) . dz2Phi_S) . Phi_S);
H := (P_H . Q) . Phi_H;
NULL;
NULL;
dh1 := simplify(((P_H . Q) . dz1Phi_H) . Phi_H, GAMMA);
dh2 := simplify(((P_H . Q) . dz2Phi_H) . Phi_H);
NULL;
L := (P_L . Q) . Phi_L;
NULL;
dl1 := simplify(((P_L . Q) . dz1Phi_L) . Phi_L, GAMMA);
dl2 := simplify(((P_L . Q) . dz2Phi_L) . Phi_L);
NULL;
I3 := (P_I . Q) . Phi_I;
di1 := simplify(((P_I . Q) . dz1Phi_I) . Phi_I, GAMMA);
di2 := simplify(((P_I . Q) . dz2Phi_I) . Phi_I);
RS1 := ds1 + (-0.0043217^z1 + ((0.125^z1*S) . I3) + (0.002^z1 + 0.0008^z1)*S);
RS2 := ds2 + (-0.0043217^z2 + ((0.125^z2*S) . I3) + (0.002^z2 + 0.0008^z2)*S);
RH1 := dh1 + (-0.535^z1 + ((0.0056^z1*H) . I3) - 0.35^z1 + (0.002^z1 + 0.0008^z1)*H);
RH2 := dh2 + (-0.535^z2 + ((0.0056^z2*H) . I3) - 0.35^z2 + (0.002^z2 + 0.0008^z2)*H);
RI1 := di1 + (-((0.125^z1*S) . I3) - ((0.0056^z1*H) . I3) - 0.029^z1*L + (0.002^z1 + 0.0008^z1 + 0.025^z1 + 0.35^z1)*I3);
RI2 := di2 + (-((0.125^z2*S) . I3) - ((0.0056^z2*H) . I3) - 0.029^z2*L + (0.002^z2 + 0.0008^z2 + 0.025^z2 + 0.35^z2)*I3);
RL1 := dl1 + (-0.025^z1*I3 + (0.002^z1 + 0.0008^z1 + 0.029^z1)*L);
RL2 := dl2 + (-0.025^z2*I3 + (0.002^z2 + 0.0008^z2 + 0.029^z2)*L);
R15 := evalf(RH1^2 + RI1^2 + RL1^2 + RS1^2, 4);

I1 := evalf(ApproximateInt(R15[1, 1], t = 0 .. 1), 4);
NULL;
h := 0.01;
B := zeros([1, 20]);
A := Matrix(1, 20, [ps0, ps1, ps2, p_i0, p_i1, p_i2, ph0, ph1, ph2, pl0, pl1, pl2, gs1, gs2, gi1, gi2, gh1, gh2, gl1, gl2]);
NULL;
for i to 20 do
    B[1, i] := evalf(diff(I1, A[1, i]), 4);
end do;
NULL;
NLPSolve(I1, {B[1, 1] = 0, B[1, 2] = 0, B[1, 3] = 0, B[1, 4] = 0, B[1, 5] = 0, B[1, 6] = 0, B[1, 7] = 0, B[1, 8] = 0, B[1, 9] = 0, B[1, 10] = 0, B[1, 11] = 0, B[1, 12] = 0, B[1, 13] = 0, B[1, 14] = 0, B[1, 15] = 0, B[1, 16] = 0, B[1, 17] = 0, B[1, 18] = 0, B[1, 19] = 0, B[1, 20] = 0});
Error, (in Optimization:-NLPSolve) could not store .2000*(4.61629127484374902+.6380*(.100000000000000019e-1*Gamma(3.)*Gamma(3.85)+.375000000000000116e-3*Gamma(4.)*Gamma(2.85))/Gamma(2.85)/Gamma(3.85))*(.6380*(.100000000000000019e-1*Gamma(3.)*D(Gamma)(3.85)+.375000000000000116e-3*D(Gamma)(4.)*Gamma(2.85)-.112350000000000015e-2*Gamma(4.)*Gamma(2.85))/Gamma(2.85)/Gamma(3.85)-.6380*(.100000000000000019e-1*Gamma(3.)*Gamma(3.85)+.375000000000000116e-3*Gamma(4.)*Gamma(2.85))/Gamma(2.85)/Gamma(3.85)^2*D(Gamma)(3.85)-.238097850000000034e-2)+.594623719800000056e-3*(.333300000000000046e-2*Gamma(3.)*Gamma(3.85)+.125000000000000030e-3*Gamma(4.)*Gamma(2.8...

Hello everybody. after running this code in Maple i get NLPSolve error (could not store...). I asked this question before and shared a picture. now i well share the code. how can i fix this error?

Equilibrium._solutions.mw

I need to obtain the equilibrium solutions when infection is present as well as analyze the local asymptotic stability when infection is absent and present in the system

I want to partially execute a large worksheet using to a point where I want to proceed with step-by-step execution . I am looking for a command that brings execution to a halt.

Threads:-Sleep(big number) and works but is there anything else that I could use.

In the same context: Execute worksheet removes all output. I would like to keep the output in the worksheet. Is that possible?

Hi,

I am building some procedures to animate trigonometric functions, but the procedure H with the tangent function does not produce the expected result. Maybe I have overlooked a detail? Thank you for your insights.AnimationCercleTrigoFonctiTrigo.mw

Hi, When I associate the function j with a random expression, its derivative does not follow. I imagine that I have forgotten something in my initial assignment. Any hints? Thanks

RandomfunctionQ.mw

I have created a several matrices for my analysis. I want to plot the determinant value of matrix vs the omega. The matrix that should be plotted is named as FINAL. I am not able to do that due to large value of matrix determinant. How to simplify the matrix and get  plot.

Download Plot_determinant_value_vs_omega.mw

Today, I trried to enter an equation into a text region using Ctrl+R (In Maple Flow 2023.1), like I have always done in previous releases of Maple Flow. Well, today it did not work. I tried to enter the equation Q=W+mCv(T2-T1). When I tried to enter the equal sign after the Q, the program would not allow me to do it. I would press the equal sign and nothing happened. I tried :=, but it would only enter the :, it would not let me enter the =. Any help will be appreciated.

Hi,

So I've just installed the 2022 version of Maple and I wanted to make a new document.

But when I'm trying to make a variable with a subscript that includes a comma, I get wrong output.

For example; I want to make the variable A_b,c as 2D input and then it shows 'A_b,c' as output.

Why are those apostrophs showing up?

Thanks for the help

I often use DocumentTools:-Tabulate or  DocumentTools:-Layout do display a vector or matrix of plots instead of plots:-display. because I find the latter less practical.
But it seems that the 'background' option is not correctly managed with DocumentTools:-Tabulate or  DocumentTools:-Layout.
The attached file shows that:

• the 'background' option is correctly managed if each "view" contains a single plot,
• but not correctly as soon as at least one "view" contains morethan one plot.

DocumentTools_and_Background.mw

How can I get with DocumentTools:-Tabulate  /  DocumentTools:-Layout the same rendering I get with plots:-display?

Hi,

I'm looking for a harmonized layout for my weighted probability trees. Are there any possibilities to improve this code?

Thanks

Arbrepondéré.mw

Good day everyone, 

How can I extract the values of x and y for plotting? 

The worksheet is attached below. Thanks

dont_get_it.mw

I couldn't locate a command to change the default frames per second (fps) of 10.  Is there an option?

The solve function is able to solve inequalities and return solutions that satisfy the inequalities. Unfortunately, it appears that for some toy problems, this command does not always perform efficiently: 
 

Download x=y=z=3.mw

So none of them is capable of giving a complete description of the region defined by or  with a time limit of "ten seconds". (I believe that completing the evaluation may take a much long time; actually, if there is no constraint on the real time allowed, the symbolic solver will be stuck.) (Note that here it is not hard to find individual instances of values of variables that satisfy them (for example, via the 3-D region boundary plot above), yet they provide less insight into all solutions that exist.)
But in my view,  are too long for these two special systems.

Can Maple solve them more quickly at present? 

> eqns__1 := {x^2+y^2+z^2+x*y*z <= 2*(y*z+z*x+x*y), 2*(x^2+y^2+z^2) < y^2*z+z^2*x+x^2*y-27, x >= 0, y >= 0, z >= 0};
> eqns__2 := {x^2+y^2+z^2+x*y*z <= 2*(y*z+z*x+x*y), 2*(x^2+y^2+z^2) <= y^2*z+z^2*x+x^2*y-27, x > 0, y > 0, z > 0};

(It seems that a potential difficult part for Maple is likely attributable to proving completeness of the solution ….)