Laplace transform using numerical integration in Python has very poor precision

Question

I have written a function to compute the Laplace transform of a function using scipy.integrate.quad. It is not a very sophisticated function and currently perfo