clear all;
close all;
Cd_sphere;
nu = 0.00000131; %m^2/s
nu = nu/(.0254*.0254); % in^2/s
rho = 64.4; %lbs/ft^3;
rho = 64.4/(12*12*12); %lbs/in^3  
rho = rho/32.174; %slugs/in^3


d = 13; %inches
U = Re*nu/d;% in/s
Drag = 0.5*rho*(pi*d*d/4)*Cd.*(U.^2);
U = U*0.0254/.515; %kts
plot(U,Drag,'bx-');
hold on;
d = 19; %inches
U = Re*nu/d;% in/s
Drag = 0.5*rho*(pi*d*d/4)*Cd.*(U.^2);
U = U*0.0254/.515; %kts
plot(U,Drag,'ro-');
xlim([0 8]);
xlabel('kts');
ylabel('lbs');

legend('13 inch diameter sphere','19 inch diameter sphere',2);

print -depsc2 dragonsphere