clear; clf;

Re = logspace(-2,6,10000);

Cd = 24./Re + 2.6*(Re/5)./(1 + (Re/5).^1.52) + ...
    0.411*(Re/263000).^(-7.94)./(1 + (Re/263000).^(-8)) + ...
    Re.^0.8/461000;

loglog(Re,Cd)

xlabel('\rho U D/\mu')
ylabel('C_d')

