varies with radius r as ϕ={ϕ 0r^2/R^3