Maple 2022 Questions and Posts

These are Posts and Questions associated with the product, Maple 2022

I am learning very basic matrix ops in Maple and I'm running into questions nearly at every step.

At the risk asking too many trivial questions, I have written my ~7 questions in the attached worksheet in red font.

Any answers will be appreciated.

matrix_ops_Qs.mw

matrix_ops_Qs.pdf

Hello

After using maple for quite a while, I am still confused by some basic concepts.   

Given the following table, how can I select the indices for which the entry is not null ([])?  

table([1 = NULL, 2 = NULL, 3 = NULL, 4 = NULL, 5 = NULL, 6 = NULL, 7 = 5, 9 = NULL, 8 = NULL, 11 = 4, 10 = NULL, 13 = NULL, 12 = NULL, 15 = 9, 14 = NULL, 18 = NULL, 19 = 8, 16 = 9, 17 = NULL, 22 = 9, 23 = NULL, 20 = NULL, 21 = 8, 27 = NULL, 26 = 8, 25 = 4, 24 = NULL, 31 = NULL, 30 = 9, 29 = NULL, 28 = 9, 36 = NULL, 37 = 9, 38 = 9, 39 = NULL, 32 = 5, 33 = NULL, 34 = NULL, 35 = NULL, 45 = NULL, 44 = NULL, 47 = NULL, 46 = NULL, 41 = 8, 40 = NULL, 43 = NULL, 42 = NULL, 54 = NULL, 55 = NULL, 52 = NULL, 53 = NULL, 50 = NULL, 51 = NULL, 48 = 5, 49 = 9, 60 = 8, 59 = NULL, 58 = 7, 57 = 7, 56 = NULL])

Many thanks

Hi, 

Maple has corrupted my file without notice. I get the error "There was a problem in the loading process, you worksheet may be incomplete."

I've tried restarting maple, checked my backup folder, and followed maples troubleshooting process (maplesoft support) without any luck.

Anyone who might help me here? 

File: Assignment2.mw

The algebraic equation in four variables is: 

restart;
expr := (a^2 + b^2 + c^2 + d^2 - 4)*(a^2 + b^2 + c^2 + d^2) + 4*d*c*b*a - (a + b + c + d - 2)*(a + b + c + d) + 6:
solve(expr = 0, useassumptions, allsolutions) assuming nonnegative;

Unfortunately, I've been waiting for a long time, and I have no more time to wait; I have to interrupt the current session by hand. However, it appears that MMA can solve it within bearable time:

Did Maple miss something here?

in my program, I keep assumptions in a set. Sometimes this is empty if no assumptions are used. This never caused a problem before (at least I do not think so, else I would have seen it) when using empty {} in assuming, except for now.

Here is one example below. Is this a known problem? I noticed when changing {} to [] the error goes away. I am not sure why, and if this is known issue. But will change from a set to a list to avoid this. 

Maple 2022.2 on windows 10

interface(version);

`Standard Worksheet Interface, Maple 2022.2, Windows 10, October 23 2022 Build ID 1657361`

Physics:-Version();

`The "Physics Updates" version in the MapleCloud is 1392 and is the same as the version installed in this computer, created 2023, February 13, 12:58 hours Pacific Time.`

restart;

ode:=diff(diff(y(x),x),x)+diff(y(x),x)^2+diff(y(x),x) = 0;
ic:=y(0) = 0;
sol:=y(x) = -ln(exp(x))+ln(-1+_C1*exp(x))-ln(-1+_C1);

diff(diff(y(x), x), x)+(diff(y(x), x))^2+diff(y(x), x) = 0

y(0) = 0

y(x) = -ln(exp(x))+ln(-1+_C1*exp(x))-ln(-1+_C1)

odetest(sol,[ode,ic]) assuming {};

Error, (in convert/multiset) too many levels of recursion

odetest(sol,[ode,ic]) assuming [];

[0, 0]

odetest(sol,[ode, ic]);

[0, 0]

 

Download feb_14_2023.mw

I am new to Maple. I'm having trouble with simple caluclations such as the one in the image. The answer (after simplificatoin and cancellations) should be (|a|^2+|b|^2)(x^2 + y^2 + z^2)  but it's not happening! The result shown is correct; only the simplification is missing. What step am I missing?

The calculation is so simple, a look at the picture should suffice for the expert,

(sorry! I still don't know how to pick-off a single tab out from the whole workbook).

Question-2: judging from the missing overbar (complex conj) on z, the "assuming" seems to only apply to z and not x and y.

if that is correct, then why?

Question-3: How can I implement the constraint (|a|^2+|b|^2)=1 (preferably)before or after this  multiplication? I tried assume((|a|^2+|b|^2)=1), and though it compiled it did not help with the simplification.

TYVM

i have on ode with some parameters. it is ok and it is solved. but when i substitute the parameter itself, it is not solved, can i use any assumption to solve this? tnx for the help.

restart

ode := diff(T(x), x, x)+q/k = 0

diff(diff(T(x), x), x)+q/k = 0

(1)

dsolve(ode)

T(x) = -(1/2)*q*x^2/k+_C1*x+_C2

(2)

ics1 := -k*(D(T))(0) = h[1]*(T[inf1]-T(0)), -k*(D(T))(0.5e-1) = h[2]*(T(0.5e-1)-T[inf2])

-k*(D(T))(0) = h[1]*(T[inf1]-T(0)), -k*(D(T))(0.5e-1) = h[2]*(T(0.5e-1)-T[inf2])

(3)

dsolve({ics1, ode})

T(x) = -(1/2)*q*x^2/k-(1/40)*h[1]*(800*k*T[inf1]*h[2]-800*k*T[inf2]*h[2]-40*k*q-q*h[2])*x/((20*k*h[1]+20*k*h[2]+h[1]*h[2])*k)+(1/40)*(800*k*T[inf1]*h[1]+800*k*T[inf2]*h[2]+40*T[inf1]*h[1]*h[2]+40*k*q+q*h[2])/(20*k*h[1]+20*k*h[2]+h[1]*h[2])

(4)

ics2 := -k*(D(T))(0) = h[1]*(T[inf1]-T(0)), -k*(D(T))(L) = h[2]*(T(L)-T[inf2])

-k*(D(T))(0) = h[1]*(T[inf1]-T(0)), -k*(D(T))(L) = h[2]*(T(L)-T[inf2])

(5)

dsolve({ics2, ode})

Error, (in dsolve) found differentiated functions with same name but depending on different arguments in the given DE system: {T(L), T(x)}

 

``

Download ExactSol.mw

My question is can the last step be equal to 1/6.

I want to output to 1/6. it outputs to sqrt(9)/18 - (See Below).

NULL"5)Square root: undefined - DNE (Limit does not exist):"

NULL

NULL

NULL

limit((sqrt(x+1)-3)/(x-8), x = 8)"(=)"1/6 

 

 

limit((sqrt(x+1)-3)/(x-8), x = 8)Limit(((x+1)^(1/2)-3)/(x-8), x = 8) = (1/18)*9^(1/2)"(=)"Limit(((x+1)^(1/2)-3)/(x-8), x = 8) = 1/6

NULL

Download limit-sqrt-5.mw

What is the logic behind this?

expr:=(c[2]+x)^3;
simplify(expr)

Which is what expected. But 

expr:=(c[2]+x)^3+a;
simplify(expr)

gives

Which does not look simpler to me. I expected it to be the same as before but with "a" added.  This is what Mathematica gives for comparison

I know I can use simplify with size option. But my question is, how did Maple decide that x^3 + 3*x^2*c[2] + 3*x*c[2]^2 + c[2]^3 + a is "simpler" than (c[2] + x)^3 + a ? It must use some logic which I am trying to understand.

Maple 2022.2 on windows 10

 The following code:

for i to np do
    pl_pp[i] := point([h_p[i], m_p[i]], color = red);
    pl_np[i] := point([h_n[i], m_n[i]], color = blue);
    pl_pn[i] := point([-h_n[i], -m_n[i]], color = red);
    pl_nn[i] := point([-h_p[i], -m_p[i]], color = blue);
end do

Gives following error:

Error, (in plottools:-point) incorrect arguments for creating points structure, try providing the dimension option
 

Dear Maple community,

I am trying to define an operator in Maple that multiplies to an expression and substitutes the indices.

The issue is that I don't want to add the expression as input of the operator/function, sth like:

p := proc(i, j, z) option operator; subs({i = j, j = i}, z); end proc

I want it to act on expression, like:

P(i,j)*X(i,j,k)=X(j,i,k)

and X can be any tensor or matrix.

Thanks in advance.

The shortcut key Ctrl+Delete is used to clear a math region or delete a cell within a group, so I can delete some output manually. However, I find that there exists some long-standing functions like PrintTemporary and Monitor in Mathematica, which can be used to print intermediate values of variables temporarily (and dynamically, if possible). The term "temporarily" means that the printed output disappears at the end of a computation / when a computation is complete / when the evaluation finishes / before returning the result of an evaluation. For instance: 

(Note that certain operations require a front end.)

Is it possible to implement such functionalities in Maple? 

IndependencePolynomial returns the independence polynomial for the graph G in the variable x.

For the following example, its calculation took over 20 minutes and still hasn't produced a result, and what's fatal is that it  has consumed 4G  memory.

with(GraphTheory):
G:=ConvertGraph("W|tNHEpCKoh`@@Po_WHB@CKC?WGO{G?KKCB`?OMG?_y_?Sn");
G1:=LineGraph(G);
IndependencePolynomial(G1, x) # be careful

I use  codes in  the link https://github.com/pernici/hobj.

It produced results quickly (It takes approximately 5 seconds.). So I think the built-in function " IndependencePolynomial " should be able to be improved. (Of course we are usually very concerned about their coefficients) 

Their coefficients of the independent polynomial of G1 are as follows.

[340649, 12329124, 68797662, 140606548, 139481127, 77027880, 25546428, 5303544, 700911, 58580, 2982, 84, 1]

It tells me the total number of independence sets with size 12 is 340649. 

I have been unable to add a directed light source to my plot3d drawings.  Perhaps I am misreading the documentation.  In the following experiment, I attempt to shine a red light from the direction (phi,theta) but changing the values of phi and theta don't seem to change the scene's lighting.  Can you change the lighting in your Maple?

restart;
with(plottools):
with(plots):
Explore(display(sphere(), style=surface, color=yellow, light=[phi,theta,1,0,0]),
	phi=0..Pi, theta=0..Pi);

First 19 20 21 22 23 24 25 Last Page 21 of 43