Hi

On the latest Maple 2022 version, using dsolve with 1st order ODE and Dirac function returns an incorrect solution

dsolve({D(y)(t) + y(t)/tau = Dirac(t)/tau, y(0) = 0})

See attached image. The correct impulse response should be :

y(t) = exp(-t/tau)*Heaviside(t)/tau

The returned solution is 0 for t = 0 and scaled by a 1/2 for t > 0. I never had such an issue with older Maple versions. Is this a bug or am I doing something wrong

Here is the expected solution followed by the Maple solution: