function v=packdensity(R,C)

v=R*C*pi/(4*(sqrt(3)/2*(R-1)+1)*(C+1/2));
end
