Two-dimensional Barenblatt solution with diffusivity proportional to concentration (source code)

= Two-dimensional Barenblatt solution with diffusivity proportional to concentration
{c}

For $C_t=k\nabla\cdot(C\nabla C)$ with total mass $2\pi M$,
$$
C(r,t)=\sqrt{\frac{M}{kt}}
\left(\frac1{\sqrt2}-\frac18\frac{r^2}{\sqrt{Mkt}}\right)_+.
$$
Its compact support has radius $r_0(t)=(32Mkt)^{1/4}$.