Maple 2016 Questions and Posts

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

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

 

Hello,

i have a Problem with maple. I have written this simple procedure

myPi_1:=proc(r)
local i,l;
l:=0;
for i from 1 to floor(r) do
    if isprime(i) and (i mod 4 = 1) then
        l:=l+1;
    end if;
end do;
return l;
end proc:    

 

When i use it like this 

k := myPi_1(80);

it works

But when i triy to plot it

plot(myPi_1(t), t = 0 .. 500, numpoints = 500, thickness = 2, color = black);

i get this message
Warning, expecting only range variable t in expression myPi_1(t) to be plotted but found name myPi_1

And no plot

Thank you for any Help

 

 

 

 

I want to calculate the following integral numerically with required precision.

First, the functions are defined:

G1:=-0.9445379894;
f:= (x) -> 0.9/abs(x-0.4)^(1/3)+0.1/abs(x-0.6)^(1/2);
U1 := unapply(-exp(-x)*(evalf(Int(f(t)*exp(t), t = 0 .. x))+G1)/2-exp(x)*(evalf(Int(f(t)*exp(-t), t = 0 .. x))+G1)/2, x);
U:= unapply(-exp(x)/2*(evalf(Int(f(t)*exp(-t),t=0..x))+G1)+exp(-x)/2*(evalf(Int(f(t)*exp(t),t=0..x))+G1), x);

Next, I calculate the integral in numerical form:

evalf(Int(U1(x)^2+U(x)^2-2*f(x)*U(x), x=0..1, digits=4, method = _Gquad));

If I specify digits=4, Maple return the answer -0.4291

If I use digits=5 or larger, Maple return someting like this

Is it possible to increase precision of calculation?


 

 

Hey there!

 

this is the formula:

diff(ca(t), t) = -3.600000000*10^20*exp(-15098.13790/(340-20*ca(t)))*ca(t), ca(0) = 2

I can easily plot it by solving DE numerically, however I can't seem to find a way to determine what 't' is at ca(t)=0.2

 

thank you so much in advance!

Kind regards, 

Gerard 

Hello People in mapleprimes,

I have a question about collecting part of expression.
Please teach me an answer to this.

I want to modify a1 to b1.


e3.mw

Thanks in advance.

First 27 28 29 30 31 32 33 Last Page 29 of 61