Maple 2016 Questions and Posts

These are Posts and Questions associated with the product, Maple 2016

What can i do in this case? I,m traying solve a system of non linear equation like this, i want to know Rs and Rsh:

I have the following fuction in Laplace domain,



Where, s1, s2, s3 are the Laplace variables.

x1:=invlaplace(u, s1, x);

This worked. But the next two doesn't work.

y1:=invlaplace(x1, s2, y);

uu:=invlaplace(y1, s3,t);

Even, I tried to plot the unevaluated invlaplace but no luck.

plot3d(uu,x =1..2, y=1..2);

Am I missing something?




I have an ODE which is based on a seperate function, and I would like to make a plot with the information

dsolve([diff(X(W), W) = (0.536000000000000e-3*(1-X(W)))*(1+X(W)), X(0) = 0], numeric)


C_A:= C_A0*(1-X(W))*(1+X(W))

which has been used as part of the ODE.

I would really like to plot C_A as a function of W. I have no problem plotting X as a function to W using odeplot. Ideally I would like to plot C_A and X vs W in the same plot.


I have a document containing notes from my past semester, which will not open properly. When i try to open up the document maple asks if i want to retrieve the content as plain text, 2D-math or maple input. The document is written in document mode and contains both text and 2D math. When I try to open it in plain text, maple gives me nothing, and with the two other options, maple is just loading forever. Can anybody help me retrieve the content og fix the document?

I am working with Euler's eqns. of motion. So want to animate the solution. I could do a simple set of spheres along the x, y, z axes  that represent the body or use an imported STL file.  Looking for some guidance on how to approech this. I can import an STL ok. How do I give it an xyz set of axes and then rotate/move these in an inertial frame? I have the equations and solutions already. A reasonable example would be adding an STL of a pendulum to  to Preben's animation response to a recent post  I placed.


how i can save results on text file from pdsolve?



PDE := diff(u(y, t), t) = diff(u(y, t), y, y)-0.9e-1*(diff(u(y, t), y, y, y, y))

diff(u(y, t), t) = diff(diff(u(y, t), y), y)-0.9e-1*(diff(diff(diff(diff(u(y, t), y), y), y), y))


IBC := {(D[1, 1](u))(0, t)-0.9e-1*(D[1, 1, 1, 1](u))(0, t) = 0, u(0, t) = 1, u(100, t) = 0, u(y, 0) = 0, (D[1, 1](u))(100, t) = 0}

{(D[1, 1](u))(0, t)-0.9e-1*(D[1, 1, 1, 1](u))(0, t) = 0, u(0, t) = 1, u(100, t) = 0, u(y, 0) = 0, (D[1, 1](u))(100, t) = 0}








pds := pdsolve(PDE, IBC, numeric, u(y, t), spacestep = 1/5, timestep = 1/5)











fname1 := "C:/Users/test/Desktop/txtop1.txt"; seq(fprintf(fname1, "%12.8f, %12.8f\n", y, u(y, .1), y = 0 .. 10, 1)); fclose(fname1)

Error, (in fprintf) file or directory does not exist


p1 := pds:-plot(t = 0):p2 := pds:-plot(t = 1/10):

plots[display]({p2}, title = `\` profile at t=0,0.1`)






Hello everyone, I use physics package and I input:


Setup(mathematicalnotation = true)


Setup(tensors = A[mu](X))


PDEtools:-declare(p0(X), pm(X), pp(X), pt(X), Un(X), Ucn(X))

pt := proc (X) options operator, arrow; Matrix(2, 2, [[p0(X), sqrt(2)*pp(X)], [sqrt(2)*pm(X), -p0(X)]]) end proc

U := proc (X) options operator, arrow; Matrix(2, 2, [[1, 0], [0, 1]])+I*pt(X)*(1/f)-pt(X)*pt(X)*((1/2)/f^2) end proc

Ucn := proc (X) options operator, arrow; Matrix(2, 2, [[1, 0], [0, 1]])-I*pt(X)*(1/f)-pt(X) . pt(X)*((1/2)/f^2) end proc

And after that I write:

Trace(-LeviCivita[mu, nu, rho, sigma] . Ucn(X) . ((1/6)*e*A[nu](X)+M5*KroneckerDelta[nu, 0]) . Matrix(2, 2, [[1, 0], [0, 1]]) . Uсn(X) . d_[rho](Un(X)) . Uсn(X) . d_[sigma](Un(X)))

But I get an error:

Error, (in unknown) invalid subscript selector

Thanks for the help.



I have a Maple code which produces an output plot for a first order differential equation,

nde := evalf(subs(npar, de));
nds := dsolve({nde, sigma(0) = -1e-8}, sigma(t), type = numeric);
acc_nds := (sig0, ae, re) -> dsolve({nde, sigma(0) = sig0}, sigma(t), type = numeric,
method = lsode[backfull], abserr = ae, relerr = re, maxfun = 0, ctrl=Ctrl);

odeplot(acc_nds(-0e-7, 1e-13, 1e-13), [t, sigma(t)], t = 0..2);

This produces the outputplot that I need for sigma(t). I need to produce a outut for d(sigmat)/dt, and how can this be done? what is the command I should use?

Additionally, how can i get the data set out of the polt in to a excel file or a text file?

I am quite new to maple, so i expect your kind support


Hello! I try to make friends physics package and matrixes. But I am faced with difficulties. To be more specific, to consider a code:


Setup(mathematicalnotation = true)


Setup(tensors = A[mu](X))

PDEtools:-declare(A(X), p0(X), pm(X), pp(X), pt(X), U(X), m5(X))

And then I create a matrix:

pt := proc (X) options operator, arrow; Matrix(2, 2, [[p0(X), sqrt(2)*pp(X)], [sqrt(2)*pm(X), -p0(X)]]) end proc;


It displays as:

Why do p0, pp, pm become function only one variable x1?

Interestingly that maple understands the next matrix:

Nevertheless even for function only one variable derivative works fine:

And I can construct more complicated functions:

But why is only one variable?



I'm wondering how to produce the graphic in the's_sequence.

my attempt:

I will try to explain with an example, lets say i have to solve 2 equations with 2 unknowns

P := solve({3*y+x, 2*x+3+y})
solve({3*y+x, 2*x+3+y})

I would then like to do P[1] and P[2] to get the values, but instead it gives me "x=-9/5" which make me unable to use it for things like

P[1] + P[2]

Hope you all understand, and thanks for taking your time to read my problem 

Hi, I'd like to know if is it possible to import questions from Maple TA into Maple and use them in quizes, for example. Your help will be highly appreciated. Thank you.

i have these commands:

f := A-B+x^y/y;

g := x^y-y;
Solve({f = 0, g = 0}, {x, y});
              / /                     y    \         \
              | | y                  x     |         |
         Solve|< x  - y = 0, A - B + -- = 0 >, {x, y}|
              | |                    y     |         |
              \ \                          /         /

restart; with(plots); with(VectorCalculus); with(LinearAlgebra); f := A-B+x^y/y;  g := x^y-y; Solve({f = 0, g = 0}, {x, y});


but I can't get it to actually solve the equations. what do I do  ?

How I can solve it for P?

P=i^(2)r+(&DifferentialD;)/(&DifferentialD; t) (1/(2)l(tetha)i^(2))+1/(2)i^(2)(&DifferentialD;l(tetha))/(&DifferentialD; theta) w

attach i(t) "corrente", l(t) "induttanza", theta "angular", w "rotary speed"graph of funtions


Hi! I have a variable polynomial expression and I want to cut off all the terms of order 3 and more, for every product of variables.

For example consider the polynomial P = x + y + a*x^2 + b*x^3 + c*x*y + d*x^2*y + e*x*y^2 + y^2, I want an operation who returns me a*x^2 + c*x*y + y^2

(terms for which the sum of exponents in x and y exceeds or less two are being deleted)

Is it possible?


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