MaplePrimes Questions

Search Questions:

Recent documents shows the last 10 files you've opened.  Any way to increase that to something like 20 or 50?

Dear everyone,

Ive been looking at this problem ive got, but i cant find out what would make it work. I can put numbers in a matrix, i can solve 3 unknowns with 3 equations. But now there is a variable, and the variable causes the solution to give either: 1 solution, no solution, or more than one solution. 

I can "brush the matrix" as they say it in Dutch, but im learning Maple so it can go a lot faster than manual labor. So is there a way to do it with Maple? 

Figure 4.1 gives from left to right: 1 solution, no solution, and an indefinite amount of solutions.

Questions: For what value of p, a. gives 1 solution? b. gives no solution? c. gives more than 1 solution?

Now i got Maple to tell me how much p x1, x2, and x3 was, but its not bringing me anywhere to where i can get to a solution. 

A manual method is given in this example. Manual is way way to slow for what i need to achieve, ive got some large big books with mechanics, dynamics, statics, fluid mechanics, concrete calcutions, chemistry etc. waiting for me. I can never achieve that manually in the same time as i could do it with maple. :S 

#example of how to do a solve of a matrix with 3 variables:

stelsel := {x1+x2+x3 = 4, 3*x1+x2-5*x3 = -10, 3*x1+2*x2-x3 = 1}

{x1+x2+x3 = 4, 3*x1+x2-5*x3 = -10, 3*x1+2*x2-x3 = 1}

(1)

solve(stelsel, {x1, x2, x3})

{x1 = -7+3*x3, x2 = 11-4*x3, x3 = x3}

(2)

#now the example number 15

stelsel := {(2-p)*x2+x3 = 2, -x1+2*x2+(2-p)*x3 = 0, (1-p)*x1+2*x2+2*x3 = p+3}

{(2-p)*x2+x3 = 2, -x1+2*x2+(2-p)*x3 = 0, (1-p)*x1+2*x2+2*x3 = p+3}

(3)

solve(stelsel, {x1, x2, x3})

{x1 = -(p^2-6)/(p^2-4*p+4), x2 = -(2*p-5)/(p^2-4*p+4), x3 = 1/(-2+p)}

(4)

``

Thank you!

Greetings,

The Function

Download Mapleprimes_Question_Book_2_Paragraph_4.2_Example_15.mw

Hello.

If I input 'floor(5.5)' in Maple 2018 I get the expected \lfloor 5.5 \rfloor thing.  I mean, with the special brackets typical for the floor notation.

However, if I try to place the floor notation onto a plot I get literally the string "floor(5.5)".  As in

  textplot([2, 2, 'floor(5.5)'])

The ' ' quotation marks ensure the delayed evaluation but I do not get the floor parentheses! 

Interestingly,

  textplot([2, 2, 'sqrt(5.5)'])

outputs the radical notation as expected.

How come I cannot produce a plot with the floor brackets notation in it?

Thanks

Minko

When running mint on a single .mpl file, I often get this warning

These names were used as global names but were not declared:  main_module

When having this code in the mpl file

o:=Object(main_module:-person_type);

The problem is that main_module is a module that lives in separate mpl file.

It is a differerent module in its own ,mpl file which inside it contains the person_type definition. (another module of type object)

I know I can ignore this. And have been. But my question is: how does one "declare" main_module inside this mpl file to avoid getting this warning message? Since main_module is its own module in separate file. I do not understand what does it mean to  "declare" it.

Here is the full code

A:=module()
   export boo:=proc()
      local o;
      o:=Object(main_module:-person_type);
      o:-foo();
   end proc;
end module;

And here is the command I used

"C:\Program Files\Maple 2022\bin.X86_64_WINDOWS\mint.exe" -i 2 A.mpl

Gives

Nested Procedure boo() on lines 2 to 6
  These names were used as global names but were not declared:  main_module

Later on, when I build my .mla library, all these modules are put inside the library. So it is not an issue when running the code, but I'd like to see if it is possible to add something to A.mpl to tell Maple that main_module is separate module and not to worry about it.

 

Hi,

How to keep the log function displayed in base b in my table, not its converted form in ln  ( f3 and f4 in my file) ?

Thanks.

DisplayLogQCM.mw

I have started to use CodeGeneration[C] instead of codegen[C] (depreciated). In the attached example not all variables and parameters of a C function are declared as double, although the default value is double.

If this is by design:

How can I force double without going back to codegen[C] or manual declaration of each item?

C_code_generation_with_CodeGeneration.mw

Pls see below worksheet. I want to compute the maximum number that appears in the Collatz series when starting with several starting numbers.

I calculate the maximum number appearing when starting with a specific number in procedure collatz. Then I want to calculate the maximum number appearing when calling different starting numbers in procedure collatz_max by calling collatz several times - collatz_max does not yet contain the handing of the results from collatz as the below error appears. If I only call collatz in collatz_max it return the right value, but in the loop or if I assign collatz to a variable within procedure collatz_max, calling collatz returns () as a result.

What do I do wrong?

Thanks for helping,

Oliver

Collatz.mws

Is it possible to write 

x^4 - 3*x^2 + 3

as a polynomial in 

x-1

? This is a math question first of all. If it is possible, how do we do it in Maple.

For context, I am investigating a particular line of inquiry into answering this question about computation of Taylor polynomial.

again, I wish Maple help can be better. A simple example of how to use listplot on matrix is all what is needed.

I have a matrix, 3 rows and 5 columns. I need to use listplot to generate 3 listplots on same graph. I can do this one row at a time. But when calling listplot(m) where m is the matrix, it gives error

Error, (in plots:-pointplot) incorrect number of coordinates in points data

Clicking on the above link, as usual sends me to page that says There is no help page available for this error

Help says

The listplot command also accepts a list or Matrix containing points data. great. But how?? I click on the link next to the above line, it sends me to page about pointplot. There is says that pointplot also accepts a matrix of size n by 2. Does this mean listplot matrix must also be two columns? It does not say this on the listplot help page.

For a workaround, I now generate plot of each row on its own, then use display to put them all on same graph. But the problem with this is that all have same color. I was hoping that if Maple did it all at once, it will automatically pick different color for each line as is the case with Mathematica (example below)

Is there a better way to do this? For reference, I'd like to generate similar plot like this

You see, the colors are automatically selected. With manual plotting of each row myself, I have to select the color myself and to make sure they are different and so on. It will be better if the system did all of this.

Please see attached worksheet

restart;

m:=Matrix([[1,2,5,6,9],[3,4,5,13,14],[1,3,4,10,11]])

Matrix(3, 5, {(1, 1) = 1, (1, 2) = 2, (1, 3) = 5, (1, 4) = 6, (1, 5) = 9, (2, 1) = 3, (2, 2) = 4, (2, 3) = 5, (2, 4) = 13, (2, 5) = 14, (3, 1) = 1, (3, 2) = 3, (3, 3) = 4, (3, 4) = 10, (3, 5) = 11})

plots:-listplot(m); #why this fail?

Error, (in plots:-pointplot) incorrect number of coordinates in points data

plots:-listplot(m[1,..])

#for a workaround, I can do this. But now all lines have same color which is not good.
map(n->plots:-listplot(m[n,..]),[$1..3]):
plots:-display(%)

 

Download how_to_list_plot_matrix.mw

What is the correct way to use listplot with matrix?

I have created some plots of inverse primes  like this example1.pdf .

The filled color-shape in the middle is what I want do do with all areas in this picture or in other pictures.

In other words the goal is to fill the differnt areas in the print with different colors.

So I need to find the points of the Polygons, as I have done by hand with that yellow Polygon.

A procdure that is ready will give give the crosspoints of the lines.

These are the line-coordinates (the 1st number ist the number of iterations)

2*L[1]=number of lines in L

L:=[14, [[1, 1], [1, 26]], [[1, 26], [26, 26]], [[26, 26], [26, 37]], [[26, 37], [37, 37]], [[37, 37], [37, 39]], [[37, 39], [39, 39]], [[39, 39], [39, 20]], [[39, 20], [20, 20]], [[20, 20], [20, 23]], [[20, 23], [23, 23]], [[23, 23], [23, 30]], [[23, 30], [30, 30]], [[30, 30], [30, 70]], [[30, 70], [70, 70]], [[70, 70], [70, 45]], [[70, 45], [45, 45]], [[45, 45], [45, 34]], [[45, 34], [34, 34]], [[34, 34], [34, 32]], [[34, 32], [32, 32]], [[32, 32], [32, 51]], [[32, 51], [51, 51]], [[51, 51], [51, 48]], [[51, 48], [48, 48]], [[48, 48], [48, 41]], [[48, 41], [41, 41]], [[41, 41], [41, 1]], [[41, 1], [1, 1]]]

These are the crosspoints:

cp := [[23, 26], [30, 37], [32, 37], [26, 30], [48, 45], [39, 34], [41, 34], [45, 41]]

To plot the pdf I used this code:

poly2 := [[32, 32], [34, 32], [34, 34], [39, 34], [39, 39], [37, 39], [37, 37], [32, 37]]

poly 2 is just an axample, how it looks like when its ready.

display(seq(line(op(L[i])), i = 2 .. 2*L[1] + 1), polygonplot([poly2], color = "Resene GoldenTainoi", axes = none, style = polygon), color = blue, thickness = 0.8);

So I hope, you can help me :)

This is a beautyfull way to paint a prime  by just printing the remainders of the recursive dividing of the inverse prime in lines.

Thanks a lot,

Arno

int(sin(x)/x*exp(-2*I*pi*f*x), x = -infinity .. infinity) gives the incorrect answer pi. The correct answer is pi for abs(f)

An important contribution to the integrability of non-linear differential equation, made by the E.S.Cheb.Terrab, is the Abel Inverse Riccati(AIR) equation. Recently I get an idea of make generalization of the equation using some technique introduced in Kamke's books. I write it in latex. I hope you guys can give me some advice. I have deleted some questions I asked before because I think it is more appropriate to create a post to share my ideas with you.Final_version.pdf

Dear Maple community, 

I'm trying to solve a complex expression by taking its first-order derivative and finding the optimal solution. As a result, I'm getting a rootOf expression, which is further not solved by using explicit or allvalues APIs provided in the maple. Can you guys help me where I'm going wrong?

I'm attaching my code file as well for
 

restart

U[d] := Zeta[n]*(2^(2/3)*(theta[n]*lambda*A*beta[n])^(2/3)/(4*beta[n])-delta[n]*`λA`+U[n]+alpha[n])+Zeta[g]*(U[g]-delta[g](1-lambda)*A+tau*A+R+(-2*A*theta[g]*(lambda-1)*beta[g])^(2/3)/(4*beta[g]))+tau*A-A

Zeta[n]*((1/4)*2^(2/3)*(theta[n]*lambda*A*beta[n])^(2/3)/beta[n]-delta[n]*`λA`+U[n]+alpha[n])+Zeta[g]*(U[g]-delta[g](1-lambda)*A+tau*A+R+(1/4)*(-2*A*theta[g]*(lambda-1)*beta[g])^(2/3)/beta[g])+tau*A-A

(1)

FOC := diff(U[d], A)

(1/6)*Zeta[n]*2^(2/3)*theta[n]*lambda/(theta[n]*lambda*A*beta[n])^(1/3)+Zeta[g]*(-delta[g](1-lambda)+tau-(1/3)*theta[g]*(lambda-1)/(-2*A*theta[g]*(lambda-1)*beta[g])^(1/3))+tau-1

(2)

evala(simplify(FOC))

(1/6)*2^(2/3)*(theta[n]*lambda*A*beta[n])^(2/3)*Zeta[n]/(A*beta[n])-(1/6)*Zeta[g]*(-2^(2/3)*(-A*theta[g]*(lambda-1)*beta[g])^(2/3)+6*delta[g](1-lambda)*A*beta[g]-6*tau*A*beta[g])/(A*beta[g])+tau-1

(3)

NULL

NULL

xyz := solve(FOC = 0, A, explicit)

RootOf(-Zeta[n]*theta[n]*lambda*(-_Z*lambda*beta[g]*theta[g]+_Z*beta[g]*theta[g])^(1/3)+3*delta[g](1-lambda)*(theta[n]*lambda*_Z*beta[n])^(1/3)*2^(1/3)*(-_Z*lambda*beta[g]*theta[g]+_Z*beta[g]*theta[g])^(1/3)*Zeta[g]-3*(theta[n]*lambda*_Z*beta[n])^(1/3)*2^(1/3)*(-_Z*lambda*beta[g]*theta[g]+_Z*beta[g]*theta[g])^(1/3)*tau*Zeta[g]+(theta[n]*lambda*_Z*beta[n])^(1/3)*lambda*Zeta[g]*theta[g]-3*tau*(theta[n]*lambda*_Z*beta[n])^(1/3)*2^(1/3)*(-_Z*lambda*beta[g]*theta[g]+_Z*beta[g]*theta[g])^(1/3)-(theta[n]*lambda*_Z*beta[n])^(1/3)*Zeta[g]*theta[g]+3*(theta[n]*lambda*_Z*beta[n])^(1/3)*2^(1/3)*(-_Z*lambda*beta[g]*theta[g]+_Z*beta[g]*theta[g])^(1/3))

(4)

allvalues(xyz)

RootOf(-Zeta[n]*theta[n]*lambda*(-_Z*lambda*beta[g]*theta[g]+_Z*beta[g]*theta[g])^(1/3)+3*delta[g](1-lambda)*(theta[n]*lambda*_Z*beta[n])^(1/3)*2^(1/3)*(-_Z*lambda*beta[g]*theta[g]+_Z*beta[g]*theta[g])^(1/3)*Zeta[g]-3*(theta[n]*lambda*_Z*beta[n])^(1/3)*2^(1/3)*(-_Z*lambda*beta[g]*theta[g]+_Z*beta[g]*theta[g])^(1/3)*tau*Zeta[g]+(theta[n]*lambda*_Z*beta[n])^(1/3)*lambda*Zeta[g]*theta[g]-3*tau*(theta[n]*lambda*_Z*beta[n])^(1/3)*2^(1/3)*(-_Z*lambda*beta[g]*theta[g]+_Z*beta[g]*theta[g])^(1/3)-(theta[n]*lambda*_Z*beta[n])^(1/3)*Zeta[g]*theta[g]+3*(theta[n]*lambda*_Z*beta[n])^(1/3)*2^(1/3)*(-_Z*lambda*beta[g]*theta[g]+_Z*beta[g]*theta[g])^(1/3))

(5)

NULL

Download RootOf_Maple2.mw

reference. 

I am not able to even trap this Maple exception.

Any suggestions what to do? It is not a problem if it can't solve it, but I need to at least be able to  trap the exception in order to go to the next one, else the whole program now stops when it hits this.

I used try..catch but this exception just ignores this and escapes to top level. This is not the first time I've seen Maple exception escape the try/catch. I do not understand why some do that and some not. I

Is this another bug?

interface(version);

`Standard Worksheet Interface, Maple 2022.1, Windows 10, May 26 2022 Build ID 1619613`

restart;

eq:=1 = -X*(-1/4*(-40*(X + x0)^(5/3) - 20/3*(X + x0)^(2/3)*Y - 20/3*(X + x0)^(2/3)*y0 - 32*A*(X + x0)^(1/3))/((X + x0)^(5/3)*(Y + y0)) + 5/12*(-15*(X + x0)^(8/3) - 4*(X + x0)^(5/3)*Y - 4*(X + x0)^(5/3)*y0 - 24*A*(X + x0)^(4/3) + 12*A^2)/((X + x0)^(8/3)*(Y + y0)))/(Y*(1/(Y + y0) + 1/4*(-15*(X + x0)^(8/3) - 4*(X + x0)^(5/3)*Y - 4*(X + x0)^(5/3)*y0 - 24*A*(X + x0)^(4/3) + 12*A^2)/((X + x0)^(5/3)*(Y + y0)^2)));

1 = -X*(-(1/4)*(-40*(X+x0)^(5/3)-(20/3)*(X+x0)^(2/3)*Y-(20/3)*(X+x0)^(2/3)*y0-32*A*(X+x0)^(1/3))/((X+x0)^(5/3)*(Y+y0))+(5/12)*(-15*(X+x0)^(8/3)-4*(X+x0)^(5/3)*Y-4*(X+x0)^(5/3)*y0-24*A*(X+x0)^(4/3)+12*A^2)/((X+x0)^(8/3)*(Y+y0)))/(Y*(1/(Y+y0)+(1/4)*(-15*(X+x0)^(8/3)-4*(X+x0)^(5/3)*Y-4*(X+x0)^(5/3)*y0-24*A*(X+x0)^(4/3)+12*A^2)/((X+x0)^(5/3)*(Y+y0)^2)))

try
   sol:=solve(identity(eq,X),[x0,y0]);
catch:
   print("trapped the error");
end try;

Error, (in anonymous procedure called from type/realcons) too many levels of recursion

 

Download how_to_trap.mw

Hi, I'm working on some chemistry and i'd like for the unit to me mol/l not mol/m^3. How do i fix this? I googled it and it seems like the unit is the universal standard for concentration but i wasn't taught this in school as we were told to use mol/l. 

I've tried googling and following the guides on the site but i honestly cant figure out what im doing wrong.

Picture for reference 

1 2 3 4 5 6 7 Last Page 1 of 2156