optical depth integral
from sage.symbolic.integration.integral import indefinite_integral t,D,b = var('t D b') indefinite_integral(exp(D*t)*t^(4*b), t)
n(gamma (3,-2),20) n(gamma (3,2),20)