a=.072;b=30/74;c=30;
a1= 1.967;b1=.023;c1=.05;


d=100;
d1=d/100;
psi=9*10^(-8);psi1=9*10^(-10)
f=a:b:c;
f1=f.^2;
f2=f.^.5;
f3=f2.^(-1);
e3=(.1*d1).*(a1.*f2+b1.*f+c1.*f3);
hc=10 .^(-e3);
kf=psi*d/.3048;kf1=psi1*d/0.3048;
xf1=kf1.*f1.*hc;
%semilogy(f,hc)
semilogy(f,xf1)
%hold on

xf=kf.*f1.*hc;
%semilogy(f,xf)
