Maple Questions and Posts Maple Questions and Posts Feed

These are Posts and Questions associated with the product, Maple


Im trying to study some questions and I'm using maple to verify my answers.

Theres a few polynomial factoring questions and linear equation questions Im trying to get

maple to show its solutions steps using showsolution() no matter where I put it  the function wont work.

Ive switched between math/text functions. Im still pretty new to maple but I can't find any information on how to do it

on the web/youtube.


Thanks in advance!

For problem simplify(abs(1-b)+abs(1+b)), I want maple to take out the abs and get results for different ranges of b.

I have double indexed functiions f[j,k] of one variable and double indexed coefficients a[j,k].

I want my print do look like a[1,1]f11+a[1,2]f12 that is, the values of a[j,k] should appear beside the functions' names, like

7f11+2f12-3f21 etc.

Thank yopu for any help


     I'm trying to solve this PDE, and Maple 2015 gives me a solution quickly. I can test the solution with pdetest() and this verifies that it works. However, when I try to verify this myself I don't get zero. Is there some trick pdetest() is using to that I am missing? Or is pdetest() wrong in this case?



eq := I*exp(-(2*I)*k*t)*k*sin(theta)*r^2*cos(theta)^3+4*exp(-(2*I)*k*t)*r*cos(theta)^3+2*(diff(Vr(t, r, theta), theta, theta))*cos(theta)*exp(-I*k*(sin(theta)*r+t))-6*(diff(Vr(t, r, theta), theta))*sin(theta)*exp(-I*k*(sin(theta)*r+t))-4*Vr(t, r, theta)*cos(theta)*exp(-I*k*(sin(theta)*r+t))-4*exp(-(2*I)*k*t)*r*cos(theta);

I*exp(-(2*I)*k*t)*k*sin(theta)*r^2*cos(theta)^3+4*exp(-(2*I)*k*t)*r*cos(theta)^3+2*(diff(diff(Vr(t, r, theta), theta), theta))*cos(theta)*exp(-I*k*(sin(theta)*r+t))-6*(diff(Vr(t, r, theta), theta))*sin(theta)*exp(-I*k*(sin(theta)*r+t))-4*Vr(t, r, theta)*cos(theta)*exp(-I*k*(sin(theta)*r+t))-4*exp(-(2*I)*k*t)*r*cos(theta)


sol := pdsolve(eq);

Vt(t, r, theta) = _F2(t, r)/cos(theta)^2+sin(theta)*_F1(t, r)/cos(theta)^2-((1/2)*I)*(cos(theta)^2*k^2*r^2-2)*exp(I*(sin(theta)*r-t)*k)/(k^3*r^2*cos(theta)^2)


pdetest(sol, eq);



eq2 := eval(eq, Vr(t,r,theta) = rhs(sol)):
eq2 := simplify(%);

-((1/2)*I)*exp(-(2*I)*k*t)*k*r^2*cos(theta)^3+2*exp(-(2*I)*k*t)*r*sin(theta)*cos(theta)-3*(diff(Vt(t, r, theta), theta))*sin(theta)*exp(-I*k*(sin(theta)*r+t))-2*Vt(t, r, theta)*cos(theta)*exp(-I*k*(sin(theta)*r+t))+(diff(diff(Vt(t, r, theta), theta), theta))*cos(theta)*exp(-I*k*(sin(theta)*r+t))


evalb(eq2 = 0);









I'm writing to ask how to equalize the coefficients of two multivariate polynomials. In particluar, I have two polynomials whose arguments are ln(E),ln(K),ln(L) (their levels, squared levels and interaction terms). The first one is:


the second one is:


I would like to know if it is possible to equalize the coefficients of the two polynomials and find the following system:

v*a*b = x_1, -v*(a-1) x_3, -v*a*(-1+b) = x_2, a*b*v*(b*rho*a-b*rho+g*(-1+b)) = x_11, v*rho*a*(a-1) = x_33, v*a*(rho*(-1+b)*a-rho*(-1+b)+b*g)*(-1+b) = x_22, -a*v*rho*(a-1)*b = x_13, -a*v*(a*rho-rho*u+g)*b*(-1+b) = x_12, a*v*u*rho*(a-1)*(-1+b) = x_23

I tried using "coeffs" and creating a sequence of values for x but then I don't know how to equalize them.

Thank you very much in advance for your time,


Hi everyone.

I'm trying to solve the following PDE


but I'm getting this error:

Error, (in pdsolve/numeric/process_IBCs) initial/boundary conditions can only contain derivatives which are normal to the boundary, got (D[1, 1](w))(x, 0)

The PDE represents the bending of a thin plate.

See File:

into the "Ask a Question" window?

Nothing to add

Actually I want to ask something else.



(Maple 2015)

For the simple ODE with initial condition
dsolve({ diff(y(t),t) = y(t)^2 - y(t)^3, y(0) = 1/10 }, y(t));

dsolve produces two different answers, almost randomly (even after restart or after closing Maple and reloading the worksheet). Namely:





but this simplifies to (2), so it's not a "true" bug.

Notice that (2) is correct but (1) is incorrect even for t=0 (the initial condition!):


Maple seems to prefer the wrong solution (1) but occasionally produces (2) e.g. in a new whorksheet!
In earlier versions it seems that only (1)  appears.

The same ODE with another IC

dsolve({ diff(y(t),t) = y(t)^2 - y(t)^3, y(0) = 1/100 }, y(t));




is always incorrect. It should be


but Lambert's function never shows up!
Let me mention that only the exact solutions are affected, numeric is ok.

Without an initial condition, dsolve always uses LambertW:

dsolve({ diff(y(t),t) = y(t)^2 - y(t)^3}, y(t));



Can you explain this behavior?



Dear Friends

My problem is related to exporting Maple code to Latex, to be more specific I want Latex output that look exactly same as it appears in Maple worksheet like like Red Maple prompt etc.

I have found another solution to this:

"In file menu go to “Export As”, then save file in Tex format. But in order to run this Tex file you need to have files like maple2e.sty, mapleenv.sty, maplestd2e.sty, mapletab.sty, mapleenv.def, mapleplots.sty, maplestyle.sty, mapleutil.sty(all part of package maplestd2e), these all file can be copied from C:drive where maple is installed, just copy all these file to Latex folder and you can run Tex file exported from Maple."

But this do not give exactly what I want.

Please see Maple file and its exported Tex version(Not possible to upload) along with pdf output for Maple worksheet.


Let A and B be two lists of monomials. I want a new list C contained that monomials of A where not divide by any monomial of B. For example if A=[x,y,x^2y,xy^2,y^2] and B=[x^2,y^3] then C=[x,y,xy^2,y^2].

I am getting the following expression when I partially differentiate an expression:

PDE11 := diff(theta(z, p), z, z, p)+2*lambda(p)*theta(z, p)*(diff(lambda(p), p))+lambda(p)^2*(diff(theta(z, p), p))+lambda(p)^2*(sin(theta(z, p))-theta(z, p))+2*p*lambda(p)*(sin(theta(z, p))-theta(z, p))*(diff(lambda(p), p))+p*lambda(p)^2*(cos(theta(z, p))*(diff(theta(z, p), p))-(diff(theta(z, p), p)))

I differentiate the above equation to get each term in the form of :table([f=......])

(table([f = 1+sum(Lambda[n](0)/factorial(n), n = 1 .. infinity)]))(p)^2

It is difficult to understand the expression. Maple does not show any error. Can you please tell me what the error is?

Hi, I am trying to find adomian's polynomial of exp(y), but after execution it shows DF(e^y), D^2F(e^y) as it is. why it does'nt show its derivative?plz help

Hi, I am working on an assignment and I am having difficulties with this Eurler method since I am working with vector.

The assignment:

If the aircraft's engine transmits power P to the air, the force of the propellar, Fp, satisfies P=Fpv. Since the force Fp acts in the direction of motion, it can be written in vector form as: (what I have so far)


v:=Vector(2,1,[vx,vy]); velocity vector

nu:=v->Norm(Vector(v),2); (length of v)

Fd:=-1/2*C*rho*A*nu*v # air resistance

Fl:=1/2*Cl*rho*A*nu*v # liftforce in y direction here v=[-vy,vx]

Fp:=(P/nu^2)*v^2 # power from aircraft

All this adds up to:

F_T:=v->Vector(evalm(Fp/nu(v)^2*v+Fd*nu(v)*v+Fl*nu(v)*zip((x,y)->x*y,v([2,1]),Vector([-1,1]))+Vector(2,1,[0,-m*g]))); # I have tried it and it works


so up til here I am good. It is this next part that I can't seam to get.

Write a routine which uses Euler's method to solve the equation of motion numerically for velocity v(t):


Test the program by trying a plane with intitial velocity [v0,0] such that it maintains this speed when the engines power is P=P0

We have been given an Esolve we can use:

Esolve:=proc(f::procedure,h,x0,y0,N::integer) # Calculates y(x0+n*h) for n=1..N given y0=y(x0) and f(x,y)

local n,y,x;

y[0]:=y0; # Start Value

for n from 1 to N do
y[n]:= y[n-1]+h*f(y[n-1],x); # Euler's formula
x:=x+h; # Next x
end do;


end proc;


If anyone could help me understand how it works and how to get this to work with vectors I would be very gratefull.

How can I plot the p stability region in H1and H2 region I already gotten the polynomial 

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