Maple 2015 Questions and Posts

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

There is a desire to explore the process of filtration combustion. To do this, you must solve a system of differential equations in partial derivatives.
I write down all the equations.
Boundary conditions in Maple 2015.0 interpreted incorrectly.
I need to write like that:

 

 

given that:


It turns out so:

 

 

ie somewhere lost derivatives

,....

 

 

As in Maple record boundary conditions correct?

Thanking you in advance.

Hi,

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!

Hi

Anyone could help me in solving the following system of equations to get constants C1, C2, C3 and C4. MALPE give me this "soution may have been lost".  The MAPLE sheet is also attached.

 

restart:

Eq1:=simplify(C3*exp(-(1/4)*(C2*(x^2-2*0)+sqrt(C2*(x^2-2*0)^2+4*M*(x^2-2*0)*w1*(x^2-2*0)))/w1)+C4*exp((1/4)*(-C2*(x^2-2*0)+sqrt(C2*(x^2-2*0)^2+4*M*(x^2-2*0)*w1*(x^2-2*0)))/w1)-U) = 0;

C3*exp(-(1/4)*(C2*x^2+(x^4*(4*M*w1+C2))^(1/2))/w1)+C4*exp(-(1/4)*(C2*x^2-(x^4*(4*M*w1+C2))^(1/2))/w1)-U = 0

(1)

Eq2:=simplify(exp(-(1/4)*(C2+sqrt(C2^2+4*M*w1))*(x^2-2*0)/w1)*C3*x+exp((1/4)*(-C2+sqrt(C2^2+4*M*w1))*(x^2-2*0)/w1)*C4*x+C2-V-z) = 0;

exp(-(1/4)*(C2+(C2^2+4*M*w1)^(1/2))*x^2/w1)*C3*x+exp(-(1/4)*(C2-(C2^2+4*M*w1)^(1/2))*x^2/w1)*C4*x+C2-V-z = 0

(2)

Eq3:=simplify((-2*w2*w5*ln(C3*exp(-(1/2)*sqrt(w2*w4*(w2*w4+w3*w6))*C2*(x^2-2*0)/(w2*w4*w5))-C4)*sqrt(w2*w4*(w2*w4+w3*w6))+w2*w5*(-w2*w4+sqrt(w2*w4*(w2*w4+w3*w6)))*ln(exp(-(1/2)*sqrt(w2*w4*(w2*w4+w3*w6))*C2*(x^2-2*0)/(w2*w4*w5)))+C1*w3*w6*sqrt(w2*w4*(w2*w4+w3*w6)))/(sqrt(w2*w4*(w2*w4+w3*w6))*w3*w6)-1)= 0;

(-ln(exp(-(1/2)*(w2*w4*(w2*w4+w3*w6))^(1/2)*C2*x^2/(w2*w4*w5)))*w2^2*w4*w5+C1*w3*w6*(w2*w4*(w2*w4+w3*w6))^(1/2)-2*w2*w5*ln(C3*exp(-(1/2)*(w2*w4*(w2*w4+w3*w6))^(1/2)*C2*x^2/(w2*w4*w5))-C4)*(w2*w4*(w2*w4+w3*w6))^(1/2)+ln(exp(-(1/2)*(w2*w4*(w2*w4+w3*w6))^(1/2)*C2*x^2/(w2*w4*w5)))*(w2*w4*(w2*w4+w3*w6))^(1/2)*w2*w5-(w2*w4*(w2*w4+w3*w6))^(1/2)*w3*w6)/((w2*w4*(w2*w4+w3*w6))^(1/2)*w3*w6) = 0

(3)

Eq4:= simplify((-C2*x^2*w2*w4-.50*C2*x^2*w3*w6+sqrt(w2*w4*(w2*w4+w3*w6))*C2*x^2+2.*w2*w4*w5*ln(w3^4*w6^2*(C3^2*exp(-1.0*sqrt(w2*w4*(w2*w4+w3*w6))*C2*x^2/(w2*w4*w5))-2*C3*exp(-.5*sqrt(w2^2*w4^2+w2*w3*w4*w6)*C2*x^2/(w2*w4*w5))*C4+C4^2)/(w2*w4*(w2*w4+w3*w6)*C2^2))-5.544000000*w2*w4*w5-w3^2*w6)/(w3^2*w6)) = 0;

(-C2*x^2*w2*w4-.5000000000*C2*x^2*w3*w6+(w2*w4*(w2*w4+w3*w6))^(1/2)*C2*x^2+2.*w2*w4*w5*ln(w3^4*w6^2*(C3^2*exp(-(w2*w4*(w2*w4+w3*w6))^(1/2)*C2*x^2/(w2*w4*w5))-2.*C3*exp(-.5*(w2*w4*(w2*w4+w3*w6))^(1/2)*C2*x^2/(w2*w4*w5))*C4+C4^2)/(w2*w4*(w2*w4+w3*w6)*C2^2))-5.544000000*w2*w4*w5-w3^2*w6)/(w3^2*w6) = 0

(4)

solve({Eq1, Eq2, Eq3,Eq4}, {C1, C2, C3,C4});

Warning, solutions may have been lost

 

``

``

Download solution_lost.mw

 

 

Hello all,

Yesterday, I upgraded to Mac OS X El Capitan.

Now, when working with Maple 2015, I feel the gui is very slow and sometimes irresponsive, when trying to scroll through my worksheets as well as through help pages (e.g., "help plot3d"). When I do the same within Maple 18, it works without any problems.

Does anybody else have the same issue?

Cheers,

Franky

Hello

 

How do I define this function? s(t) = 5*t1/2

This function shows the pace of a particle. I have to decide the time when the pace of the particle is 2 m/s

 

help 

Why does the collect command work for some expressions and not for others. Here is a screen shot

I assume the collect command is supposed to rewrite the expression in terms of the variable descending order.

p := expand((a^2+2*x)*(a^2+2*x));
                        4      2        2   2
                       a  + 4 a  x + 4 x

collect(p, x);
                        4      2        2   2
                       a  + 4 a  x + 4 x

Does not work.

But if you look at the screenshot , it works for other expressions.

Hello

Previously today I had to do some chemistry assignments, which I did in maple, my long time favorite for such assignments.

However I have now encountered a severe problem. I cannot open my worksheet file, .mw. When I try to open it, it asks me whether to open the file as maple text, plain text or maple input??? 

Whatever one I choose, it either opens a blank page in maple, or maple just crashes.

The file in question should be attached. If you can open it, don't mind the language, it is danish

Any help will be greatly appreciated, best regards, Jens

Maple Worksheet - Error

Failed to load the worksheet /maplenet/convert/Opgaver_fra_bog.mw .

Download Opgaver_fra_bog.mw

Here is a serious achievement of the Roots command:


Student[Calculus1]:-Roots(2^x+3^x+6^x-x^2);

                              [-1]

plot(2^x+3^x+6^x-x^2, x = -6 .. 2, gridlines = false);





Download Roots.mw

The solve command also does the job here:

sol := solve(2^x+3^x+6^x-x^2);



allvalues(sol);


evalf(%):

The RealDomain:-solve command fails here.

I wonder how Maple solves it. It would be kind of Maple developers and experts to explain that.

PS. I tried printlevel:=10, but understood the output a little.

 

 

Hi, i got a litlle problem with fromats, i m realy bad at this.
i got procedure 'getone(z::list)' i was execute it like this  
>getone([seq(lineP[i], i = 1 .. pm)]::list)
where lineP is 

>pm := 3;
>rollx := rand(98.0 .. 102.0); rolly := rand(-1.0 .. 5);
>rolll := rand(.5 .. 1.0); rollf := rand(0. .. evalf(2*Pi));
>rollm := rand(0. .. 20.0);

>for i to pm do lineP[i] := [rollx(), rolly(), rolll(), rollf(), rollm()] end do;

it's pm times of z:=[rollx(), rolly(), rolll(), rollf(), rollm()] ::list

and i need make optimization work for any pm=1..100; with rabdomly generated iteration point as above.

The main queston is how make minimize as below
Minimize(z,variables=[z],   initialpoint={z=[seq(lineP[i], i = 1 .. pm)]},iterationlimit=1000,optimalitytolerance=0.01)

work with z::list kind of 

[[100.7798614, 1.109653266, 0.9840371500, 4.257816686, 9.737110573][100.0135459, 0.887539481, 0.9144164697, 3.980093624, 7.343851161][100.0661308, 3.724268229, 0.5020544909, 2.052134822, 14.37408543]]








I find fsolve() to be curiously unstable. Some of the behaviour I can guess at, but other parts are not as clear.

 

Define a function.

P := x->fsolve(exp(-(1/100)*t^2)*cos(2*t),t=x):

plot(P,4..5,smartview=false);

I set up a worksheet with embedded components in Maple 18. I am now trying to use the worksheet with Maple 2015. The GetProperty( ) procedure in DocumentTools cannot find a number of kinds of the embedded components (perhaps all kinds).

Why would this happen, and how can I fix it?

Hello

I have a simple list:

xlist := [150, 250, 500, 800, 1300, 2500, 5000]:

ylist := [.3, .5, .8, 1.0, 1.2, 1.4, 1.6]:

where after i said x:=plot(xlist,ylist);

What I want to do is ask maple" At what x value does the graph intersect with 1.5 on the y axis" or "at what y value does the graph intersect with 3000 on the x axis". 

The data points are just that, I've been looking for commands to ask these questions but have been without luck. You don't need to give me the answer if you can point me to somewhere where this information is written, that would be very helpful!

 

Thank you.

Hi everybody,

i'm trying to do an elliptic regularization but i don't know how to proceed ?

Is someone know how to achieve to do that  with an example ?

thanks a lot !

 

PS: i know only how to do a linear regularisation.

 

 

I am considering the following PDE and I am getting an error, please suggest a better numerical method than the default one used in maple:

 

the PDE is:

u_{xx}u^3 - sin(xt)u_{tt} = u(x,t)

u(x, 0) = sin(x), (D[2](u))(x, 0) = cos(x), u(0, t) = cos(t), (D[1](u))(0, t) = sin(t)

Please suggest me a method that will also work for the following PDEs:

u^m* u_{xx} - sin(xt)u_{tt} = u^n

for m,n =0,1,2,3,... for the cases m=n and m not equal n

Here's the code:

 

pde := u(x, t)^3*(diff(u(x, t), x, x))-sin(x*t)*(diff(u(x, t), t, t)) = u(x, t);

u(x, t)^3*(diff(diff(u(x, t), x), x))-sin(x*t)*(diff(diff(u(x, t), t), t)) = u(x, t)

(1)

ibc := u(x, 0) = sin(x), (D[2](u))(x, 0) = cos(x), u(0, t) = cos(t), (D[1](u))(0, t) = sin(t);

u(x, 0) = sin(x), (D[2](u))(x, 0) = cos(x), u(0, t) = cos(t), (D[1](u))(0, t) = sin(t)

(2)

pds := pdsolve(pde, [ibc], numeric, time = t, range = 0 .. 1, spacestep = 0.1e-1)

module () local INFO; export plot, plot3d, animate, value, settings; option `Copyright (c) 2001 by Waterloo Maple Inc. All rights reserved.`; end module

(3)

pds:-plot3d(u(x, t), t = 0 .. 1, x = 0 .. 1, labels = [t, x, u(x, t)], labelfont = [times, bold, 20], axesfont = [times, bold, 16])

Error, (in pdsolve/numeric/plot3d) unable to compute solution for t>HFloat(0.0):
Newton iteration is not converging

 

``

 

Download nonlinear_hyperbolic_PDE.mw

First 11 12 13 14 15 16 17 Last Page 13 of 71