Maple 2016 Questions and Posts

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

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

• 

If the solve command does not find any solutions, then if the second argument is a name or set of names, then the empty sequence (NULL) is returned; if the second argument is a list, then the empty list is returned. This means that there are no solutions, or the solve command cannot find the solutions. In the second case, a warning is issued, and the global variable_SolutionsMayBeLost is set to true.

 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?

hi

please help me for fsolve algebric equations...

thanks

AGM.mw
 

restart:

F(eta):=sum(a[i]*eta^i,i=0..5):

theta(eta):=sum(b[i]*eta^i,i=0..5):K(eta):=sum(c[i]*eta^i,i=0..5):Omega(eta):=sum(d[i]*eta^i,i=0..5):

``

U1:=diff(theta(eta), eta, eta)-3*Omega(eta)*(F(eta)*(diff(theta(eta), eta))-theta(eta)*(diff(F(eta), eta)))/(2*K(eta))+((diff(K(eta), eta))/K(eta)-(diff(Omega(eta), eta))/Omega(eta))*(diff(theta(eta), eta)) = 0:U2:= diff(F(eta), eta, eta, eta)+Omega(eta)*(3*F(eta)*(diff(F(eta), eta, eta))-(diff(F(eta), eta))^2)/(2*K(eta))+((diff(K(eta), eta))/K(eta)-(diff(Omega(eta), eta))/Omega(eta))*(diff(F(eta), eta, eta))+Omega(eta)/K(eta) = 0:U3:= diff(K(eta), eta, eta)+Omega(eta)*(1.5*F(eta)*(diff(K(eta), eta))-K(eta)*(diff(F(eta), eta)))/K(eta)+((diff(K(eta), eta))/K(eta)-(diff(Omega(eta), eta))/Omega(eta))*(diff(K(eta), eta))+(diff(F(eta), eta, eta))^2-Omega(eta)^2 = 0:U4:= diff(Omega(eta), eta, eta)+Omega(eta)*(3*F(eta)*(diff(Omega(eta), eta))+Omega(eta)*(diff(F(eta), eta)))/(2*K(eta))+((diff(K(eta), eta))/K(eta)-(diff(Omega(eta), eta))/Omega(eta))*(diff(Omega(eta), eta))+Omega(eta)*(diff(F(eta), eta, eta))^2/K(eta)-Omega(eta)^3/K(eta) = 0:

F(eta):=unapply(F(eta),eta):

theta(eta):=unapply(theta(eta),eta):K(eta):=unapply(K(eta),eta):Omega(eta):=unapply(Omega(eta),eta):

U1:=unapply(U1,eta):U2:=unapply(U2,eta):U3:=unapply(U3,eta):U4:=unapply(U4,eta):

 

s1:=F(eta)(0) = 0:

s2:=K(eta)(0) = 0:s3:=Omega(eta)(0) = 0:s4:=theta(eta)(0) = 1:s5:=theta(eta)(1) = 0:s6:=(D(F(eta)))(0) = 0:s7:=(D(K(eta)))(1) = 0:s8:=(D(Omega(eta)))(1) = 0:s9:=((D@@2)(F(eta)))(1) = 0

20*a[5]+12*a[4]+6*a[3]+2*a[2] = 0

(1)

s10:=U1(0):s11:=U2(0):s12:=U3(0):s13:=U4(0):        s14:=U1(1):s15:=U2(1):s16:=U3(1):s17:=U4(1):    s18:=D(U1)(0):s19:=D(U2)(0):s20:=D(U3)(0):s21:=D(U4)(0):     s22:=D(U1)(1):s23:=D(U2)(1):s24:=D(U3)(1):

 

 

Q:=fsolve([s1,s2,s3,s4,s5,s6,s7,s8,s9,s10,s11,s12,s13,s14,s15,s16,s17,s18,s19,s20,s21,s22,s23,s24],{a[0],a[1],a[2],a[3],a[4],a[5],b[0],b[1],b[2],b[3],b[4],b[5],c[0],c[1],c[2],c[3],c[4],c[5],d[0],d[1],d[2],d[3],d[4],d[5]}):

F(eta):=eval(sum(a[i]*eta^i,i=0..5),Q):

Error, invalid input: eval received S, which is not valid for its 2nd argument, eqns

 

theta(eta):=eval(sum(b[i]*eta^i,i=0..5),Q):K(eta):=eval(sum(c[i]*eta^i,i=0..5),Q):Omega(eta):=eval(sum(d[i]*eta^i,i=0..5),Q):

Error, invalid input: eval received S, which is not valid for its 2nd argument, eqns

 

plot(g(x),x=0..1,axes=boxed,color=green,thickness=2,labels=[x,g]):

Warning, unable to evaluate the function to numeric values in the region; see the plotting command's help page to ensure the calling sequence is correct

 

plot(f(x),x=0..1,axes=boxed,color=blue,thickness=2,labels=[x,f]):

Warning, unable to evaluate the function to numeric values in the region; see the plotting command's help page to ensure the calling sequence is correct

 

 


 

Download AGM.mw

 

How to plot f(x)=x in the interval -2<x<2 with a period 4?

The plot should look like this


Hello all, first time posting on here so if I make any mistakes please let me know. So I have created a function called "g" and plotted it on a graph. After looking at the two graphs I have plotted (one with discont=false and one with discont=true), I can see 3 discontinuities are present. However, attempting to solve for these discontinuities by setting the denominator of "g" = 0 only returns 1 of the points. I'm not sure why this is happening and how would I go by solving for the other 2 points? Thank you.

f := exp(x-188.5)-4*x^2+1507.964*x-142122.3+cos(x)^2;

exp(x-188.5)-4*x^2+1507.964*x-142122.3+cos(x)^2

(1)

g := (-f+x^2-376.991*x+35530.6)/(f+x^2-376.991*x+35530.6);

(-exp(x-188.5)+5*x^2-1884.955*x+177652.9-cos(x)^2)/(exp(x-188.5)-3*x^2+1130.973*x-106591.7+cos(x)^2)

(2)

``

plot(g, x = 186 .. 194, y = -25 .. 25,discont=false)

 

plot(g, x = 186 .. 194, y = -25 .. 25,discont=true)

 

vasymp := fsolve( denom(g)=0);

192.2026966

(3)

``

After changing the time and date on my computer, I get a message saying “Invalid or missing license file.” whenever I try to open Maple. I have tried activating it with my license code. It then gives me a message saying “Activation succesful! Please restart Maple for the new license file to take effect.”. After opening it again, I will once again get the original message. I have also tried reinstalling it twice. One time with the file my school provides, and with a file from Maples website. I get the same error message no matter what I do. 

I haven't been able to access Maple for several weeks now, so I'm getting a bit behind in math class.

How should I go about fixing this?

Hi all,
I have a question about editing strings in Maple2016.
When I make a variable "DataName" which is as below in Maple worksheet,
DataName:=Array(["Data1","Data2","Data3"]):

I want to make the variable "DataSet" whose elements are as below.
DataSet:=[1="Data1",2="Data2",3="Data3"]

Of course I can type the elements of "DataSet", but if the number of DataName increases, it requires a lot of time to do that, so I want to make the "DataSet" by using general command so that I can obtain "DataSet" automatically.
I tried to use the command "cat", "||", "map", and so on... but I could not find how to do that.
Is there anyone who can solve this problem? Please let me know how to solve this problem. I appreciate any idea you may have.

Always thanks for the nice answer

First 42 43 44 45 46 47 48 Last Page 44 of 61