Single particle diffusion in a 2D unit square.
Source of N particles per unit time wiht a 2D unit square.
... with Robin boundary conditions; analytic solution via eigenfunction expansion. Density plots at multiple times. Parameters D (diffusion constant) and K (reflective ratio) are tunable.