[[1,1,["ab","xy"]], [0,1,["xy"]]]
 

how to sort by first number , second number , and list of string ?

how to count and group by third list of strings having first two number 0,0 or 0,1 or 1,0 or 1,1 

for example 

["xy"] at most count 4  , at least count 1

would like to find which is [1,1,... ] and only count 1

Is Prüfer's algorithm available in Maple (given a tree, produce the code, and vice versa)? If not, has anyone written the code for it?

Hello

I need to check if the solution of a polynomial system (for instance a set of polynomial equations in y and z) using two different approaches is the same (equal or symmetric).  I thought if I use simplify plus abs I could solve the problem, but that is not the case.   Here is an example;

The first method returns the following solution:

aa := {{y = -2*X1*X2*alpha[1, 8]*alpha[2, 6]/(sqrt((X1^4*alpha[1, 8]^2*alpha[2, 4]^2 + 2*X1*((-2*X1*X2*alpha[3, 6] - 2*X2*alpha[2, 2] + 2*X3)*alpha[2, 6] + X1*X2*alpha[2, 4]*(alpha[2, 8] + alpha[3, 9]))*alpha[1, 8] + X2^2*(alpha[2, 8] + alpha[3, 9])^2)*alpha[1, 8]^2) + alpha[1, 8]^2*alpha[2, 4]*X1^2 + X2*alpha[1, 8]*(alpha[2, 8] + alpha[3, 9])), z = (-sqrt((X1^4*alpha[1, 8]^2*alpha[2, 4]^2 + 2*X1*((-2*X1*X2*alpha[3, 6] - 2*X2*alpha[2, 2] + 2*X3)*alpha[2, 6] + X1*X2*alpha[2, 4]*(alpha[2, 8] + alpha[3, 9]))*alpha[1, 8] + X2^2*(alpha[2, 8] + alpha[3, 9])^2)*alpha[1, 8]^2) - alpha[1, 8]^2*alpha[2, 4]*X1^2 + (-alpha[3, 9] - alpha[2, 8])*X2*alpha[1, 8])/(2*alpha[1, 8]^2*alpha[2, 6]*X1)}, {y = -2*X2*alpha[1, 8]*alpha[2, 6]*X1/(-sqrt((X1^4*alpha[1, 8]^2*alpha[2, 4]^2 + 2*X1*((-2*X1*X2*alpha[3, 6] - 2*X2*alpha[2, 2] + 2*X3)*alpha[2, 6] + X1*X2*alpha[2, 4]*(alpha[2, 8] + alpha[3, 9]))*alpha[1, 8] + X2^2*(alpha[2, 8] + alpha[3, 9])^2)*alpha[1, 8]^2) + alpha[1, 8]^2*alpha[2, 4]*X1^2 + X2*alpha[1, 8]*(alpha[2, 8] + alpha[3, 9])), z = (sqrt((X1^4*alpha[1, 8]^2*alpha[2, 4]^2 + 2*X1*((-2*X1*X2*alpha[3, 6] - 2*X2*alpha[2, 2] + 2*X3)*alpha[2, 6] + X1*X2*alpha[2, 4]*(alpha[2, 8] + alpha[3, 9]))*alpha[1, 8] + X2^2*(alpha[2, 8] + alpha[3, 9])^2)*alpha[1, 8]^2) - alpha[1, 8]^2*alpha[2, 4]*X1^2 + (-alpha[3, 9] - alpha[2, 8])*X2*alpha[1, 8])/(2*alpha[1, 8]^2*alpha[2, 6]*X1)}}

and the second Method:

bb := {{y = -2*X2*alpha[2, 6]*X1/(-sqrt(X1^4*alpha[1, 8]^2*alpha[2, 4]^2 + 2*X2*alpha[1, 8]*((alpha[2, 8] + alpha[3, 9])*alpha[2, 4] - 2*alpha[2, 6]*alpha[3, 6])*X1^2 - 4*alpha[1, 8]*alpha[2, 6]*(X2*alpha[2, 2] - X3)*X1 + X2^2*(alpha[2, 8] + alpha[3, 9])^2) + (alpha[2, 8] + alpha[3, 9])*X2 + alpha[1, 8]*alpha[2, 4]*X1^2), z = (sqrt(X1^4*alpha[1, 8]^2*alpha[2, 4]^2 + 2*X2*alpha[1, 8]*((alpha[2, 8] + alpha[3, 9])*alpha[2, 4] - 2*alpha[2, 6]*alpha[3, 6])*X1^2 - 4*alpha[1, 8]*alpha[2, 6]*(X2*alpha[2, 2] - X3)*X1 + X2^2*(alpha[2, 8] + alpha[3, 9])^2) - alpha[1, 8]*alpha[2, 4]*X1^2 + (-alpha[3, 9] - alpha[2, 8])*X2)/(2*alpha[1, 8]*alpha[2, 6]*X1)}, {y = -2*X2*alpha[2, 6]*X1/(sqrt(X1^4*alpha[1, 8]^2*alpha[2, 4]^2 + 2*X2*alpha[1, 8]*((alpha[2, 8] + alpha[3, 9])*alpha[2, 4] - 2*alpha[2, 6]*alpha[3, 6])*X1^2 - 4*alpha[1, 8]*alpha[2, 6]*(X2*alpha[2, 2] - X3)*X1 + X2^2*(alpha[2, 8] + alpha[3, 9])^2) + alpha[1, 8]*alpha[2, 4]*X1^2 + (alpha[2, 8] + alpha[3, 9])*X2), z = (-sqrt(X1^4*alpha[1, 8]^2*alpha[2, 4]^2 + 2*X2*alpha[1, 8]*((alpha[2, 8] + alpha[3, 9])*alpha[2, 4] - 2*alpha[2, 6]*alpha[3, 6])*X1^2 - 4*alpha[1, 8]*alpha[2, 6]*(X2*alpha[2, 2] - X3)*X1 + X2^2*(alpha[2, 8] + alpha[3, 9])^2) - alpha[1, 8]*alpha[2, 4]*X1^2 + (-alpha[3, 9] - alpha[2, 8])*X2)/(2*alpha[1, 8]*alpha[2, 6]*X1)}}

Notice (if I am not mistaken) that the first pair of the first solution is equal to the second pair of the second solution.   If I compare them using evalb(simplify(aa[1,1])=simplify(bb[2,1])), Maple returns false.  Again, if I am not mistaken I think they are the same.

a) How can the solutions be compared?

b) I also need to determine if there are symmetric roots in a set of solutions (either in aa or in bb) and a procedure that returns just one solution.  Something like:

func:=(auxsolsx,varsx)->`if`(nops(map(v->op(map(w->abs(subs(w[ListTools:-Search(v,varsx)],v)),auxsolsx)),varsx))=2,ifelse(nops(auxsolsx)=1,auxsolsx,{auxsolsx[1]}),NULL):

Many thanks

Hi,

I try to represent the successive derivatives of a function using a tabulation . How do I insert titles in the seq command options?

Thanks

QuestionTabulation.mw

Apparently when raising an error exception, Maple stops further code.

Is there any possibility to just show a popup message, and then continue with the code after the message?

ErrorGoOn.mw

hello i startet using subs. and i got a problem, my expression is only finished and valid if it stands in the factorised version. i am calculating some koefficients and then with subs i am trying to insert it back into the original expression. but it always expand and simplyfy. looking at it, i can see it is a oneway arrow => after the simplicifaction, so it is impossible to rearrange.

what i would like it to happen is the expression i am asking subs to insert a koefficient into is insertet only, and not evaluated and rearrange. can i use a tecknique er something to stop the evaluation and only perform an insertion.

my equation are

a*(2*A + 2*B) + b*(C) = c*(2*A + B) + d*(B*2) + e*(C + 4*B)

i have calculated the koefficients as

{a = 6, b = 1, c = 6, d = 1, e = 1}

i would like it to look like this when i use subs

6*(2*A + 2*B) + 1(C) = 6*(2*A + B) + 1*(B*2) + 1(C + 4*B)      

but it looks like this

12*A + 12*B + C = 12*A + 12*B + C

How to determine populations y for which y ' = 0? 

Have no idea about solving Newton's law of cooling in maple. 

Is anyone can help me? Thanks! 

Hello!

I have a function, say f(x,y), and I want to solve for x for a range of values of y. The function is complicated, and thus, I gather I need to use fsolve to find the set of solutions numerically. I know how to use fsolve to find x for a given (single) value of y, but I have never used it to generate a list of x values given a range of y values. Then, I would l like to plot x as function of y, but I will leave that for later.

Thank you!

Hello!
Is there a way i can plot these two functions in one graph?

Suppose we have a 200-litre tank. Suppose 10 litres per minute of salt water (brine) flows into the tank at the top, it mixes perfectly with the contents of the tank, and 10 litres per minute of the mixture flows out of the tank at the bottom. Assume for simplicity that the salt water in the tank is stirred so that its concentration is uniform in the tank. Let S(t) be the amount of salt, in grams, in the tank at t time minutes. Suppose the salt water flowing into the tank has concentration 80 grams per litre.

Find the differential equation to model the change in S(t).

Assuming that at t = 0 time the concentration of salt in the tank is 10 grams per litre, solve the differential equation using Maple.

What is S(infinite)? That is, how much salt is in the tank after a long time?

Now graph S(t) for a suitable domain so we can see the function approaching the value S(infinite)

---------------------

Have some difficulty to understand the problem. I mean, such modelling problem rather than pure math computation.

Thank you so much!

I have a set of curves  :   plot({seq((6*x-2*t)/x^2, t = 1 .. 3)}, x = -1 .. 5, y = -1 .. 6)

and   a function           :   plot(3/x, x = 0 .. 5, y = -1 .. 6).

How is the calling sequence to plot the set and the function in the same graph ?

Is there a way to get a feedback from Maple, if it doesn't find any corresponding proc: to a procedure call?

Right now it takes very much time to bugtrack issues where there might be a spelling error in a call.

See example.

ProcNameWrong.mw

Hello, I'm new to Maple, transitioning over from Mathematica.  I have a system of equations I'm trying to solve, but am trying to substitute values for some of the variables.  I don't think I'm doing it right, because I don't see the value in place of the variable when it echos the expression. Here are my equations, and how I'm trying to solve them:

wres := sqrt(C1+C2)/sqrt)C1*C2*L1)

vores := -C1/C2*vin

z3 = (C1^2 *R1- j*sqrt(C1*C2*L1*(C1+C2)))/C2^2

sys :={wres=13.56e6*2*pi, C1=C2, z3=40,abs(vores/vin)>=1}

solve(sys, {C1, C2, L1, R1})

(I apologize; I coulnd't figure out how to enter the Maple code here, directly).

When I try this, solve returns nothing.  I then try forcing R1 to be 2000, for ex.  I tried modifying sys to be

sys:={wres=13.56e6*2*pi, C1=C2, R1=2000,  z3=40,abs(vores/vin)>=1}, but this does't work.

then I tried, R1:=2000, exepcting to see R1 substituted with 2000 in the echo from sys, but it stills shows R1....

How do I do this?  

Also, how do I go about investigating why solve is having trouble solving a system of equations?

Thanks

Let's say you make a tabulate of plots (using tabulate, because you don't want equal space for your plots as in array-plot etc). Then how should one save or export the final result as one eps or pdf file? When I right click and choose export, it only save one cell of the tabulate alone.