q=.2, t = .4, a=.6, solve (q)/(4 pi t r^2)(1 - e^(-2 r /a)(2(r/a)^2+2(r/a)+1)) = 1000