Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

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? 

plot3d of procedure Sievert correctly displays the constant curvature Sievert surface, but the procedure uses the deprecated command evalm.

What Maple 2016 statement(s) would create the same value of X in Sievert?

Sievert := proc (B)

local a, b, denom, m, X;

a := sinh(B)*u; b := cosh(B)*v;

denom := sinh(B)*((cosh(2*a)-cos(2*b))*cosh(2*B)+2+cosh(2*a)+cos(2*b));

m := cosh(B)*[sinh(a), sin(b)*cos(v), sin(b)*sin(v)]+[0, -cos(b)*sin(v), cos(b)*cos(v)];

X := evalm([u, 0, 0]-8*cosh(B)*cosh(a)*m/denom);

end proc:

plot3d(Sievert(.75), u = -2.5 .. 2.5, v = -10.5 .. 10.5, scaling = constrained, grid = [30, 100], style = patch, shading = xy, lightmodel = light3, orientation = [-3, 140], title = "Sievert's surface", titlefont = [Courier, bold, 14]);

Can any one help me with an example of how to use module in embedded components where a (action:-
 

 

action := module () export arithmatic, times;  arithmatic := proc (a, b) options operator, arrow; a+b end proc; arithmatic := proc (a, b) options operator, arrow; a-b end proc; times := proc (a, b) options operator, arrow; `mod`(a*b, 5) end proc end module

_m1580648227872

(1)

action:-arithmatic(2, 5)

-3

(2)

action:-times(2, 6)

2

(3)

``

I want a procedure (say plot math container in a plot area) in place of arithmatic above.

OR I want a number of calculations under arithmatic ( a+b and a-b)

What is the method I should adopt? How do i insert another module or procedure inside the above module?

Which help page will help me?

Thanks.

Ramki

NULL


 

Download ModuleExample.mwModuleExample.mw

get command is used)?

I attach a module i made for edit and comment.

Thanks.

Ramakrishnan V

Why doesn't Maple provide enough white space when user reaches end of document ? As you can see below, I have reached almost to the end of the page. I believe the default process would be to bring the cursor to eye-level (middle of the screen). Now in this case, I have to hit RETURN enough times to bring the cursor to my eye level.

 

Anyone else find this as an issue ?

 

 


I'm back from presenting work in the "23rd Conference on Applications of Computer Algebra -2017" . It was a very interesting event. This fifth presentation, about "The Appell doubly hypergeometric functions", describes a very recent project I've been working at Maple, i.e. the very first complete computational implementation of the Appell doubly hypergeometric functions. This work appeared in Maple 2017. These functions have a tremendous potential in that, at the same time, they have a myriad of properties, and include as particular cases most of the existing mathematical language, and so they have obvious applications in integration, differential equations, and applied mathematics all around. I think these will be the functions of this XXI century, analogously to what happened with hypergeometric functions in the previous century.

At the end, there is a link to the presentation worksheet, with which one could open the sections and reproduce the presentation examples.
 

The four double-hypergeometric Appell functions,

a complete implementation in a computer algebra system

 

Edgardo S. Cheb-Terrab

Physics, Differential Equations and Mathematical Functions, Maplesoft

 

Abstract:
The four multi-parameter Appell functions, AppellF1 , AppellF2 , AppellF3  and AppellF4  are doubly hypergeometric functions that include as particular cases the 2F1 hypergeometric  and some cases of the MeijerG  function, and with them most of the known functions of mathematical physics. Appell functions have been popping up with increasing frequency in applications in quantum mechanics, molecular physics, and general relativity. In this talk, a full implementation of these functions in the Maple computer algebra system, including, for the first time, their numerical evaluation over the whole complex plane, is presented, with details about the symbolic and numerical strategies used.

Appell Functions (symbolic)

 

 

The main references:

• 

P. Appel, J.Kamke de Feriet, "Fonctions hypergeometriques et Hyperspheriques", 1926

• 

H. Srivastava, P.W. Karlsson, "Multiple Gaussian Hypergeometric Series", 1985

• 

24 papers in the literature, ranging from 1882 to 2015

 

Definition and Symmetries

   

Polynomial and Singular Cases

   

Single Power Series with Hypergeometric Coefficients

   

Analytic Extension from the Appell Series to the Appell Functions

   

Euler-Type and Contiguity Identities

   

Appell Differential Equations

   

Putting all together

   

Problem: some formulas in the literature are wrong or miss the conditions indicating when are they valid (exchange with the Mathematics director of the DLMF - NIST)

   

Appell Functions (numeric)

 

 

Goals

 

• 

Compute these Appell functions over the whole complex plane

• 

Considering that this is a research problem, implement different methods and flexible optional arguments to allow for:

a) comparison between methods (both performance and correctness),

b) investigation of a single method in different circumstances.

• 

Develop a computational structure that can be reused with other special functions (abstract code and provide the main options), and that could also be translated to C (so: only one numerical implementation, not 100 special function numerical implementations)

Limitation: the Maple original evalf command does not accept optional arguments

 

The cost of numerically evaluating an Appell function

 

• 

If it is a special hypergeometric case, then between 1 to 2 hypergeometric functions

• 

Next simplest case (series/recurrence below) 3 to 4 hypergeometric functions plus adding somewhat large formulas that involve only arithmetic operations up to 20,000 times (frequently less than 100 times)

• 

Next simplest case: the formulas themselves are power series with hypergeometric function coefficients; these cases frequently converge rapidly but may involve the numerical evaluation of up to hundreds of hypergeometric functions to get the value of a single Appell function.

 

Strategy for the numerical evaluation of Appell functions (or other functions ...)

 

 

The numerical evaluation flows orderly according to:

1) check whether it is a singular case

2) check whether it is a special value

3) compute the value using a series derived from a recurrence related to the underlying ODE

4) perform an sum using an infinite sum formula, checking for convergence

5) perform the numerical integration of the ODE underlying the given Appell function

6) perform a sequence of concatenated Taylor series expansions

Examples

   

Series/recurrence

   

Numerical integration of an underlying differential equation (ODEs and dsolve/numeric)

   

Concatenated Taylor series expansions covering the whole complex plane

   

Subproducts

 

Improvements in the numerical evaluation of hypergeometric functions

   

Evalf: an organized structure to implement the numerical evaluation of special functions in general

   

To be done

   


 

Download Appell_Functions.mw   
Download Appell_Functions.pdf

Edgardo S. Cheb-Terrab
Physics, Differential Equations and Mathematical Functions, Maplesoft

I'm back from presenting work in the "23rd Conference on Applications of Computer Algebra -2017" . It was a very interesting event. This fourth presentation, about "The FunctionAdvisor: extending information on mathematical functions with computer algebra algorithms", describes the FunctionAdvisor project at Maple, a project I started working during 1998, where the key idea I am trying to explore is that we do not need to collect a gazillion of formulas but just core blocks of mathematical information surrounded by clouds of algorithms able to derive extended information from them. In this sense this is also unique piece of software: it can derive properties for rather general algebraic expressions, not just well known tabulated functions. The examples illustrate the idea.

At the end, there is a link to the presentation worksheet, with which one could open the sections and reproduce the presentation examples.
 

The FunctionAdvisor: extending information on mathematical functions

with computer algebra algorithms

 

Edgardo S. Cheb-Terrab

Physics, Differential Equations and Mathematical Functions, Maplesoft

 

Abstract:

A shift in paradigm is happening, from: encoding information into a database, to: encoding essential blocks of information together with algorithms within a computer algebra system. Then, the information is not only searchable but can also be recreated in many different ways and actually used to compute. This talk focuses on this shift in paradigm over a real case example: the digitizing of information regarding mathematical functions as the FunctionAdvisor project of the Maple computer algebra system.

The FunctionAdvisor (basic)

   

Beyond the concept of a database

 
  

" Mathematical functions, are defined by algebraic expressions. So consider algebraic expressions in general ..."

Formal power series for algebraic expressions

   

Differential polynomial forms for algebraic expressions

   

Branch cuts for algebraic expressions

   

The nth derivative problem for algebraic expressions

   

Conversion network for mathematical and algebraic expressions

   

References

   


 

Download FunctionAdvisor.mw

Download FunctionAdvisor.pdf

Edgardo S. Cheb-Terrab
Physics, Differential Equations and Mathematical Functions, Maplesoft

I'm back from presenting work in the "23rd Conference on Applications of Computer Algebra -2017" . It was a very interesting event. This third presentation, about "Computer Algebra in Theoretical Physics", describes the Physics project at Maplesoft, also my first research project at University, that evolved into the now well-known Maple Physics package. This is a unique piece of software and perhaps the project I most enjoy working.

At the end, there is a link to the presentation worksheet, with which one could open the sections and reproduce the presentation examples.
 

 

 

Computer Algebra in Theoretical Physics

 

Edgardo S. Cheb-Terrab

Physics, Differential Equations and Mathematical Functions, Maplesoft

 

Abstract:

 

Generally speaking, physicists still experience that computing with paper and pencil is in most cases simpler than computing on a Computer Algebra worksheet. On the other hand, recent developments in the Maple system have implemented most of the mathematical objects and mathematics used in theoretical physics computations, and have dramatically approximated the notation used in the computer to the one used with paper and pencil, diminishing the learning gap and computer-syntax distraction to a strict minimum.

 

In this talk, the Physics project at Maplesoft is presented and the resulting Physics package is illustrated by tackling problems in classical and quantum mechanics, using tensor and Dirac's Bra-Ket notation, general relativity, including the equivalence problem, and classical field theory, deriving field equations using variational principles.

 

 

 

 

... and why computer algebra?

 

We can concentrate more on the ideas instead of on the algebraic manipulations

 

We can extend results with ease

 

We can explore the mathematics surrounding a problem

 

We can share results in a reproducible way

 

Representation issues that were preventing the use of computer algebra in Physics

   

Classical Mechanics

 

*Inertia tensor for a triatomic molecule

   

Quantum mechanics

 

*The quantum operator components of  `#mover(mi("L",mathcolor = "olive"),mo("→",fontstyle = "italic"))` satisfy "[L[j],L[k]][-]=i `ε`[j,k,m] L[m]"

   

*Unitary Operators in Quantum Mechanics

 

*Eigenvalues of an unitary operator and exponential of Hermitian operators

   

*Properties of unitary operators

 

 

Consider two set of kets " | a[n] >" and "| b[n] >", each of them constituting a complete orthonormal basis of the same space.

*Verify that "U=(&sum;) | b[k] >< a[k] |" , maps one basis to the other, i.e.: "| b[n] >=U | a[n] >"

   

*Show that "U=(&sum;) | b[k] > < a[k] | "is unitary

   

*Show that the matrix elements of U in the "| a[n] >" and  "| b[n] >" basis are equal

   

Show that A and `&Ascr;` = U*A*`#msup(mi("U"),mo("&dagger;"))`have the same spectrum (eigenvalues)

   

Schrödinger equation and unitary transform

   

Translation operators using Dirac notation

   

*Quantization of the energy of a particle in a magnetic field

   

Classical Field Theory

 

The field equations for the lambda*Phi^4 model

   

*Maxwell equations departing from the 4-dimensional Action for Electrodynamics

   

*The Gross-Pitaevskii field equations for a quantum system of identical particles

   

General Relativity

 

Exact Solutions to Einstein's Equations  Lambda*g[mu, nu]+G[mu, nu] = 8*Pi*T[mu, nu]

   

*"Physical Review D" 87, 044053 (2013)

   

The Equivalence problem between two metrics

   

*On the 3+1 split of the 4D Einstein equations

   

Tetrads and Weyl scalars in canonical form

   

 

 


 

Download Physics.mw

Download Physics.pdf

Edgardo S. Cheb-Terrab
Physics, Differential Equations and Mathematical Functions, Maplesoft

First 924 925 926 927 928 929 930 Last Page 926 of 2215