MaplePrimes Questions

equations are as folllows

1) x>0

2) x>2(y-a-1)

3)( x+a+1)>y>(x+a-1)

how to plot in Maple the region of solution for these equations in x and y plane keeping "a" as it is. 

Helo,

how can i export big equations to MS WORD? I trying copy by clicking RMB on equation and paste to WORD but i see only "null". What i can doo now?

I have a complicated set of first-order differential equations, which Maple seems capable of solving. So far so good. The obtained solution is fed into another complicated set of first-order differential equations, which, again, Maple seems capable of solving.

But this final (combined) solution cannot possibly be the complete one as it does not contain a specific (sub)solution which DOES satisfy both sets of equations (checked in Maple by evaluating these equations for that specific solution). And there is no warning raised concerning any solutions that may have been lost.

The problem seems to be that the integrational constants associated with solving (in turn) the two sets of equations become messed up together, or erroneously 'reused', in some mysterious way. For if I rename the integration constants of the first solution before feeding it into the second set of equations, then this specific subsolution IS contained in the final solution.

Therefore the following questions:

1.) Could there possibly be some erroneous 'reuse' of integration constants going on? It seems to me that in solving the second set of equations, PDEtools:-Solve would have to make sure that it numbers any new integration constants in such a way that there would be no conflict with the ones obtained in solving the first set of equations.

2.) Is there some (global) variable that determines what letter is assigned for the integration constants? I could use any such to easily switch names between the first and second solving.

PS: I think, it makes little sense to upload any code/worksheet here.

what is the code to output every symbol assigned to a unique individual operator that is accepted by the maple interface, example what code can i assign to the function:

List_of_all_Operators_Available();

which will produce output in a list of single character strings

["*", "-","+","=",....]

Hello,

I am trying to plot this map

y:=y[t+1] +e(x[t]) mod N
x:=x[t]+y[t+1] mod N

where e(x[t]) is a piecewise function defined by

e(x[t])=1 if 0<=x< floor(N/2) OR

e(x[t])=-1 if floor(N/2) <=x<N
for any initial conditions [x,y] and N (but must be integers).
What I did is two separate functions: 1) for e(x[t]) functiion is given by

E := proc (x, N) global e; if 0 <= x and x <= floor((1/2)*N) then e := 1 else e := -1 end if end proc

and 2) for the main function

S := proc (z) local i, x, y, c, e; global N, E; y := `mod`(z[2]-E, N); x := `mod`(z[1]+z[2]+E, N); return [x, y] end proc

My orbit procedure (so that the output is list of integers [x,y]) is given by:

This DataPlot does not give the list output of integers. I do not know whether this function is correctly written. 
Can somebody help me producing the list of the orbit [x,y]?
Thank you very much for helping!

Hello every body i need help with my program. how can i plot r&s versus each other?



 

restart

with(plots):

with(PDEtools, casesplit, declare)

[casesplit, declare]

(1)

with(DEtools, gensys)

[gensys]

(2)

declare((Omega, H)(z));

Omega(z)*`will now be displayed as`*Omega

 

H(z)*`will now be displayed as`*H

(3)

`&epsilon;` := -1

-1

(4)

h := 70

70

(5)

r := 0.3e-3

0.3e-3

(6)

n := 42

42

(7)

``

``

ode1 := -(diff(Omega(z), z))+Omega(z)*((3*(1-2*`&epsilon;`*h*sqrt(r)/H(z)-Omega(z)+2*Omega(z)^(3/2)/(3*n)))*(1+h*`&epsilon;`*sqrt(r)/sqrt(h^2*r+H(z)^2*(1-2*`&epsilon;`*h*sqrt(r)/H(z))))-2*sqrt(Omega(z))/n)/(1+z) = 0:

ode2 := -(diff(H(z), z))-(3*H(z)*(1/2))*(1-2*`&epsilon;`*h*sqrt(r)/H(z)-Omega(z)+2*Omega(z)^(3/2)/(3*n))*(1+h*`&epsilon;`*sqrt(r)/sqrt(h^2*r+H(z)^2*(1-2*`&epsilon;`*h*sqrt(r)/H(z))))/(1+z) = 0:

q := -1+(3/2*(1-2*`&epsilon;`*h*sqrt(r)/H(z)-Omega(z)+2*Omega(z)^(3/2)/(3*n)))*(1+h*`&epsilon;`*sqrt(r)/sqrt(h^2*r+H(z)^2*(1-2*`&epsilon;`*h*sqrt(r)/H(z)))):

``

sys := {ode1, ode2}


ics := {H(0) = h, Omega(0) = .77}:
sol := dsolve(`union`(sys, ics), numeric, output = listprocedure, stiff = true):``

 

 

 

 

 

r2 := q+2*q^2+(1+z)*(diff(q, z)):

``

``

s := (r2-1)/(3*(q-1/2)):

 

 

 

 



 

Download rs.mw
 

 

 



 

Download rs.mw

 

Hello

I am trying to plot a harmonic oscillation, but maple keeps rewriting my defintion so when I plot, it gives me a wrong graph...

How do I stop it from doing this, I am assuming it is this that causes the problem, and not the plot command?

Here is what I get:

And here is what I want - made with geogebra.

If it isn't the definition that is the problem, what plot command should I be using then?

With gratitude
Toby

Maple can calculate one Inverse Laplace Transform and give the result.

But I find this integral diverges because this term exp(alpha1*(Dp*alpha1*t-lh+x)) is infinite in one interval.

Hey. I have the following Diff. Equation:

diffeq := diff(y(t),t,t,t,t)+10*diff(y(t),t,t)+169*y(t)=0

and I have the conditions y(0) = 0, y'(0) = 1, and y(t) -> 0 for t -> infinity 

I know how to do with the first two conditions, but how do you do it with a boundary? I've read the posts made earlier on here, but I can't seem to figure it out. I've tried:

dsolve({diffeq,y(0)=0,D(y)(0)=1,D(y)(infinity)=0},y(t))
 

Returns nothing.

Any help would be appreciated.

 

Hi everyone,

 

I am trying to plot the standard map equations which is defined by

p_{n+1}=p_{n}+k sin x_{n}  (mod 2pi)
x_{n+1}=x_{n}+p_{n}+k sin x_{n} (mod 2pi).

 

Is there is any simple code for plotting this map in Maple? 

Thank you in advance.

I would like to publish a technical paper about a renewable energy with you. I use Maple 2016. What i need?. Thanks!

Hello again,

I had to download an add-on package to make this graph in the first place, so it might not be possible to help me.
Additionally, I have no idea how to translate "sumkurve" from my native language (danish) to English, which is also why I had to download the add-on package in the first place....

As you can see in the picture below, I have a table at the top, with a matrix underneath, then with the "plotSumkurve(A)" I got this beautiful graph that also shows the quartiles (kvartiler).

But what if I wanted to find the Y value that corresponds to X=12, any idea on how to do this would be great appreciated.

With gratitude,
Toby

#first_question :how can i solve set of nonlinear ODEs,faster or using any packages ?
#second_question :what can be some boundary conditions for this type of nonlinear ODEs? how many BCs are required for this set of nonlinear ODEs? ( to use numeric solution)

 


 

restart:with(DEtools):with(DifferentialAlgebra):

eq[1]:=diff(N(r),r$2)+2/r*diff(N(r),r)+diff(phi(r),r)/phi(r)*diff(N(r),r)-mu^2/(32*phi(r))*N(r);

diff(diff(N(r), r), r)+2*(diff(N(r), r))/r+(diff(phi(r), r))*(diff(N(r), r))/phi(r)-(1/32)*mu^2*N(r)/phi(r)

(1)

eq[2]:=diff(phi(r),r$2)+2/r*diff(phi(r),r)-1/2*diff(phi(r),r)^2/phi(r)-8*diff(N(r),r)^2/(omega*(1-2*G*M/r))*phi(r);

diff(diff(phi(r), r), r)+2*(diff(phi(r), r))/r-(1/2)*(diff(phi(r), r))^2/phi(r)-8*(diff(N(r), r))^2*phi(r)/(omega*(1-2*G*M/r))

(2)

dsolve({eq[1],eq[2]});

``


 

Download nonlinear_ODE.mw

I would like to have my calculus results in its simplified form; below are some of the scenarios that can explain this.

 

 

As you can see the results are much simplified and in reduced form compared to what Maple gives. This has been tested for all integrals with roots and integrals with trigonometry functions. Is there any workaround in Maple that I would use to get it like the textbook result.

Thanks

 

As I work a lot with lists (need to merge them frequently), I would like to redefine the `union` operator in such a way that it will merge two lists together. It would makes thing more efficient than writing 

list1 := [a,b,c];
list2 := [d,e,f];

newlist = [op(list1), op(list2)];

I've already tried

`&union` := proc(list1, list2)
  [op(list1), op(list2)];
end proc:

[a,b,c] &union [d,e,f]

but it's not what I'm looking for. In fact I could have used any name after the `&`... and really don't like having to type the & at the beginning (is there a way to define an infix operator without having to use the `&`?)

Maybe creating a module to override  the original definition?

my_module := module()
  export `union`:

  `union` := proc(foo,bar)
      if some_trigger then
         # return something
      else
         # use global `union` definition
      end if
   end proc:

end module;

Any suggestion?

P.S.

- I generally don't use the original set union operator, so redefining it is not an issue;

- I have to use lists 

4 5 6 7 8 9 10 Last Page 6 of 1488