function [K]=GlabCorrectionFactor(b_dRatio);
%[K]=GlabCorrectionFactor(b_dRatio);
% we can write the correction factor for G as:
% K = (1 - b^3/dp^3)^-1 w/ dp = sqrt(4d^2+b^2)
% or, in terms of b/d as:
% K = 1./(1 - (1./sqrt(1+4./b_dRatio.^2)).^3) w/ b_dRatio = b/d

K = 1./(1 - (1./sqrt(1+4./b_dRatio.^2)).^3)

clf
plot(b_dRatio,K);


return