Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

I'm willing to draw phase portrait for the equation of ode fifth order

u'''''+6(u'u'''+(u'')^2)-u'''(a+2b)=0

where a and b are arbitrary constant not equal to zero and (') represents derivative w.r.t  'x'.

We have just run into a problem (loss of data), when 2 users have opened the same workbook.

Does Maple support Multiuser access to a workbook, and if not - why isn't there any lock mechanism in place to prevent that?

I've asked about this 5 years ago

 

https://www.mapleprimes.com/questions/203628-Maple-2015-No-Longer-Autoscroll-As

And according to the reply there at the time, this is supposed to have been fixed in Maple 2015.1:

It appears that the scrolling issue has been fixed in the next version of Maple

But ever since then, whenever I use Maple, the worksheet do not auto-scroll as the program runs and prints outout to the screen. I have to keep using the mouse to move the bar down to see the latest prints.

Not only that, the bar itself (on the right side) shrinks in size as more outpout is printed.

It becomes smaller and smaller with time, and becomes very hard to grab it with the mouse since it become so thin.  Once all the output is deleted, the bar become long again.

So I think there is an option to enable autoscrolling somewhere.

But I am not able to find it. I looked at options, display, general and interface and see nothing there to turn it on.

Sometimes I get the feeling that the folks at Maplesoft do not use Maple interface the way we users use it, becuase this is such a basic interface problem and any one who uses Maple interface for more than one day must see these problems all the time.

So how come these are still not fixed?  Here is a screen shot showing the bar on the right after long print out, and showing how it was before.

I use only worksheet mode, not document mode. On windows 10.

 

 

 

 

Hi, I want to use Maple to derive the eigenvectors of a 3x3 matrix A. I can obtian the 3 eigenvalues correctly by using eigenvalues(A), and the results are "u, -a+u, a+u". But after inputing eigenvectors(A), I can only get "[-a+u, 1, {r}], [a+u, 1, {r}], [u, 1, {r}]". All the 3 eigenvectors are written as {r}, I dont know why...and how can I tranform {r} into visiable results? 

I would be very appreciate if anyone could give me some help!

How to draw a circle transfomed by orthogonal affinity to oblique line and k ration ? Thank you.

 Help me please.  How can i simulate murmuration in maple.

Dear All.

I hope we are all staying safe.

Please I want to solve Sine Gordon Equation using a numerical method I constructed (non-classical), I need to compare the result of the method with the exact solution to generate the errors. However, since the exact solution has two variables, x, and t, I really don't know how to accommodate the two in my coding.

Can someone be of help in this regard?

Thank you and kind regards

 


Download Discretization_of_Sine_Gordon_Equation.mw

Download Sine_Gordon_Implementation_Trial.mw

 

Hi there,

Can someone explain to me why is it happening? If smaller equal than coloured well, if smaller then just the equal line appears. Why?

with(plots);
dm1 := .23; pm12 := .7; dm12 := .23;
inequal(pm1*(dm12+am21)-pm12*(dm1+am21) <= 0, pm1 = 0 .. 1, am21 = 0 .. 1, color = "Nautical 1");


 

inequal(pm1*(dm12+am21)-pm12*(dm1+am21) < 0, pm1 = 0 .. 1, am21 = 0 .. 1, color = "Nautical 1")

Thank you!

Hi , 

i want to find simple equation of rotate Ellipse  ( El1,EL2,El3 in the  worksheet), with calssic formula 

Ideas? Thanks

QuestionConiqueRotation.mw

Etude d'un cas particulier a := 5: b := 7: k := 9: A := [a, 0]: B := [0, b]: #A et B fixes P := [t, 0]: Q := [0, k/t]:#P et Q 2 points mobiles cir := -a*x-b*y+x^2+y^2 = 0: sol := solve(subs(y = 5, cir), x): cen := [solve(diff(cir, x)), solve(diff(cir, y))]: x0 := sol[1]: y0 := 5: M := [x0, y0]: R := sqrt(cen[1]^2+cen[2]^2): beta := arctan(diff(solve(EQ(M, cen), y), x)): Recherche des valeurs de t pour que les 2 droites soient perpendiculaires eq := t^2*(y0-b)+t*(a*b-a*y0+b*x0-k)-x0*(a*b-k) = 0; sol := solve(eq, t); t := sol[1]; tp := sol[2]; P1 := [t, 0]; Q1 := [0, k/t]; PQ1 := simplify(x*(-a*b+b*t+k)+y*t*(t-a)-t*(-a*b+b*t+k)) = 0:#1ere tangente PQ2 := simplify(x*(-a*b+b*tp+k)+y*tp*(tp-a)-tp*(-a*b+b*tp+k)) = 0:#2ième tangente P2 := [tp, 0]; Q2 := [0, k/tp]; CIR := implicitplot(cir, x = -4 .. 8, y = -4 .. 12, color = red); Fig := proc (alpha) local Dr1, DR1, Dr2, DR2, N, u0, v0, Po, t, tp, sol; global a, b, k, cen, R; u0 := cen[1]+R*cos(alpha); v0 := cen[2]+R*sin(alpha); N := [u0, v0]; sol := solve(t^2*(v0-b)+t*(b*u0-a*v0+a*b-k)-u0*(a*b-k) = 0, t); t := sol[1]; tp := sol[2]; Dr1 := simplify(x*(-a*b+b*t+k)+y*t*(t-a)-t*(-a*b+b*t+k)) = 0; DR1 := implicitplot(Dr1, x = -4 .. 8, y = -4 .. 12, color = brown); Dr2 := simplify(x*(-a*b+b*tp+k)+y*tp*(tp-a)-tp*(-a*b+b*tp+k)) = 0; DR2 := implicitplot(Dr2, x = -4 .. 8, y = -4 .. 12, color = pink); Po := pointplot([N[]], symbol = solidcircle, color = [black], symbolsize = 8); display([Po, DR1, DR2]) end proc; DrPQ1 := implicitplot(PQ1, x = -4 .. 22, y = -4 .. 12, color = blue); DrPQ2 := implicitplot(PQ2, x = -4 .. 22, y = -4 .. 12, color = blue); Points := pointplot([A[], B[], M[], P1[], P2[], Q1[], Q2[], cen[]], symbol = solidcircle, color = [green], symbolsize = 10); T := plots:-textplot([[A[], "A"], [B[], "B"], [M[], "M"], [P1[], "P1"], [P2[], "P2"], [Q1[], "Q1"], [Q2[], "Q2"], [cen[], "cen"]], font = [times, 10], align = {below, left}); n := 19; display([seq(Fig(2*i*Pi/n), i = 0 .. n), Fig(beta), CIR, DrPQ1, DrPQ2, Points, T], scaling = constrained, size = [500, 500]); I would find out the focus of the ellipse. Thank you.
Hello, dear forum users!

Does anyone use the method of homotopy analysis (HAM) and the NOPH package in their work. (moderator: link)
 


It seems to me that only HAM can help.

I ask for help if someone has already mastered.

There are no developments, as I do not own the NOPH package.

When I type "with(combinat)" in Maple Linux, I got the following erros message:

"Error, invalid input: with expects its 1st argument, pname, to be of type {`module`, package}, but received combinat"

So how to make it work? 

 

Many thanks.

code_(1).mw
 

 

 

Hello! hope you all be fine. i m running a code for `Anew Technique of intial Boundary Value Problem Using ADM `(Elaf*Jaafar*Ali) but my code is for positive integer and i want to run it for fractional can anyone help me please code file is attached 

The pdf shows the derivation procedure of Image from 1.2.1 to 1.2.5. However, when we obtain the (7) in pdf, which is the 

 

expansion of 1.2.4, we hope to see the hamiltonian operator in the output.  We obtain the (10) in pdf (Maple). Actually, it is the expansion of 1.2.5. However, it is too cockamamie compared with that (1.2.5) in the textbook.

 

Hence, I hope to show the hamiltonian operator in the (10) in PDF. I hope to see the simplest type of equation.

How could I do? Appreciate the kind help from anyone!

The Maple code is also uploaded.
 

 

 

 

restart*with(Physics[Vectors])

Step(mathematicalnotation = true)*SetCoordinates('cartesian'[i, j, k])

Step(mathematicalnotation = true)*SetCoordinates(cartesian[i, j, k])

(1)

``

NULL``

"  UU(x,y,z):=ux(x,y,z)*(_i)+uy(x,y,z)*_j+uz(x,y,z)*_k; "

proc (x, y, z) options operator, arrow, function_assign; Physics:-Vectors:-`+`(Physics:-Vectors:-`+`(ux(x, y, z)*_i, uy(x, y, z)*_j), uz(x, y, z)*_k) end proc

(2)

"&nabla;(pressure(x,y,z));  Laplacian(UU(x,y,z));  BB(x,y,z):=0*_i+0*_j+10*_k"

proc (x, y, z) options operator, arrow, function_assign; Physics:-Vectors:-`+`(Physics:-Vectors:-`+`(0*_i, 0*_j), 10*_k) end proc

(3)

-VectorCalculus[Nabla](pressure(x, y, z))+Laplacian(mu*UU(x, y, z))+BB(x, y, z) = 0*_i+0*_j+0*_k

_i*(-(diff(pressure(x, y, z), x))+mu*(diff(diff(ux(x, y, z), x), x))+mu*(diff(diff(ux(x, y, z), y), y))+mu*(diff(diff(ux(x, y, z), z), z)))+_j*(-(diff(pressure(x, y, z), y))+mu*(diff(diff(uy(x, y, z), x), x))+mu*(diff(diff(uy(x, y, z), y), y))+mu*(diff(diff(uy(x, y, z), z), z)))+_k*(-(diff(pressure(x, y, z), z))+mu*(diff(diff(uz(x, y, z), x), x))+mu*(diff(diff(uz(x, y, z), y), y))+mu*(diff(diff(uz(x, y, z), z), z))+10) = 0

(4)

``

Divergence(-VectorCalculus[Nabla](pressure(x, y, z))+Laplacian(mu*UU(x, y, z))+BB(x, y, z) = 0*_i+0*_j+0*_k)

-(diff(diff(pressure(x, y, z), x), x))+mu*(diff(diff(diff(ux(x, y, z), x), x), x))+mu*(diff(diff(diff(ux(x, y, z), x), y), y))+mu*(diff(diff(diff(ux(x, y, z), x), z), z))-(diff(diff(pressure(x, y, z), y), y))+mu*(diff(diff(diff(uy(x, y, z), x), x), y))+mu*(diff(diff(diff(uy(x, y, z), y), y), y))+mu*(diff(diff(diff(uy(x, y, z), y), z), z))-(diff(diff(pressure(x, y, z), z), z))+mu*(diff(diff(diff(uz(x, y, z), x), x), z))+mu*(diff(diff(diff(uz(x, y, z), y), y), z))+mu*(diff(diff(diff(uz(x, y, z), z), z), z)) = 0

(5)

 

algsubs(Divergence(UU(x, y, z)) = 0, %)

-(diff(diff(pressure(x, y, z), x), x))+mu*(-(diff(diff(diff(uy(x, y, z), x), x), y))-(diff(diff(diff(uz(x, y, z), x), x), z)))+mu*(-(diff(diff(diff(uy(x, y, z), y), y), y))-(diff(diff(diff(uz(x, y, z), y), y), z)))+mu*(-(diff(diff(diff(uy(x, y, z), y), z), z))-(diff(diff(diff(uz(x, y, z), z), z), z)))-(diff(diff(pressure(x, y, z), y), y))+mu*(diff(diff(diff(uy(x, y, z), x), x), y))+mu*(diff(diff(diff(uy(x, y, z), y), y), y))+mu*(diff(diff(diff(uy(x, y, z), y), z), z))-(diff(diff(pressure(x, y, z), z), z))+mu*(diff(diff(diff(uz(x, y, z), x), x), z))+mu*(diff(diff(diff(uz(x, y, z), y), y), z))+mu*(diff(diff(diff(uz(x, y, z), z), z), z)) = 0

(6)

simplify(%)

-(diff(diff(pressure(x, y, z), x), x))-(diff(diff(pressure(x, y, z), y), y))-(diff(diff(pressure(x, y, z), z), z)) = 0

(7)

Laplacian(-VectorCalculus[Nabla](pressure(x, y, z))+Laplacian(mu*UU(x, y, z))+BB(x, y, z) = 0*_i+0*_j+0*_k)

(mu*(diff(diff(diff(diff(ux(x, y, z), x), x), x), x))+2*mu*(diff(diff(diff(diff(ux(x, y, z), x), x), y), y))+2*mu*(diff(diff(diff(diff(ux(x, y, z), x), x), z), z))+mu*(diff(diff(diff(diff(ux(x, y, z), y), y), y), y))+2*mu*(diff(diff(diff(diff(ux(x, y, z), y), y), z), z))+mu*(diff(diff(diff(diff(ux(x, y, z), z), z), z), z))-(diff(diff(diff(pressure(x, y, z), x), x), x))-(diff(diff(diff(pressure(x, y, z), x), y), y))-(diff(diff(diff(pressure(x, y, z), x), z), z)))*_i+(mu*(diff(diff(diff(diff(uy(x, y, z), x), x), x), x))+2*mu*(diff(diff(diff(diff(uy(x, y, z), x), x), y), y))+2*mu*(diff(diff(diff(diff(uy(x, y, z), x), x), z), z))+mu*(diff(diff(diff(diff(uy(x, y, z), y), y), y), y))+2*mu*(diff(diff(diff(diff(uy(x, y, z), y), y), z), z))+mu*(diff(diff(diff(diff(uy(x, y, z), z), z), z), z))-(diff(diff(diff(pressure(x, y, z), x), x), y))-(diff(diff(diff(pressure(x, y, z), y), y), y))-(diff(diff(diff(pressure(x, y, z), y), z), z)))*_j+(mu*(diff(diff(diff(diff(uz(x, y, z), x), x), x), x))+2*mu*(diff(diff(diff(diff(uz(x, y, z), x), x), y), y))+2*mu*(diff(diff(diff(diff(uz(x, y, z), x), x), z), z))+mu*(diff(diff(diff(diff(uz(x, y, z), y), y), y), y))+2*mu*(diff(diff(diff(diff(uz(x, y, z), y), y), z), z))+mu*(diff(diff(diff(diff(uz(x, y, z), z), z), z), z))-(diff(diff(diff(pressure(x, y, z), x), x), z))-(diff(diff(diff(pressure(x, y, z), y), y), z))-(diff(diff(diff(pressure(x, y, z), z), z), z)))*_k = 0

(8)

algsubs(-(diff(pressure(x, y, z), x, x))-(diff(pressure(x, y, z), y, y))-(diff(pressure(x, y, z), z, z)) = 0, %)

(mu*(diff(diff(diff(diff(ux(x, y, z), x), x), x), x))+2*mu*(diff(diff(diff(diff(ux(x, y, z), x), x), y), y))+2*mu*(diff(diff(diff(diff(ux(x, y, z), x), x), z), z))+mu*(diff(diff(diff(diff(ux(x, y, z), y), y), y), y))+2*mu*(diff(diff(diff(diff(ux(x, y, z), y), y), z), z))+mu*(diff(diff(diff(diff(ux(x, y, z), z), z), z), z)))*_i+(mu*(diff(diff(diff(diff(uy(x, y, z), x), x), x), x))+2*mu*(diff(diff(diff(diff(uy(x, y, z), x), x), y), y))+2*mu*(diff(diff(diff(diff(uy(x, y, z), x), x), z), z))+mu*(diff(diff(diff(diff(uy(x, y, z), y), y), y), y))+2*mu*(diff(diff(diff(diff(uy(x, y, z), y), y), z), z))+mu*(diff(diff(diff(diff(uy(x, y, z), z), z), z), z)))*_j+(mu*(diff(diff(diff(diff(uz(x, y, z), x), x), x), x))+2*mu*(diff(diff(diff(diff(uz(x, y, z), x), x), y), y))+2*mu*(diff(diff(diff(diff(uz(x, y, z), x), x), z), z))+mu*(diff(diff(diff(diff(uz(x, y, z), y), y), y), y))+2*mu*(diff(diff(diff(diff(uz(x, y, z), y), y), z), z))+mu*(diff(diff(diff(diff(uz(x, y, z), z), z), z), z)))*_k = 0

(9)

collect(%, mu)

((diff(diff(diff(diff(ux(x, y, z), x), x), x), x)+2*(diff(diff(diff(diff(ux(x, y, z), x), x), y), y))+2*(diff(diff(diff(diff(ux(x, y, z), x), x), z), z))+diff(diff(diff(diff(ux(x, y, z), y), y), y), y)+2*(diff(diff(diff(diff(ux(x, y, z), y), y), z), z))+diff(diff(diff(diff(ux(x, y, z), z), z), z), z))*_i+(diff(diff(diff(diff(uy(x, y, z), x), x), x), x)+2*(diff(diff(diff(diff(uy(x, y, z), x), x), y), y))+2*(diff(diff(diff(diff(uy(x, y, z), x), x), z), z))+diff(diff(diff(diff(uy(x, y, z), y), y), y), y)+2*(diff(diff(diff(diff(uy(x, y, z), y), y), z), z))+diff(diff(diff(diff(uy(x, y, z), z), z), z), z))*_j+(diff(diff(diff(diff(uz(x, y, z), x), x), x), x)+2*(diff(diff(diff(diff(uz(x, y, z), x), x), y), y))+2*(diff(diff(diff(diff(uz(x, y, z), x), x), z), z))+diff(diff(diff(diff(uz(x, y, z), y), y), y), y)+2*(diff(diff(diff(diff(uz(x, y, z), y), y), z), z))+diff(diff(diff(diff(uz(x, y, z), z), z), z), z))*_k)*mu = 0

(10)

``

``

 


``

 

Download 20200619.mw

 

 

Download 20200619-1.mw

20200619.pdf

 

Error, (in Compiler:-Setup) unable to write to initialization file, C:\Program Files\Maple 2020\bin.X86_64_WINDOWS/launch.ini, possibly due to insufficient permissions
NULL;

when choose bat file

but i see sufficient permission

any w should i s

choose a particular compiler or even visual studio i have in my system

 

where can i find if anything else

First 85 86 87 88 89 90 91 Last Page 87 of 2218