clear all
close all

di = .028;
do = .03:.005:.06;

g = 9.81;  %m/s
rho_w = 1025;  %kg/m^3
sg = 0.7:.1:.9;
rho_s = sg*rho_w;

A = pi*(do.^2-di.^2)/4;


Fb = -(rho_s-rho_w)'*A;  %kg/m


plot(do*1000,Fb)

xlabel('do (mm)');
ylabel('Fb (kg/m)');
legend('SG = .7','SG = .8','SG = .9');


print -depsc2 FloatationRequirements