quadraturesize = posint

This command computes a numerical discrete inverse Laplace transform using a method of Peter den Iseger involving only simple linear algebra and InverseFFT and $\mathrm{O}\left(M\right)$ evaluations of the function f.

For a function given by f, this command outputs a vector of M equally spaced numerical values of the inverse Laplace transform of $f\left(t\right)$, each Delta apart, starting at $t=0$. The algorithm is expected to be accurate for $f$ with a smooth inverse transform.

Intermediate calculations are done using hardware double precision floating-point numbers, so raising the value of $\mathrm{Digits}$ will not increase the precision of the results.

$\mathrm{with}\left(\mathrm{inttrans}\right)\:$

$\mathrm{invlaplace}\left(\frac{s^2}{{\left(s^2+2\right)}^{\frac{3}{2}}}\,16\,0.01\,'\mathrm{numeric}'\right)$

$F≔\mathrm{unapply}\left(\mathrm{invlaplace}\left(\frac{s^2}{{\left(s^2+2\right)}^{\frac{3}{2}}}\,s\,t\right)\,t\right)$

$F≔t\mapsto -\sqrt{2}\cdot t\cdot \mathrm{BesselJ}\left(1\,\sqrt{2}\cdot t\right)+\mathrm{BesselJ}\left(0\,\sqrt{2}\cdot t\right)$

$\mathrm{Vector}\left(16\,i\mapsto \mathrm{evalf}\left(F\left(\left(i-1\right)\cdot 0.01\right)\right)\right)$

$\mathrm{invlaplace}\left(s\mapsto \frac{\mathrm{\Gamma }\left(\frac{4}{3}\right)}{s^{\raisebox{1ex}{$4$}\!\left/ \!\raisebox{-1ex}{$3$}\right.}}\,16\,0.01\,'\mathrm{numeric}'\right)$

den Iseger, Peter. "Numerical Transform Inversion Using Gaussian Quadrature. Probability" in The Engineering and Informational Sciences, Vol. 20, pp. 1-44, 2006.

The inttrans/invlaplace/numeric command was introduced in Maple 2025.

For more information on Maple 2025 changes, see Updates in Maple 2025.