Integral Transforms - Maple Help
For the best experience, we recommend viewing online help using Google Chrome or Microsoft Edge.

Online Help

All Products    Maple    MapleSim

Home : Support : Online Help : System : Information : Updates : Maple 2020 : Integral Transforms

Integral Transforms

Integral transforms are special integrals that appear frequently in mathematical physics and that have remarkable properties.  For Maple 2020, the implementation of these integrals has been extended in several ways, making them more useful for a variety of applications, including computing integrals, finding exact solutions to PDEs with boundary conditions (see also what's new in Maple 2020 for PDE & Boundary Conditions).  Improvements include the option to compute derivatives, numeric support, an alternate definition of the Hankel transform, and the ability to compute more transforms.


As background, these integrals are represented by the commands of the inttrans package:




Three of these commands, addtable, savetable, and setup (new as of Maple 2020) are "administrative" commands, while the others are computational representations for definite integrals. For example,

FunctionAdvisorintegral_form, fourier

a,b,z=aⅇIbzⅆb,with no restrictions on a,b,z


FunctionAdvisorintegral_form, mellin

a,b,z=0abz1ⅆb,with no restrictions on a,b,z


For all integral transform commands, the first argument is the integrand, the second one is the dummy integration variable of a definite integral and the third one is the evaluation point (also called transform variable).


Compute derivatives: yes or no.

Numerical evaluation

Two Hankel transform definitions

More integral transform results

Compute derivatives: yes or no.

Generally speaking, in computer algebra, "given a function fx, the input ⅆⅆxfx should return the derivative of fx". For the purpose of using integral transforms to solve differential equations, however, the implementation in previous Maple releases worked in the opposite direction: if you were to input the result of the derivative, you would receive the derivative representation. For example, to the input laplacetft,t,s you would receive ddslaplaceft,t,s. To extend the capabilities of the system and have the best result in any scenario, a new command, setup, has been added to the package, so that you can set whether or not to compute derivatives, and the default has been changed to computederivatives = true while the previous behavior is obtained only if you input setupcomputederivatives=false. To query about the status of this new setting enter





and so differentiating returns the derivative computed

%diff = difflaplaceft,t,s,s



while changing this setting to work as in previous releases you have this computation reversed: you input the output (1.2) and you get the corresponding input

setupcomputederivatives = false





Reset the value of computederivatives

setupcomputederivatives = true





In summary: by default, derivatives of all the 11 integral transforms are now computed. This setting can be changed any time you want within a Maple session without restarting, and changing it does not have any impact on the performance of intsolve, dsolve and pdsolve to solve differential equations using integral transforms.


Numerical evaluation

In previous releases, integral transforms had no numerical evaluation implemented. This is in the process of changing. So, for example, to numerically evaluate the inverse laplace transform (invlaplace command), three different algorithms have been implemented: Gaver-Stehfest, Talbot and Euler, following the presentation by Abate and Whitt, "Unified Framework for Numerically Inverting Laplace Transforms", INFORMS Journal on Computing 18(4), pp. 408–421, 2006.


For example, consider the exact solution to this partial differential equation subject to initial and boundary conditions


iv  ux, 0 = 0, u0, t = 1:


Note that these two conditions are not entirely compatible: the solution returned cannot be valid for x=0 and t=0 simultaneously. However, a solution discarding that point does exist and is given by

sol  pdsolvepde, iv



Verifying the solution, one condition remains to be tested




Since we now have numerical evaluation rules, we can test that what looks different from 0 in the above is actually 0.

zero  1



Add a small number to the initial value of t to skip the point t=0:


The default method used is the method of Euler sums and the numerical evaluation is performed as usual using the evalf command. For example, consider





The Laplace transform of F is given by

LT  laplaceF,t,s



and the inverse Laplace transform of LT in inert form is

ILT  %invlaplaceLT, s, t



Since by construction ILT = F, at t=1 we have

evalILT=F, t=1



This inert form on the left-hand side is numerically evaluated next, directly, without symbolically resolving it to be equal to sin2. The result of both sides is consistent:




In addition to the standard use of evalf to numerically evaluate inverse Laplace transforms, one can invoke each of the three different methods implemented using the MathematicalFunctions:-Evalf command




Evalflhs, method = Talbot



MathematicalFunctions:-Evalflhs, method = GaverStehfest



MathematicalFunctions:-Evalflhs, method = Euler



Regarding the method we use by default: from a numerical experiment with varied problems, we have concluded that our implementation of the Euler (sums) method is faster and more accurate than the other two.


Two Hankel transform definitions

In previous Maple releases, the definition of the Hankel transform was given by


where Jνst is the BesselJnu, st function. This definition, sometimes called alternative definition of the Hankel transform, has the advantage of having a large table of transforms presented in the literature (ref.[1]) but has the inconvenience of the square root st in the integrand, complicating the form of the Hankel transform for the Laplacian in cylindrical coordinates. On the other hand, the other definition in the literature,


has the advantage that, with it, the Hankel transform of 2r2ur,t+rur,tr+2t2ur,t is given by the simple ODE form d2dt2 ℋk,tk2ℋk,t. Several other transforms also acquire a simpler form. So for Maple 2020 we added an algorithm that takes formulas for one definition and transforms them into formulas for the other definition, and have aligned Maple with this simpler definition (no st in the integrand), while keeping the previous definition as an alternative.


Hence, by default, when you load the inttrans package, the new definition in use for the Hankel transform is

hankelft,t,s,nu: % = convert%,Int



You can change this default so that Maple works with the alternative definition as in previous releases.  For that purpose, use the new inttrans:-setup command (which you can also use to query about the definition in use at any moment):




The change in the default definition implemented in Maple 2020 is automatically taken into account by other parts of the Maple library using the Hankel transform. For example, the differentiation rule with the new definition is

%diff = diffhankelft,t,z,nu,z



This differentiation rule resembles (is connected to) the differentiation rule for BesselJ, and this is another advantage of the new definition.

%diff = diffBesselJnu,z,z



Furthermore, several transforms have acquired a simpler form, as for example:

%hankel &equals; hankelexpI a rr&comma; r&comma; k&comma; 0 assuming a &gt; 0&comma; k < a


Let's compare: make the definition be as in previous releases.

setupalternativehankeldefinition &equals; true



hankelft&comma;t&comma;s&comma;nu&colon; % &equals; convert%&comma;Int



The differentiation rule with the previous (alternative) definition was not as simple:

%diff &equals; diffhankelft&comma;t&comma;s&comma;nu&comma;s



And the transform (3.5) was also not so simple:

%hankel &equals; hankelexpI a rr&comma; r&comma; k&comma; 0 assuming a &gt; 0&comma; k < a



Reset to the new default value of the definition.

setupalternativehankeldefinition &equals; false



hankelft&comma;t&comma;s&comma;nu&colon; % &equals; convert%&comma;Int




[1] Oberhettinger Fritz, Tables of Bessel Transforms. Springer Verlag, 1972.


More integral transform results


Maple is now able to compute more transforms. For example, for the Hankel transform, consider the operators



Being able to transform these operators into algebraic expressions or differential equations of lower order is key for solving ODE and PDE problems with Boundary Conditions. To illustrate that, set computederivatives&equals;false, a more convenient setting to use transforms to simplify differential equations, and insert some values in formula_minus

setupcomputederivatives &equals; false&semi;



simplifyevalformula_minus&comma; nu &equals; 6&comma; m &equals; 3



Its Hankel transform is




An example with formula_plus:

simplifyevalformula_plus&comma; nu &equals; 7&comma; m &equals; 4






Not just differential operators but also several new transforms are now computable











fouriersinlna&plus;xxb&comma;x&comma;y assuming 0<a&comma;0<b



fouriersinlna&plus;xxax2&comma;x&comma;y assuming 0<a



fouriersinlna2&plus;x2b2&plus;x2x&comma;x&comma;y assuming 0<a&comma;0<b



fouriercoslna&plus;xxax&comma;x&comma;y assuming 0<a