t = symbols('t', real=True, positive=True) a = symbols('a', real=True) s = symbols('s') laplace_transform(1 + 2*Sum((-1)**n*Heaviside(t - n*a), (n, 1, oo)), t, s)

returns: (2*LaplaceTransform(Sum((-1)**n*Heaviside(-a*n + t), (n, 1, oo)), t, s) + 1/s, 0, True)

Expected result should be: 1/s + 2*Sum( (-1)**n*exp(-s*n*a)/s, (n, 1, oo))

So it just commutes with the sum?

We have moved issues to GitHub https://github.com/sympy/sympy/issues.

Migrated to http://github.com/sympy/sympy/issues/7177