Maple 2015 Questions and Posts

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

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.

Elast2.mw

 




 

 

Here the potential of maple 2015 to the quantitative study of the decomposition of a vector table is shown in two dimensions. Application for the exclusive use of engineering students, which was implemented with embedded components.

Atte.

Lenin Araujo Castillo

Archivo Corregido:  Decomposición_Vectorial_Corregido.mw

Hello!

 

Is it possible convert some equations in R² (ellipse, hyperbola) to polar coordinates using Maple? And a regular  object in R³ like a sphere, a cylinder or a cube could be converted from cartesian coordinates to polar coordinates?

 

Thank you so much.

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

HI there,

I m getting an error message . Could someone help me

 

v1 := int(cos(tau)*g(tau), tau = t0 .. t);
int(cos(tau) g(tau), tau = t0 .. t)
v2 := int(-sin(tau)*g(tau), tau = t0 .. t);
int(-sin(tau) g(tau), tau = t0 .. t)
soln := C1*y1+C2*y2+v1*y1+v2*y2;
C1 y1 + C2 y2 + (int(cos(tau) g(tau), tau = t0 .. t)) y1

+ (int(-sin(tau) g(tau), tau = t0 .. t)) y2
soln := combine(soln);
(int(-y2 sin(tau) g(tau) + y1 cos(tau) g(tau), tau = t0 .. t))

+ C1 y1 + C2 y2
eval(soln, t = t0) = 0, eval(diff(soln, t), t = t0);
C1 y1 + C2 y2 = 0, -y2 sin(t0) g(t0) + y1 cos(t0) g(t0)
solve({%}, {C1, C2});
soln := eval(soln, %);
Error, invalid input: eval received (C1*y1+C2*y2 = 0, -y2*sin(t0)*g(t0)+y1*cos(t0)*g(t0)), which is not valid for its 2nd argument, eqns

 

Is there a maple function that'll allow me to generate the matrix M(5 by 5) below given that matrix A and B are 3 by 3 M:=LinearAlgebra:-ZeroMatrix(5) A:=Matrix([[a11,a12,a13],[a21,a22,a33],[a31,a32,a33]]) B:=Matrix([[b11,b12,b13],[b21,b22,b33],[b31,b32,b33]]) Matrix(5, 5, {(1, 1) = a[11], (1, 2) = a[12], (1, 3) = a[13], (1, 4) = 0, (1, 5) = 0, (2, 1) = a[12], (2, 2) = a[22], (2, 3) = a[23], (2, 4) = 0, (2, 5) = 0, (3, 1) = a[13], (3, 2) = a[23], (3, 3) = a[33]+b[11], (3, 4) = b[12], (3, 5) = b[13], (4, 1) = 0, (4, 2) = 0, (4, 3) = b[12], (4, 4) = b[22], (4, 5) = b[23], (5, 1) = 0, (5, 2) = 0, (5, 3) = b[13], (5, 4) = b[23], (5, 5) = b[33]});

Just purchased Maple 2015 and playing with it for the first time.  I'm running the 64-bit version on Win 8.1.

Anytime the program generates a pop-up dialog box, the pop-up seems to get stuck behind the main program window.  I can't alt+tab to get to the pop-up window and I can't click on anything in the main program window because it's frozen while waiting for me to aknowledge the pop-up dialog box.  So I have to kill everything from the task manager and lose anything that's unsaved.

I'm talking about pop-ups for things like "Error, (in @@) invalid arguments"... I click on the link and see a brief flash while that pop-up quickly flies behind the main window.  Same thing happens when I click on a link that asks what web browser I want to use.

 

I can't be the only one that has this problem, so is there a fix or workaround for this somewhere?  It's not really useable like this.  Thanks!

I want to solve numerically the PDE:

u_xx + u_yy= = u^{1/2}+(u_x)^2/(u)^{3/2}

 

My assumptions are that  |sqrt(2)u_x/u|<<1 (but I cannot neglect the first term since its in my first order approximation of another PDE.

 

So I tried solving by using pdsolve in maple, but to no cigar.

 

Here's the maple file:

 nonlinear.mw

PDE := diff(diff(u(x, y), x), x)+diff(diff(u(x, y), y), y) = u^(1/2)+(diff(u(x, y), x))^2/u^(3/2); IBC := {D[1](u)*(1, t) = 0, D[2](u)*(x, 1) = 0, u(0, t) = 1, u(x, 0) = 1}; pds := pdsolve(PDE, IBC, type = numeric); pds:-plot3d(t = 0 .. 1, x = 0 .. 1, axes = boxed, orientation = [-120, 40], color = [0, 0, u])

diff(diff(u(x, y), x), x)+diff(diff(u(x, y), y), y) = u^(1/2)+(diff(u(x, y), x))^2/u^(3/2)

 

{D[1](u)*(1, t) = 0, D[2](u)*(x, 1) = 0, u(0, t) = 1, u(x, 0) = 1}

 

Error, (in pdsolve/numeric/process_PDEs) all dependent variables in PDE must have dependencies explicitly declared, got {u}

 

Error, `pds` does not evaluate to a module

 

``

 

Download nonlinear.mw

Hello All,

I looked through the Maple help on PDE systems and pdsolve and the physics problems that appear there. THere are a number of single-PDE cases with initial / boundary conditions; but I couldn't find PDE systems with ics/bcs.

 

Would you have a (simple) example of a PDE system with its initial / boundary conditions? I am attempting to build understanding of the syntax and different options of "pdsolve". Examples seem to be a great way to learn how to solve PDE systems. One can then pdsolve the PDE system without ics; then add them, try different options etc.

 

Thank you!

 

First 64 65 66 67 68 69 70 Page 66 of 72