Update to iterations, now use proper diffusion equations.  Current run shows a transition to a linear density function, as would be expected between constant source and perfect sink boundary conditions.