Why Deplot3d does not work in Maple 2016?

For example, 

with(DEtools);

DEplot3d({diff(y(x), x) = y(x)-z(x), diff(z(x), x) = z(x)-2*y(x)}, {y(x), z(x)}, x = 0 .. 3, [[y(0) = 1.638, z(0) = 2.31]], y = 0 .. 2, z = -4 .. 4, scene = [x, z(x), y(x)], linecolor = COLOR(HUE, .5));

There is nothing shown in Maple 2016.

Let us consider

with(Statistics):
X1 := RandomVariable(Normal(0, 1)):
X2 := RandomVariable(Normal(0, 1)):
X3 := RandomVariable(Uniform(0, 1)): 
X4 := RandomVariable(Uniform(0, 1)):
Z := max(X1, X2, X3, X4); CDF(Z, t);

int((1/2)*(_t0*Heaviside(_t0-1)-_t0*Heaviside(_t0)-Heaviside(1-_t0)*Heaviside(-_t0)+Heaviside(-_t0)+Heaviside(1-_t0)-1)*(1+erf((1/2)*_t0*2^(1/2)))*(2^(1/2)*Heaviside(_t0-1)*exp(-(1/2)*_t0^2)*_t0-2^(1/2)*Heaviside(_t0)*exp(-(1/2)*_t0^2)*_t0-2^(1/2)*Heaviside(-_t0)*Heaviside(1-_t0)*exp(-(1/2)*_t0^2)-Pi^(1/2)*undefined*erf((1/2)*_t0*2^(1/2))*Dirac(_t0)-Pi^(1/2)*undefined*erf((1/2)*_t0*2^(1/2))*Dirac(_t0-1)+2^(1/2)*Heaviside(-_t0)*exp(-(1/2)*_t0^2)+2^(1/2)*Heaviside(1-_t0)*exp(-(1/2)*_t0^2)-Pi^(1/2)*undefined*Dirac(_t0)-Pi^(1/2)*undefined*Dirac(_t0-1)+Pi^(1/2)*Heaviside(_t0-1)*erf((1/2)*_t0*2^(1/2))-Pi^(1/2)*Heaviside(_t0)*erf((1/2)*_t0*2^(1/2))-exp(-(1/2)*_t0^2)*2^(1/2)+Pi^(1/2)*Heaviside(_t0-1)-Pi^(1/2)*Heaviside(_t0))/Pi^(1/2), _t0 = -infinity .. t)

whereas Mma 11 produces the correct piecewise expression (see that here screen15.11.16.docx).

Edit. Mma output.

I'm trying solve a problem with NLPSolve using procedures like constraints, but it say: Error, (in Optimization: -NLPSolve) constraints must be specified as a set or list of procedures. Theses constraints are inequality and equality.  How i can get procedures as a list or set?. Is there anyone to give a example of it?. Thanks!

Hello

I have a procedure that builds an ideal from a specific set of polynomials and then calls the Groebner basis package to eliminate some of the variables.  Even though the procedure is running on a machine with 2 processors, 24 cores and 72 GB of ram, only one core has been used (and is always on a 100% usage).  Would the Grid Computing Toolbox be of some hope in this case?  If so,  how to insert the Grid commands so that Maple sends the calculations to the other cores (I find the document rather confusing)?   If I am talking non sense,  please let me know.

Many thanks

Ed

fieldplot is a wonderful tool for plotting vector fields. The option 'fieldstrength' is very useful to scale the arrows so that one can better visualize the field. I often use fieldstrength=log.

However, if one includes the RealDomain library, the fieldstrength=log option fails. I don't see why it should. Can someone enlighten me before I report this as a bug.

I am having difficulty to solve the following ODE

ode:=diff(R(x),x$6)=Dirac(x-y);

with initial and boundary conditions

bcs := ((D@@3)(R))(1) = 0, R(0)-((D@@5)(R))(0) = 0, ((D@@2)(R))(0)-((D@@3)(R))(0) = 0, ((D@@5)(R))(1) = 0, (D(R))(0)+((D@@4)(R))(0) = 0, ((D@@4)(R))(1) = 0;

I want to solve the above ODE for two cases 1) x>y and 2) x<=y

Thanks

Hello,

I need a bimodal distribution. Since I could not find any among the ones provided by Maple, I created a simple one:

with(Statistics):
U := Distribution(PDF = (proc (t) options operator, arrow; piecewise(t < -5, 0, t < 5, -(1/2000)*t^4+(9/1000)*t^2+7/80, 0) end proc)):
X := RandomVariable(U):

#Plotting PDF and CDF works fine:
plot(PDF(X, t), t = -infinity .. infinity);

plot(CDF(X, t), t = -infinity .. infinity)

However, plotting the quantile function does not work:

plot(Quantile(X, z), z = 0 .. 1);

it has a decreasing part for z<1/2 and a discontinuity at z=1/2.
I can plot it correctly as
plot('Quantile'(X, z), z = 0 .. 1);

but I wonder why the first option does not work for such a simple distribution.

Up to http://www.maplesoft.com/support/help/Maple/view.aspx?path=solve&term=solve

 Let us consider 

solve({x > -Pi, (tan(x)-tan(x)^2)^2-cos(x+4*tan(x)) = -1, x < Pi}, [x]);
                               []

We see the command omits the solution x=0 without any warning. It should be noticed that Mathematica solves it, outputting

{{x -> 0}, {x -> 0}}

and the warning

Solve::incs: Warning: Solve was unable to prove that the solution set found is complete.

One may draw a conclusion on her/his own.

Hey, this is not the I've had this encounter. I want to open this saved document but when I open it and Maple starts up it just hits me with "A problem was encountered while opening the workbook. Database is not opened". How can I get to open it properly and see my math notes?

How can this be prevented?

When I try to upload the file in this message it says "Cant open a null file"

Any help?

Jacob

Hey, I recently switched to Maple 2016, but now whenever I write equations in the text field it calculates it like if it was a math problem it had to solve, but I just want it to be text, nothing else.

Do you know of any way to turn this off? As in Maple 2015 I could write equations in the text field, and when I ran the worksheet through it wouldn't try to calculate text fields, but 2016 calculates the text fields and shows the solution as a blue text... So annoying.

Hi,

I want to solve 2 linear equations in p[1] and p[2] 

eq3 = -(1/8)*(x^2+y^2+((x^2+y^2)^2+2*(omega[1]-omega[2])*(x^2-y^2)+(omega[1]-omega[2])^2)^(1/2)+omega[1]-omega[2])*(x^2*p[1]+2*x*y*p[2]-y^2*p[1]-((x^2+y^2)^2+2*(omega[1]-omega[2])*(x^2-y^2)+(omega[1]-omega[2])^2)^(1/2)*p[1]+omega[1]*p[1]-omega[2]*p[1])/(y^2*x*(omega[1]-omega[2])) - P[1];
eq4 = -(1/8)*(x^2*p[1]+2*x*y*p[2]-y^2*p[1]+((x^2+y^2)^2+2*(omega[1]-omega[2])*(x^2-y^2)+(omega[1]-omega[2])^2)^(1/2)*p[1]+omega[1]*p[1]-omega[2]*p[1])*(x^2+y^2-((x^2+y^2)^2+2*(omega[1]-omega[2])*(x^2-y^2)+(omega[1]-omega[2])^2)^(1/2)+omega[1]-omega[2])/(y^2*(omega[1]-omega[2])*x) - P[2]

solve({eq3,eq4},{p[1],p[2]});

I don't receive any answer. Why?

Hello people in mapleprimes,

I hope that you will give me an answer about the calculation surrounding D.
D((x+y)^2) can be calculated with maple, but D((x+y)^a), where a is a constant, cannot at least directly.
And, I have the code, which Mr. Carl Love kindly gave me at
http://www.mapleprimes.com/questions/218045-Logarithmic-Derivative
, which is

h:= F-> expand(evalindets(D(F)/F, specfunc(D), d-> op(d)*'h'(op(d))));

This works for (x+y)^2, but does not work for (x+y)^a.
And, I think this reason is that the behavior of D to (x+y)^2 is not the same
 as that to (x+y)^a.
And, I want to modify this code.
But, I can't .
So, I am asking this question. How should I modify the above code so that it works for (x+y)^a.
I'm so sorry for a lot of impoliteness.
I will be very glad if you give me an answer.

Take care.

taro

Please help me to solve the system of 1st order singular O.D.E  (see uploaded file)....New_Microsoft_Office_Word_Document.docx
 

Hello,

I tried to define the Poisson bracket in this way

restart:
Ps:=proc(f,g)
local P;
P:=simplify(diff(f,P[1])*diff(g,X)-diff(f,X)*diff(g,P[1])+
diff(f,P[2])*diff(g,Y)-diff(f,Y)*diff(g,P[2]));
end;

The problem is that, if I ask Ps(P[1],X) I got 0 which is clearly wrong.

Where is the error?

https://dl.dropboxusercontent.com/u/51122291/Scrreshot.jpeg

Thanks Nicola 

I'd like to bring up all the Maple 2016 applications from the application center.

How do I do that?