Thanks for that!
To calculate the integral, I wrote a program in Matlab.
I assign the constans A, B, C, D, E values between 0.1 and 1.0, and use the following command:
y = quadgk(@(tau)FsfaTauZ(tau,z).*exp(tau.*t),-Inf,Inf)
Where FsfaTauZ contains the the function F as you wrote it, with...