integral (a-A exp(-lambda s))^(2*m+1) exp(-epsilon*(t-s))ds