@Arny 

After these corrections youwill have the problem formulated in a way the computer understands it. Regarding your use of DotProduct and CrossProduct it made me think that it may be convenient for you to give a further look at the help page ?Physics,Vectors and pages linked therein to see what are the names of the commands representing the operations.

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

@Arny 

To get the right display, enter Setup(mathematicalnotation = true). But mostly: you need a more modern Maple.

Time ago, these integrals were computed as undefined. Nowadays they are computed as a distribution (Dirac). I don't remember in which version I introduced Diracs but probably Maple 2015, or Maple 18. Note as well that there are updates for the Physics library distributed in the Maplesoft R&D Physics webpage. Physics is almost duplicating its size at every release till 2016 ... and it is still growing in many directions.

These updates are posted basically every week, sometimes several times per week. These updates started being distributed with Maple 18 and there is a link for updating Maple 18 in that R&D page - I don't remember if it contains this change that makes int return Diracs. Anyway there is no update for Maple 17, that is by now somewhat old.

PS: a note on my previous reply, you can save some typing using [x, y, z] =~ -infinity..infinity instead of repeating the integration range three times as I did in the reply worksheet.

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

@Arny 

Download VectorIntegral.mw

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

@Louis Lamarche,  tomleslie

I misread here! Sorry. The problem you mention was indeed an issue. In Physics, the order for spherical coordinates is [r, theta, phi] while in VectorCalculus is [r, phi, theta], this is taken into account in one of the conversions (to Physics) but was not taken into account in the other conversion (to VectorCalculus).

This issue in the conversion to VectorCalculus of a vector of the Physics package in spherical coordinates is fixed now. The fix is available to everybody within the updated Physics library distributed in the Maplesoft R&D Physics webpage.

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

Very nice everything!

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

@_Maxim_ , Axel Vogt,

Big oops. You are both right, this conversion of MeijerG to hypergeom has a problem for abs(z) > 1 in general, regardless of Re(z) < 0, case p=q and m+n<=p. Have you noticed any other case, so that when fixing this conversion I can check them as well. Thanks

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

@_Maxim_ 

The plots, as well as the original example of this post, are correct for Re(z) < 0 or abs(z) < 1, not just the latter. The problem is in rather ancient code from 1996 ... I will give this a further look.

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

@Markiyan Hirnyk 

The improvements mentioned by Eithne are in the computation of symbolic (so exact, not numerical) solutions to "PDEs & Boundary Conditions" problems. I just posted a set of examples illustrating these improvements.

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

@MrMarc 

The article behind this map of topics, linked in the post, is what gives context to the map. It says right after in its title:  "Explore the deepest mysteries at the frontier of fundamental physics, and the most promising ideas put forth to solve them". The theories mapped, thus, refer to fundamental physics only, as the map itself shows. This map is indeed an excellent summary of the ideas being explored in the area, that not surprisingly revolve around general relativity in one way or another. Mainly on how to unify this theory of spacetime with quantum mechanics, the topic of Nima's presentation, Quantum Mechanics and Spacetime in the 21st Century, also linked in the post.

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

@Pascal4QM 

:) Thanks. This problem you indirectly focused, however, after some discussion, we agreed at this end that it is indeed unexpected behavior in the kernel. Not OK. In the example mentioned F[1,2](x,y) should return f[1,2](x,y). The issue is now tracked to be fixed in the kernel. Meantime, and only for Physics tensors, the library (not kernel) tweak I introduced yesterday works around the issue.

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

@_Maxim_ ,  Pascal4QM,  vv

Yes, your suggestion works, but there is a consistency problem. Consider

with(Physics);
Define(F[mu, nu] = Matrix(4, symbol = f));

Now, without defining f as a procedure, just enter

F[1,2]
                  f[1,2]

Good. Next:

F[1,2](x,y)
                 F[1,2](x,y)

I was expecting f[1,2](x,y). This is the problem I referred to in my previous reply, perhaps not clearly. I expect F[1,2] to be evaluated, so become f[1,2], then f[1,2] to be applied to (x,y).

Indeed this problem is not visible when you define f itself as a procedure, what you _Maxim_ and PascaQM are saying, eg one that makes use of its indices for the purpose of Pascal's example (see vv's reply, many Maple procs are defined that way, for instance, Physics:-d_),

So why my fuzz? Perhaps because these things, if let in place, for me they attempt against the intuition on how the language works.

Anyway, I fixed this now, as explained in the previous reply, so that we now receive f[1,2](x,y) even when f is not a procedure. The adjustment in the code is available to everybody now, as usual in the Maplesoft R&D Physics webpage.

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

This is indeed a typo in the code, not in the rewriting in terms of abs  (conversion, I dislike this word ...) but in the rewriting of arctan(y, x) in terms of argument. The network works avoiding repetition of formulas when possible, and arctan goes to argument, and from there to abs. Indeed, in addition, the formula is only valid if both x and y are reals. I will fix this and see if this can still be included in the upcomming 2017.3.

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

@wndrslope 

First of all on how it works: the coordinates used by the Physics:-Vectors commands are fixed in that, for cylindrical, they are always [rho, phi, z] and cannot be changed. For spherical they are [r, theta, phi], for cartesian they are [x, y, z] and they cannot be changed too. Their definitions are on the pages ?Physics:-Vectors and ?Vectors,ChangeBasis. The corresponding unit vectors also cannot be changed. So, even if you enter Coordinates(cylindrical) followed by Coordinates(X = [r,theat,z,t]), the first call sets an alias(X = (rho, phi, z, t)) and the second one resets that alias to alias(X = (r, theta, z, t)), and regardless of your calls to Coordinates, for the command Vectors:-Nabla, the cylindrical coordinates are always rho, phi and z.

Now, why would you want to use [r, theta, z] with the geometrical meaning (and the corresponding mathematical relationships) of cylindrical coordinates instead of [rho, phi, z]? It is only a difference in 'name' and 'look'. Well, the answer could still be any, say you find more comfortable to write r instead of rho, or the book you are using uses r, theta instead of rho, phi. So here is what you could do to make Vectors:-Nabla and all the commands of Physics:-Vectors work with a different name and look.
 

Download ChangingDefaultLookOfCoordinates.mw

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

@wndrslope 

Yes it is possible, Setup(coordinates = cylindrical) or Coordinates(cylindrical) will both set X = (rho, phi, z, t), or Coordinates(Y = cylindrical) will do the same with Y. And you can always set Coordinates(X = [a,b,c,d]) where a, b, c, d can be any object of type 'name' (see ?type,name). The advantage of using rho, phi, z is that the Physics:-Vectors use those with the exact geometrical meaning of cylindrical coordinates (this meaning is automatically taken into account, for instance, if you Setup(geometricdifferentiation = true)). The documentation of the package has many pages (the Physics environment/package has many commands and some conventions) but all this is explained in the respective pages.

One thing to the side, sometimes relevant: for the signatures +++- or ---+, the position of the time component is 4. If however you want to work with time in position 1, you need to change the signature to either -+++ or +---; for that you use for instance Setup(signature = `-+++`).

And a note about Setup : this is a command with a myriad of options, all explained in ?Physics,Setup, as a mimicry of what we do with paper and pencil: we set a number of conventions in our mind, then proceed computing. On the worksheet, you set these conventions using Setup. What people sometimes miss is that if you enter Setup(); so without any arguments and press enter, an applet pops up displaying all  the settings valid at that moment, and you can change them directly from the applet. This makes setting things simpler: you can rapidly see all the current settings, organized, and change them sometimes with a few clicks.

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

@_Maxim_ 

The FunctionAdvisor command is for retrieving mathematical information in a computationally active way (not just display it), and it does that in two different ways

1. it return mathematical information in an organized way, grabbing it from all around the Maple library
2. it extends that functionality to algebraic expressions.

For an example of 1. enter, for instance, FunctionAdvisor(GAMMA). For an example of 2. see, for instance, ?FunctionAdvisor[branch_cuts].

Now, when it is not the case of "2.", e.g. when computing just what a Maple command already computes naturally, for example to compute asympt(hypergeom([1, 1], [2, 2], -z), z) , the idea is to use the command. Of course the FunctionAdvisor - that for asymptotic expansions only uses asympt (ie there is no "2." functionality in this case, asympt is already the command for that purpose) could also just forward the stuff with "-z" argument to asympt, of course with no gain whatsoever. Either way the point remains that what the FunctionAdvisor does regarding asymptotics is just what asympt does, no more, because it just forwards your input to asympt.