clear all;
close all;


f = 0.035;

dj = .1:.1:1.0;

do = 1.5;  %inches

lc = 1.0;
th = atan2(do-dj/2,lc);

k = 0.5*(1-(dj/do).^2)*(sin(th).^.5)./(dj/do).^4;
