function p = per(h, r)
    p = (h + 2*r)*(pi/4)^(8*h*r/(2*r + h)^2) - 1;
end

