Generally speaking, in computer algebra, "given a function $f\left(x\right)$, the input $\frac{\ⅆ}{\ⅆx}f\left(x\right)$ should return the derivative of $f\left(x\right)$". 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 $\mathrm{laplace}\left(-t\cdot f\left(t\right)\,t\,s\right)$ you would receive $\frac{d}{\mathrm{ds}}\mathrm{laplace}\left(f\left(t\right)\,t\,s\right)$. 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 $\mathrm{setup}\left(\mathrm{computederivatives}=\mathrm{false}\right)$. To query about the status of this new setting enter

and so differentiating returns the derivative computed

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

Reset the value of computederivatives

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.

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

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

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.

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

and the inverse Laplace transform of LT in inert form is

Since by construction $\mathrm{ILT} \= F$, at $t\=1$ we have

This inert form on the left-hand side is numerically evaluated next, directly, without symbolically resolving it to be equal to $\sin \left(\sqrt{2}\right)$. 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

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.

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

$\mathrm{hankel}\left(f\left(t\right)\,t\,s\,\mathrm{\nu }\right)={\int }_0^{\mathrm{\infty }}f\left(t\right)\sqrt{st}{J}_{\mathrm{\nu }}\left(st\right)\ⅆt$

where $J_{\mathrm{\nu }}\left(st\right)$ is the $\mathrm{BesselJ}\left(\mathrm{nu}\, s\cdot t\right)$ 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 $\sqrt{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,

 $\mathrm{hankel}\left(f\left(t\right)\,t\,s\,\mathrm{\nu }\right)={\int }_0^{\mathrm{\infty }}f\left(t\right)t{J}_{\mathrm{\nu }}\left(st\right)\ⅆt$

has the advantage that, with it, the Hankel transform of $\frac{{\partial }^2}{\partial r^2}u\left(r\,t\right)+\frac{\frac{\partial }{\partial r}u\left(r\,t\right)}{r}+\frac{{\partial }^2}{\partial t^2}u\left(r\,t\right)$ is given by the simple ODE form $\frac{d^2}{{\mathrm{dt}}^2} \mathrm{\ℋ}\left(k\,t\right)-k^2\mathrm{\ℋ}\left(k\,t\right)$. 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 $\sqrt{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

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

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

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

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

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

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

Reset to the new default value of the definition.

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

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 $\mathrm{computederivatives}\=\mathrm{false}$, a more convenient setting to use transforms to simplify differential equations, and insert some values in formula_minus

Its Hankel transform is

An example with formula_plus:

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