Maple 2016 Questions and Posts

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

Hi all....

I wanted to run the bat file from maple worksheet, if anyone knows please help me.

Thanks

Experts.

In view of the worksheet:

VRP_permute.mw

What I need is
a) a better way to display the Tour_Distances (from 1 to M) and
b) a way to select tour elements where each component is <= to a specified distance (say <=H) and the corresponding tour(s) (from Tour2).
 

My question is: Use the laplace transform to solve the system.

dx/dt + d^2y/dt^2 = 5e^(2t)

dx/dt - x - dy/dt + y = 8e^(2t)

x(0) = 2, y(0) = 1, y'(0) = 1

What I've done in Maple:

with(inttrans);
with(DEtools);
eq5 := (diff(x(t), t)+diff(y(t), t$2) = 5*exp(2*t), t, s);

eq5s := laplace(%, t, s);

eq6 := (diff(x(t), t)-x-(diff(y(t), t))+y = 8*exp(2*t), t, s);

eq6s := laplace(%, t, s);

solve({eq5s, eq6s}, {laplace(x(t), t, s), laplace(y(t), t, s)});

subs({x(0) = 2, y(0) = 1, (D(y))(0) = 1}, %);

eq3 := invlaplace(%, s, t);

How do I simplify?  If you plug it into maple I come up with an answer that has x and y on each side.  I guess I'm just wondering how I can set them equal to each other to solve and get rid of the variable x and y.  I know answer is correct as I've also ran it through ODEtest.  Please help.

However you figure out getting rid of the variables I assume will help me also in solving the next problem:

Use the Laplace Transform to solve the system

dx/dt = 7x - y + 6z

dy/dt = -10x + 4y - 12z

dz/dt = -2x + y - z

x(0) = 5, y(0) = 7, z(0) = 2

I have attempted the second problem much like the first.  Thank you for your time.

Experts,

This may sound like a dumb question, but i'm seeking a procedure to do it better.
 

``

 

with(combinat, setpartition) :
P := [$2..5] :

Tours := setpartition(P);M:=nops(Tours)

[[[5], [2, 3, 4]], [[2], [5], [3, 4]], [[3], [5], [2, 4]], [[4], [5], [2, 3]], [[2], [3], [4], [5]], [[2, 3, 4, 5]], [[2, 5], [3, 4]], [[2], [3, 4, 5]], [[2, 4], [3, 5]], [[3], [2, 4, 5]], [[2, 3], [4, 5]], [[4], [2, 3, 5]], [[3], [4], [2, 5]], [[2], [4], [3, 5]], [[2], [3], [4, 5]]]

 

15

(1)

 

number of elements in each 'group'

seq(nops(Tours[i]),i=1..nops(Tours))

2, 3, 3, 3, 4, 1, 2, 2, 2, 2, 2, 2, 3, 3, 3

(2)

 

i need to add 1 to each 'subgroup' : These are the first two:

[[[1,op(Tours[1,1])],[1,op(Tours[1,2])]],[[1,op(Tours[2,1])],[1,op(Tours[2,2])],[1,op(Tours[2,3])]]]

[[[1, 5], [1, 2, 3, 4]], [[1, 2], [1, 5], [1, 3, 4]]]

(3)

 

I need to add 1 to each 'subgroup' in a more automatic way.

``


 

Download add_1.mw

 

hi
I want to solve a pde equation:
 

equa1 := diff(u(x,y), x, x)-y(1+x) = 0;

# with codition:

con:=u(0,y) = 0, (D(u[x]))(0,y) = 0;

the anwer must be :    u(x,y)= y(x2/2  + x3/6)
How can i solve that with maple?

Please excuse my bad English
thanks

I am plotting a pair of lists of points using

pointplot(Listap, Listbp, symbol = point, symbolsize = 1, size = [1200, 1200])

 

How could is do this with plot so I can add colours? Along the lines Listap(i)^2+Listbp(i)^2 =R, R is in the range 0..1,then colour =R*256 or any other imaginative way of adding colour.
 

restart

with(plots):

with(plottools):

``

``

NULL

``

z := (m+I*n)/(p+I*q)

(m+I*n)/(p+I*q)

(1)

g := proc (z) options operator, arrow; (z-I)/(z+I) end proc;

proc (z) options operator, arrow; (z-I)/(z+I) end proc

(2)

bz := simplify(evalc(Im(z)));

(-m*q+n*p)/(p^2+q^2)

(3)

a := simplify(evalc(Re(g(z))));

(m^2+n^2-p^2-q^2)/(m^2-2*m*q+n^2+2*n*p+p^2+q^2)

(4)

b := simplify(evalc(Im(g(z))));

(-2*m*p-2*n*q)/(m^2-2*m*q+n^2+2*n*p+p^2+q^2)

(5)

``
"  r:=15;   Lista:=Vector();  Listb:=Vector();  j:=1;  for m from -r to r do   for n from -r to r do   for p  from -r to r do   for q from -r to r do  if p<>0 and q<>0 and m^2-2 m q+n^2+2 n p+p^2+q^2<>0 and bz>=0 then  Lista(j):=a; Listb(j):=b;  j:=j+1;  end if;  end do:  end do;  end do;  end do:  j; :"

435713

(6)

``

``

``

``

``

pointplot(Lista, Listb, symbol = plottools:-point, symbolsize = 1, size = [1200, 1200])

 

NULL

``

``


 

Download 096-Chayley_transform_for_MP_question.mw

i tried to solve a nonlinear ode with numerical method but maple can't solve it and this error occur:

Error, (in dsolve/numeric/bvp) initial Newton iteration is not converging

my maple codes are attached below:

numeriacal_sol.mw

can any help me?

As I understand it, Maple will detect and use the available cores in a system, if the calculation is suitable for multi-core use.

As I am installing Maple on a multi-user cluster, using a scheduler to run maple scripts, I want to ensure the maple jobs only use the number of cores allocated to the job.  

Is it possible to set the number of cores used ? 

If I have misunderstood how Maple works (I am new to it), or if there is a section in the documentation which explains this, please point me in the right direction.  I haven't found this info so far.

if I use the command below, I will get the right answer.

`assuming`([simplify((2*I)*deltas*Pi*BW*Ei(1, (2*I)*deltas*Pi*BW)-(2*I)*deltas*Pi*BW*Ei(1, -(2*I)*deltas*Pi*BW))], [real, BW > 0, deltas > 0, ts > 0, deltas <= ts, 2*BW*ts >= 1])

and the answer is :

2*BW*Pi*deltas*(Pi-2*Si(2*BW*Pi*deltas))

while I use the same command for the compliacated Ei function, it doesn't work any more, I don't know the reason.

 `assuming`([simplify(Ei(1, I*BW*Pi*deltas-(5*I)*BW*Pi*ts)-Ei(1, -I*BW*Pi*deltas+(5*I)*BW*Pi*ts), Ei)], [real, BW > 0, deltas > 0, ts > 0, deltas <= ts, 2*BW*ts >= 1])

 and it's reault is :

Ei(1,I*BW*Pi*deltas-5*I*BW*Pi*ts)-Ei(1,-I*BW*Pi*deltas+5*I*BW*Pi*ts)

 

Dear Maple users

I was using the animation command to create an animation consisting of a running sine curve with a ball following the curve, as shown below. I wonder if it is possible to create it with the Explore command too? Or is this command more limited? I need to be able to add more animations, like I can with the display command ...

Regards, Erik.

Hello people in mapleprimes,

I want to modify the inequality of e assuming d.
But, e/(1-a) assuming(d) does not work well.
How can I do for this?

d:=(a>0,a<1,b>0,c>0);
                   0 < a, a < 1, 0 < b, 0 < c
e:=(1-a)*b>c;
                         c < (1 - a) b
e/(1-a)assuming(d);
                         `*`(c<(1-a)b,1/(1-a))

Thanks in advance.

inequality.mw

 

 

I have made an animate-plot with trace function enabled and I was wondering if there was something I could do so it would show the value of the variable when trace reached frame x and y and plot it in the plot or maybe somewhere else.

Thanks in advance.

Hello,

when trying to solve my equation I always get this error-message:

Error, (in DEtools/convertsys) invalid input: degree expects 1 or 2 arguments, but received 3

To solve, I use this one:

Anfangsbed := phi(0) = 0,  D(phi)(0) = 10;

Einst :=             range = 0..10,            #Integrationsgrenzen
                     relerr = 1e-3,            #Toleranz zur Bestimmung der Integrationsschrittweite
                     optimize = true,          #Optimierung der DGL
                     maxfun = 0,               #Anzahl der Integrationsschritte -- 0 = keine Obergrenze
                     output = operator,        #Lösungsfunktionausgabe: damit Definition als Funktion möglich 
                     stiff = true:

 

Lösung_Bewgl := dsolve({Bewegungsgl_mit_Para, Anfangsbed}, numeric, Einst);

The equation to solve is this one: http://imgur.com/a/SFxI8

What am I doing wrong?

Edit:

This is the worksheet, problem is the last point.

WS.mw

Hi,

Please help me to solve numerically and visualise the following time-delayed integro-differential equation (see also the file attached)

diff(x(t), t$2) = .8/(x(t)^3*exp((2*0.1e-1)*(int(1/x(t), t = 0 .. t-tau))))-1/x(t)^2 , 

with the initial conditions x(0) = 1, (D(x))(0) = 0, and tau=0.1 (variable parameter).

The expected solution is an oscillatory function.

Thanks for your time and advice!

IDE.mw

 

hi--how i can solve following equation?

thanks

 

Eq.mw


Maple Worksheet - Error

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

Download Eq.mw

 

First 26 27 28 29 30 31 32 Last Page 28 of 60