clear all
close all



I = zeros(3);
I(1,1) = 1613.25; I(2,2) = 1615.53; I(3,3) = 1613.25;

Buoy = RigidBody;
Buoy = set(Buoy,'mass',3220.13,'S',3.14159,'disp',3.14159,'g',9.81);
Buoy = set(Buoy,'xG',[0 0 -.25]);
Buoy = set(Buoy,'xB',[0 0 -0.5]);
Buoy = set(Buoy,'S11',0.79,'S22',0.79);  %m^4
Buoy = set(Buoy,'I',I);


x = [-.1:.01:.1];


for k = 1:length(x)
 Buoy = set(Buoy,'sv',x(k)*[0 0 1 0 0 0 0 0 0 0 0 0]);
 F_HS(k,:) = F_Hydrostatic(Buoy,k);  
end


plot(x,F_HS);
legend('x1','x2','x3','x4','x5','x6');

