expand(integral exp(-(k+1)T/tau)*exp(t/tau)*(x1 + (t-kT)*(x2-x1)/T)/tau dt from kT to (k+1)T)