Maple 2016 Questions and Posts

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

I have the following PDE system steaming from Flash Photolysis:

pdesys := [diff(J(x, t), x) = varepsilon*J(x, t)*C(x, t), diff(C(x, t), t) = phi*varepsilon*J(x, t)*C(x, t)]

when I use pdsolve(pdesys,[J,C]) I get:

{C(x, t) = 0, J(x, t) = _F1(t)}, {C(x, t) = _F1(x)*_C1, J(x, t) = 0}

The solution appears to be either C(x,t) = 0 or J(x,t) = 0. This are the obvious solutions (0 = 0). I have the analytical solution to this PDE system where neither C(x,t) nor J(x,t) are 0.

How to solve this system in maple? Thanks.


Example 6 from Maple Help.




LPSolve(2*x+5*y, {3*x-y = 1, x-y <= 5}, assume = {integer, nonnegative})


Kernel crash the same with Maple 2015.


Does anyone can confirm?

Mariusz Iwaniuk

Maple evaluates:

but when I ask:

Maple will verify sin(u)/cos(u)=tan(u) but beyond that it does not seem to work too well. I suspect this has something to do with verify being a "semi-boolean" function but searching the docs for what that mean returns no results. It may mean "I am unable to decidee". But how should I know? :-(

This is Maple 2016.1 on Windows. There is a post on verifying very simple identities.

Edit As usual some very helpful answers for which many thanks. It is poor that Maple does not document the meaning of FAIL.It seems that for trig identities the expandverfier is the best.

Edit Maple documents FAIL under ? FAIL. I am not sure how I got a blank page when I searched. convert root of to explicit form.

w is a imaginary..


restart; w := (1/2)*(2*d-5+I*sqrt(4*d-9))/(d-2)




{Q1 = RootOf((2*I)*(4*d-9)^(1/2)*_Z*d-((2*d-5+I*(4*d-9)^(1/2))/(d-2).(d*(I*(4*d-9)^(1/2)+1)*_Z/((d-2)*b)))*b*d+2*((2*d-5+I*(4*d-9)^(1/2))/(d-2).(d*(I*(4*d-9)^(1/2)+1)*_Z/((d-2)*b)))*b-2*d*_Z), Q2 = -(1/2)*d*(I*(4*d-9)^(1/2)+1)*RootOf((2*I)*(4*d-9)^(1/2)*_Z*d-((2*d-5+I*(4*d-9)^(1/2))/(d-2).(d*(I*(4*d-9)^(1/2)+1)*_Z/((d-2)*b)))*b*d+2*((2*d-5+I*(4*d-9)^(1/2))/(d-2).(d*(I*(4*d-9)^(1/2)+1)*_Z/((d-2)*b)))*b-2*d*_Z)/((d-2)*b)}







I'm new to the physics package - wondering if i can tweak it a bit to look like things i'm used to:


is there a way to make Christoffel symbols print as upper case gamma, instead of  'G'?

KroneckerDelta print as lower case delta, instead of 'd'?


can i make the Schwarzschild metric look like it does in Hartle, Carroll, and others:


-(1 - 2M/r)dt^2 + (1-2M/r)^-1 + r^2(dtheta^2 + sin(theta)^2 dphi^2)


i know about setting the signature in Setup.

i have tried the 'Coordinates' command, but when i give it X=[t,r,theta,phi] i always seem to get back



i am running maple 2016


many thanks,



Hello everybody.

I have a function:

f(x,y)=GAMMA(y, -ln(x))/GAMMA(y)

seq(sum(f(x, y), y = 0 .. 1), x = 0 .. 5)


and I got a error message:

Error, (in ln) numeric exception: division by zero ??
This is normal behavior in seq function or Bug?


but  when I'm first calculate the sum sol := sum(f(x, y), y = 0 .. 1) -> x,

and evalf([seq(sol, x = 0 .. 5)]) ->[0., 1., 2., 3., 4., 5.] works fine.

Mariusz Iwaniuk

My old harddisk recently died and I am currently using a new one. I still have the license code for Maple 2016. Is it still possible for me to install Maple 2016 or 1 license = 1 installation?

Maplet file failed to open by double click in windows 7 x64 file browser

There is no response when double click a Maple file

But it runs correctly by using command line

Anyone knows why?

I attempted to show that two lines are parallel.  I started with a problem in Geometry for which I do not have the solution.

I tried several ways with Maple to show this to be true.  Most of the time, I ended when maple could not determine if a-b = c-d, etc.

brg_proof.txt contains a statement of the problem and my latest maple code.

Question: How should I approach the proof, by the compass and straight edge method?  Is this possible in maple?

intersection in the geometry package does not seem to recognize assume.

restart: with(geometry):

assume(p[1]<>0, p[2]<>0, p[3]<>0);
assume(q[1]<>0, q[2]<>0, q[3]<>0);



point(A,eval(x,sol),eval(y,sol));  ## the intersection exists
intersection(xA,AD,OT); ## fails

I have several plots generated with geometry objects:

plt1 := draw(...); plt2 := draw(...)

I want to display all them with a reasonable window.  I have been looping and and using a rather dumb proc to find the window size. I use the view=[...] option with the values calculated to set the plot view.

I wanted to post the proc, but, really messed that up!  I will try the proc later.


Some of the plots I am working with have views like

> for itm in [relBrgLinePlot,relTgtPositPlot,relTgtLinePlot] do
end do;
                     VIEW(0. .. 95., -61.58179461 .. 16.)
                 VIEW(5.48225506 .. 95., -61.58179461 .. 16.)
                 VIEW(5.48225506 .. 95., -61.58179461 .. 16.)

> for itm in [geoSensPositPlot,geoSensLinePlot,\
            geoBrgLinePlot,geoTgtPositPlot,geoTgtLinePLot] do
end do;
                  VIEW(0. .. 35.35331852, 0. .. 87.13244438)
                  VIEW(0. .. 35.35331852, 0. .. 87.13244438)
                     VIEW(0. .. 95., 0. .. 87.13244438)
                 VIEW(40.83557358 .. 95., 16. .. 25.55064977)
                 VIEW(40.83557358 .. 95., 16. .. 25.55064977)

Is there an easier way to do this?

Linux.  I want to put grid lines on an x11 device.

The x11 device is better for plotting because if the window is resized, the plot resizes with it.

plotsetup(x11);plot(sin(x),x=-1..1,gridlines = true); ## no gridlines

plot(sin(x),x=-1..1,gridlines = true); ## have gridlines

Any way to get gridlines on x11 device?

Tom Dean

I have a Maple result presented in the first (top) formula which I would like to format as in the second (bottom) formula:

Can this be done? This is Maple 2016.1


I have experienced that maple does not save all of the varibles. But some it does.

I calculate with units, could that be the reason?

I have allso been thinking that it has something to do with saving the document online in onenote. But that works like the file is saved on the Pc's harddrive.

Are there anybody else that has experienced this?

I calculate with units, but as the varible does not appear in the calculation with units, I make one varible with the same result, to get the next calculation to work.   



Why won't this procedure Compile? All I get is a cryptic and ungrammatical error message.

    step:= proc(
          XYZ::Matrix(datatype= float[8]),
          E::Vector(datatype=integer[4]), F::Vector(datatype=integer[4]),
          W::Matrix(datatype= float[8]), #3x2 scratch matrix
     option autocompile;
          i::integer[4], j::integer[4],
          p::float[8], t::float[8]       
          to mu do
               for i to n do
                    ed:= 0;  fd:= 0;
                    for j to 3 do
                          p:= XYZ[i,j];
                          t:= XYZ[E[i],j] - p;
                          ed:= ed+t^2;
                          W[j,1]:= t;
                          t:= XYZ[F[i],j] - p;
                          fd:= fd+t^2;
                          W[j,2]:= t
                    ed:= sqrt(ed) + .01;  fd:= sqrt(fd) + .01;
                    for j to 3 do
                         XYZ[i,j]:= 0.995*XYZ[i,j] - 0.01*W[j,1]/ed + 0.02*W[j,2]/fd
     end proc;

Error, (in Compiler:-Compile1) In memory compilation failed

First 28 29 30 31 32 33 34 Page 30 of 35