Integrate[(1 - r / R)^2 r^3, {r, 0, R}]