Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

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

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

As you can see this app performs the trace of a given path r (t), then locate the position vector in a specific time. It also graphs the velocity vector, acceleration, Tangential and Normal unit vectors, along with the Binormal. Very good app developed entirely in Maple for our engineering students.

Plot_of_Position_Vector_UPDATED.mw

https://youtu.be/OzAwShHHXq8

Lenin Araujo Castillo

Ambassador of Maple

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

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 

It seems a large number of people, when initially using maple, wrongly deduce that for example sin(60) is the sin of 60 degrees and not the sin of 1/3 Pi radians.  I believe mathematica's default is degrees.  When a student compares an expression to another but forgets to realize a value is read as radians and not degrees they are perplexed when Maple returns false and Mathematica returns true.

As a suggestion, under tools->options allow a user to be able to change how maple reads values within trigonometric funtions as either radians or degrees.

Most times when someone computes the sin(60) what they really mean in Maple..
is sin(convert(60 degrees, radians))

 

Hello,

I am not sure how exactly to word my question, so forgive the mess ahead.

Say I have two points, A and B, and I want to make them variables in Maple so whener I refer to A, it refers to the point A.

Currenty what I do, is I write the line "A:=(x,y,z)" enter and then the same for B, but what I would like to do, is maybe write both the variables in one statement - something like "A,B:=(x,y,z),(x,y,z)", so when I am solving problems with 5-10 different points, I won't be filling up 1-2 pages with statements on variables and their values.

I am assuming this is a possiblity, but I haven't been able to find anything on-the-line that could show me how to...

Example picture - where A and B should have had been in the same statement, so when I find the vector between them it'll still give the correct answer.

Any help is appreciated here :)

With gratitude:
Toby

What method is implemented in the procedure int with option numeric? Does the method depend on the integrand? Where can I find this information? 

First 933 934 935 936 937 938 939 Last Page 935 of 2224