I’m trying to implement a convolution integral in Python with limits, but I cannot get it to work.

I have this convolution integral:

Where t0 is the start year and t is the end year. E is a list of values (one value per year) and G is defined as:

Where A, alpha and tau are known constants.

Now I need to calculate the result from the convolution integral with the t0 to t limits, but I cannot get it to work. Scipy seems to not have the option to enter the limits.

Anyone that can help?