MaplePrimes Questions

On my journey of discovery through the world of Maple, I would like to ask for help again:
How are the properties of variables and the indexing of sequences handled/determined? For this purpose, I chose an old problem from a challenging MO as an exercise and tried it in the attached file, which, of course, failed.

restart

a[i] := 2*cos(t/2^i)-1

2*cos(t/2^i)-1

(1)

b[n] := product(a[i], i = 1 .. n)

product(2*cos(t/2^i)-1, i = 1 .. n)

(2)

limit(b[n], n = infinity)

limit(product(2*cos(t/2^i)-1, i = 1 .. n), n = infinity)

(3)

simplify(limit(b[n], n = infinity))

limit(product(2*cos(t*2^(-i))-1, i = 1 .. n), n = infinity)

(4)

NULL

Download test.mw

Regarding to this equation i want to do the same as paper did, but i am not sure how i can determine the both  point in plot and how thus point show effect of system of equation and  how find jacobian in the two different point  i think he find it in general but i think we have to find it for each one  of equalibriom point  or we just find in general .
thus plot also are emazing did he use special code in matlab or mathematica? there is some app for plot thus kind of phase portrait but they not like this this is must have a special code for plotting in matlab or mathematica

restart

with(PDEtools)

undeclare(prime, quiet)

with(LinearAlgebra)

declare(u(x, t), quiet); declare(U(xi), quiet); declare(V(xi), quiet)

NULL

DU := V

V

(1)

DV := U^2*beta/(delta*k^2)+(alpha*w/(delta*k^3)+1/(delta*k^2)-w^2/(delta*k^4))*U

U^2*beta/(delta*k^2)+(alpha*w/(delta*k^3)+1/(delta*k^2)-w^2/(delta*k^4))*U

(2)

var := U, V

U, V

(3)

J := Student:-MultivariateCalculus:-Jacobian([DU, DV], [var], 'output' = 'matrix')

Matrix(%id = 36893490902008029780)

(4)

Eigenvalues(J, implicit, output = 'list')

[RootOf(_Z^2*delta*k^4-2*U*beta*k^2-alpha*k*w-k^2+w^2, index = 1), RootOf(_Z^2*delta*k^4-2*U*beta*k^2-alpha*k*w-k^2+w^2, index = 2)]

(5)

Eigenvalues(J)

Vector[column](%id = 36893490902007993756)

(6)

NULL

NULL

DU1 := V

V

(7)

V = solve(DU1, V)

V = 0

(8)

DV1 := U^2*F[2]+U*F[1]

U^2*F[2]+U*F[1]

(9)

U = solve(DV1, U)

U = (0, -F[1]/F[2])

(10)

so*we*have*three*equalibriom*point*(0, 0)

so*we*have*three*equalibriom*point*(0, 0)

(11)

J1 := Student:-MultivariateCalculus:-Jacobian([DU1, DV1], [var], 'output' = 'matrix')

Matrix(%id = 36893490902007256836)

(12)

Eigenvalues(J1, implicit, output = 'list')

[RootOf(-2*U*F[2]+_Z^2-F[1], index = 1), RootOf(-2*U*F[2]+_Z^2-F[1], index = 2)]

(13)

Eigenvalues(J1)

Vector[column](%id = 36893490901938203092)

(14)
 

NULL

Download bi-1.mw

The attached problem is from a 1988 MO. It can be solved using complete induction, paper, and pencil, and with some effort, yields a simple answer. It's quite challenging to do by hand, but with "derive", it only takes three lines and a fraction of a second. Mow test.mw

restart

"puzzle(n):=(∑)(cos(k*Pi/(2*n+1))^())^(4);"

Download test.mw

I can't do it with Maple because I'm doing something wrong again. Therefore, I'm asking for help.

Recently, I was surprised by the fact that it is not possible to create a table, where the values of entries are tables.

Toy example:

T1:=table(["1"=table(["first"=12, "second"=15]), "2"= table(["first"=9, "second"=7])]);

The code

T1["1"];

doesn't return table(["first"=12, "second"=15])

If table has lists in right-hand side

T1:=table(["1"=[12, 15], "2"= [9, 7]]);

all works properly.

 Based on the type of coefficients in the linear equations, SolveTools[Linear] provide several method including method = Rational, Polynomial etc. 
The Polynomial method of SolveTools[Linear] however, cannot be directly called by SolveTools[Linear]. The interpreter complains that no such methd called "Polynomial".

This could only be solved if I import the SolveTools package in advance and call Linear.

I hope Maple could solve this issue. 

I'm trying to reproduce a manual asymptotic analysis (see the attached pdf file) in Maple for a two-soliton solution. Specifically, I want to evaluate the limit of a function (e.g., r[2]r[2]r[2] or ∂xq[2]\partial_x q[2]∂x​q[2]). How can I properly perform the limit a2→+−∞ as t​→+−∞ in Maple, either by substitution or by reparametrization, in order to study the asymptotic behavior of a multi-variable expression symbolically? 

restart

with(Student[Calculus1])

lambda1 := I*mu1; lambda2 := I*mu2; a1 := -2*x/mu1+mu1*t; a2 := -2*x/mu2+mu2*t

numer_r := lambda2*cosh(a1)-lambda1*cosh(a2)

denom_r := (lambda1^2+lambda2^2)*cosh(a1)*cosh(a2)-2*lambda1*lambda2*(1+sinh(a1)*sinh(a2))

r2 := I*(-lambda1^2+lambda2^2)*numer_r/denom_r

numer_dq := (sinh(a1)-sinh(a2))^2; denom_dq := ((lambda1^2+lambda2^2)*cosh(a1)*cosh(a2)-2*lambda1*lambda2*(1+sinh(a1)*sinh(a2)))^2

dq2 := 1-2*(lambda1^2-lambda2^2)^2*numer_dq/denom_dq

limit_r2 := limit(r2, a2 = -infinity); simplify(limit_r2)

Error, invalid input: limit expects its 2nd argument, p, to be of type Or(name = algebraic,set(name = algebraic),list(name = algebraic)), but received -2*x/mu2+mu2*t = -infinity

 

limit_r2

(1)

limit_r2_pos := limit(r2, a2 = infinity); simplify(limit_r2_pos)

Error, invalid input: limit expects its 2nd argument, p, to be of type Or(name = algebraic,set(name = algebraic),list(name = algebraic)), but received -2*x/mu2+mu2*t = infinity

 

limit_r2_pos

(2)

NULL

Download asymptotic.mw cooocp_(2).pdf

Any suggestions for reformulating the limit or change of variables would be appreciated. 

I have solution of this problem but how I can plot simulation on Maple?Thanks, In attach my problem

Good day, everyone

I am trying to code HPM but it's giving me the error code "Error, invalid subscript selector
" once I increase the number of iterations above 2. 

Attached below is the code. 

HPM.mw

Thanks 

In a recent MaplePrimes question: https://mapleprimes.com/questions/240543-Maple-20251-On-Windows-11-How-To-Resize 

a discussion of interrupting a computation also came up.

I gave some comments. In one of these I wanted to include a link to a screenshot. 
I tried several times. It seemed that I only had the choice between the full image or a thumb nail.
So I chose the latter. That appeared to me as worthless. I tried a couple of times more, but decided to delete all 3.

Then I wrote a comment without a screenshot. 

It looks silly to have these empty (and anonymous) "comments" still around.
Is there any way to remove these?

Hello, when I have a complex conjugate in the denominator of a fraction, the conjugate bar symbol and the division line of the fraction overlap. When this happens, it is not apparent that the variable is a conjugate, so I have had to just type the command conjugate() in order to make this clear.

Reference the linked sample document that illustrates this: Conjugate.flow

I would appreciate any suggestion on how to use the conjugate symbol in this case. I would prefer the symbol wherever possible rather than the text command for clarity in the flow document.

Thanks

Hello, I have a procedure "timing" in the code editor that I am calling from the main canvas of a flow document. I would like to print the procedure description text string from the Description option in the procedure using the Describe command. However, Describe(timing) does not return anything.

The flow document is at this link: Procedure-Describe.flow

I would appreciate any suggestions.

Thanks

I am new to Windows 11 (edit: I was forced to migrate from 10).

I am looking for an option to resize the Window that the interupt button becomes accessible (see red box)

Also: there are no window controls (maximize minimize) visble . (edit: the title bar is clipped off).

Window key + arrow key does not work.

I have also tried the Snap Layout option to resize but cannot find Maple 2025.1 under the selectable tasks.

What else can I try?

I checked the ConsistencyTest of the system of equations but no output with 'true' or 'False'. Is it not work in 'DEtools'? Download consistency.mw

I have a piecewise density function (f).  I am trying to find the median value.

I have tried the Median function and the Percentile function, but neither work for me.  I am not sure why.

I have also tried to integrate the density function on the (0,x) interval such that the area under the curve is 1/2 and then solve for x.  This works for simple problems, but not the one attached.

There is something very simple that I'm not doing, but I am not sure what.

restart

kernelopts(version)

`Maple 2023.2, X86 64 WINDOWS, Nov 24 2023, Build ID 1762575`

(1)

with(Statistics)

 

(2)

f := 2*t*piecewise(t <= 0, 0, t <= 1, t^2+Pi-4*t, t <= sqrt(2), -(sqrt(t^2-1)*t^2+2*sqrt(t^2-1)*arcsin((t^2-2)/t^2)-4*t^2+2*sqrt(t^2-1)+4)/sqrt(t^2-1), sqrt(2) < t, 0)

2*t*piecewise(t <= 0, 0, t <= 1, t^2+Pi-4*t, t <= 2^(1/2), -((t^2-1)^(1/2)*t^2+2*(t^2-1)^(1/2)*arcsin((t^2-2)/t^2)-4*t^2+2*(t^2-1)^(1/2)+4)/(t^2-1)^(1/2), 2^(1/2) < t, 0)

(3)

plot(f, t = 0 .. sqrt(2))

 

Median(f)

2*t*piecewise(t <= 0, 0, t <= 1, t^2+Pi-4*t, t <= 2^(1/2), -((t^2-1)^(1/2)*t^2+2*(t^2-1)^(1/2)*arcsin((t^2-2)/t^2)-4*t^2+2*(t^2-1)^(1/2)+4)/(t^2-1)^(1/2), 2^(1/2) < t, 0)

(4)

Percentile(f, 50)

2*t*piecewise(t <= 0, 0, t <= 1, t^2+Pi-4*t, t <= 2^(1/2), -((t^2-1)^(1/2)*t^2+2*(t^2-1)^(1/2)*arcsin((t^2-2)/t^2)-4*t^2+2*(t^2-1)^(1/2)+4)/(t^2-1)^(1/2), 2^(1/2) < t, 0)

(5)

myMedian := solve(int(f, t = 0 .. x) = 1/2, x)

Warning, solutions may have been lost

 

(6)

NULL

Download Median.mw

When calculating limits of real-valued functions, sometimes (especially in competitions) tricky approaches are taken using pen and paper. I repeatedly encountered the simple conclusion that, for example, for natural k, the value sin(k*pi) = 0. Thus, the function value is determined logically without specifying a specific number. There are numerous other examples of this that can easily be constructed.
My question after unsuccessful attempts using "assume" is:
How, for example, does Maple determine the value of sin(k*pi) from the assumption "k is natural" alone? Are such prominent values ​​implemented in tables?

4 5 6 7 8 9 10 Last Page 6 of 2427