A method for solving the time-dependent diffusion and transport equations is described in which the flux Φ(r, t) at position r and time t takes the approximate form α(t)ψ(r, t), where α(t) depends solely on time. The treatment includes a heat transfer model, thus enabling temperature and expansion feedback effects to be incorporated into the solution.