Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

it took a while to track this. When I run these 3 ode's in this sequence, Maple gives 

Error, (in Engine:-Dispatch) too many levels of recursion

and it has to be this sequence.  I also found I might have to run the code 2 or 3 times to get this error if it does not show first time. 

Will show a movie below also.

Here is the code

restart;
ode1:=cos(y(x)) - x*sin(y(x))*diff(y(x), x) = sec(x)^2:
ode2:=y(x)*sin(x/y(x)) + x*cos(x/y(x)) - 1 + (x*sin(x/y(x)) - x^2*cos(x/y(x))/y(x))*diff(y(x), x) = 0:
ode3:=x/(x^2 + y(x)^2) + y(x)/x^2 + (y(x)/(x^2 + y(x)^2) - 1/x)*diff(y(x), x) = 0:

func:=y(x);
the_ic:={};
the_time_limit:=60;

advisor:=timelimit(the_time_limit,DEtools:-odeadvisor(ode1)):
solution_by_CAS := timelimit(the_time_limit,:-dsolve([ode1,y(0)=0],func,singsol=all)):
result_of_ode_test := timelimit(the_time_limit,:-odetest~([solution_by_CAS],ode1)):

advisor:=timelimit(the_time_limit,DEtools:-odeadvisor(ode2)):
solution_by_CAS := timelimit(the_time_limit,:-dsolve([ode2,y(0)=0],func,singsol=all)):
result_of_ode_test := timelimit(the_time_limit,:-odetest~([solution_by_CAS],ode2)):

advisor:=timelimit(the_time_limit,DEtools:-odeadvisor(ode3)):
solution_by_CAS := timelimit(the_time_limit,:-dsolve([ode3,y(1)=0],func,singsol=all)):
result_of_ode_test := timelimit(the_time_limit,:-odetest~([solution_by_CAS],ode3));


And here is a movie (this time it happend on second run) sometimes it happens on first run. so it is random.

 

Maple 2020.1

Why does it happen and is there a work around?


 

Noticed few errors generated on some ode's. I do not know if these are expected or known, I thought may be someone might want to check them.

restart;
ode1:=(a^2*x+y(x)*(x^2-y(x)^2))*diff(y(x),x)+x*(x^2-y(x)^2) = a^2*y(x);
ode2:=x^7*y(x)^2*diff(y(x),x)^3+(1-3*x^6*y(x)^3)*diff(y(x),x)^2+3*x^5*y(x)^4*diff(y(x),x)-x^4*y(x)^5 = 0;
ode3:=diff(y(x),x)^4+x*diff(y(x),x)-3*y(x) = 0;

sol:=timelimit(60,:-dsolve(ode1,y(x))):

Error, (in dsolve) numeric exception: division by zero
 

sol:=timelimit(60,:-dsolve(ode2,y(x))):
timelimit(60,:-odetest([sol][5],ode2))

Error, (in IntegrationTools:-Indefinite:-AlgebraicFunction) argument is not an algebraic
 

sol:=timelimit(60,:-dsolve(ode3,y(x))):
timelimit(60,:-odetest(sol[1],ode3));
timelimit(60,:-odetest(sol[2],ode3))

Error, (in odetest) received not an ODE w.r.t x(_T)
Error, (in odetest) received not an ODE w.r.t y(_T)
 

This last one is only due to how the solution is given. So I do not think it is a real problem.

Maple 2020.1, Physics 708

windows 10



Maple Worksheet - Error
Failed to load the worksheet /maplenet/convert/check.mw .
 

Download check.mw

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

Hello,

I try import a big Excel file with the first row as column labels for performing statistical etc.

I use this command :

data := Import(FileTools:-JoinPath(["d:", "directory", "subset2.xlsx"], platform = "windows"))

i need a dataframe for performing statistical with column label as first row, but i have this error :

Error, (in Import) invalid input: MatrixToDataFrame uses a 3rd argument, columns, which is missing

I can not understand what is this error...so, what i can do ?

Thanks, regards


This is my notebook

 

restart

with(DifferentialGeometry)

DGsetup([t, rho, theta, phi, `#mover(mi("&theta;",fontstyle = "normal"),mo("&uminus0;"))`, `#mover(mi("&phi;",fontstyle = "normal"),mo("&uminus0;"))`, `#mover(mi("&theta;",fontstyle = "normal"),mo("&tilde;"))`, `#mover(mi("&phi;",fontstyle = "normal"),mo("&tilde;"))`], M)

M > 

``

M > 

evalDG(`&w`(dt, `d&rho;`))

_DG([["form", M, 1], [[[1], `d&rho;`]]])

(1)
M > 

evalDG(`&s`(dt, `d&rho;`))

_DG([["tensor", M, [["cov_bas"], []]], [[[1], `d&rho;`]]])

(2)
M > 

 

 

 

``

I don't understand why computing the wedge product of two differential forms gives me an output which prints the product of both. Why  don't it just return d rho ^ dt?

Dear Maple Users @acer @Carl Love @Kitonum @Preben Alsholm ,

Greetings.

How to plot a function "Am" for various values of "kt"(eg: kt=-1..1) at a point x=0.

Am vs kt.

The code has provided below.

Waiting for your response.

ktplot.mw

 

Hi experts

with(plots):
#epsilon:=5*Pi/2:
gammaa:=0.4:
gammao:=0.1:
gammab:=0.002:
omegab:=100:
X:=60:
Omega:=6*Pi:
g:=10:
sigmag:=5*gammaa:
a:=1/(sigmag*sqrt(2*Pi)):

#deltao:=1:
D7:=evalf(a*int(g*Omega*(gammao/2-I*Delta)*exp(-0.5*((Delta-deltao)/sigmag)^2)/(gammao^2/4+Delta^2),Delta=-infinity..infinity)):
D8:=evalf(0.5*a*int(g^2*gammao*exp(-0.5*((Delta-deltao)/sigmag)^2)/(gammao^2/4+Delta^2), Delta=-infinity..infinity)):
D9:=evalf(a*int(g^2*Delta*exp(-0.5*((Delta-deltao)/sigmag)^2)/(gammao^2/4+Delta^2), Delta=-infinity..infinity)):
b1 := (deltao-2*X^2*omegab*Y/(gammab^2+omegab^2)-D9)^2:

 f:=Y*((gammaa+D8)^2+(b1))=Im(D7)^2+(epsilon-Re(D7))^2:
#P1:=implicitplot(f,deltao=-70..50,Y=0..1,numpoints=100,axes=boxed,thickness=2,color=black,font=[1,1,20],tickmarks=[4, 3],linestyle=1);

implicitplot3d(f,Y = 0 .. 1,deltao=-40..40,epsilon=0..0.4,numpoints=100,labels=[Y,deltao,E],tickmarks=[3,3,3], style = surface, color = "Niagara Azure");

The above code to plot the implicitplot3d ,

first the figures are not smooth

second thing is how to make contourplot at epsilon= 5*Pi/2.

I appreciated any comments


 

First 11 12 13 14 15 16 17 Last Page 13 of 1777