wclear all;
clf;         %Figure Window lschen
clc;         %Eintrge in Command Window lschen 
%hold off;   %Einschalten des normalen Lschmodus
%hold on;    %hold on:  "halten" der letzten dargestellten Kurve 
             %>> berlagerung der einzelnen Kurven

%% Hauptabmessungen
Hb=75;
Tb=33.75;
Rb=57.5;
rho=1025; %[kg/m]
g=9.81;
a=18; %[m] Abstand von CL zur Kette auf der BL
alpha=pi/4; %Seilwinkel zum Boden
A_WL=pi*Rb^2;   %Wasserlinienflche
zb=-Tb/2;
zg=-Hb/2;
b=-zg;

%% Hydrodyn. Massen
A22=rho*A_WL*Tb;
A24=-rho*A_WL*Tb*(zb-zg);
A33=2.5*4/3*rho*Rb^3;
%A33=4/3*rho*Rb^3;
A44=rho*A_WL*((Tb^3)/12+Tb*(zb-zg)^2);
A42=A24;
%% Masse der verdrngten Flssigkeit
M_0=rho*pi*Rb^2*Tb;
%% Masse Zylinder
M_0_zyl=296500000;
%% Trgheitsmoment Zylinder
I44=0.5*M_0_zyl*Rb^2;
%% Gesamtmasse
M22=M_0+A22;
M33=M_0+A33;         
M44=I44+A44;
%% Rckstellkrfte
C22=cos(alpha)*x;
C33=rho*g*A_WL+sin(alpha2*x);      
C44=rho*g*(A_WL*Tb*(zb-zg)+I44)+sin(alpha)*x*a-cos(alpha)*x*b;
%% Eigenkreisfrequenz
w22_0=sqrt(C22/(M22_ges+A24));
w33_0=sqrt(C33/M33_ges);
w44_0=sqrt(C44/(M44_ges+A42));
%% Dmpfung
delta=0.1;               %Dmpfungsverhltnis
B22=delta*w22_0*2*M22_ges; 
B33=delta*w33_0*2*M33_ges; 
B44=delta*w44_0*2*M44_ges; 
%% Ausgabe im Command Window
disp(sprintf('M22=  %g',M_0));
disp(sprintf('M33=  %g',M_0));
disp(sprintf('I44=  %g\n',I44));
disp(sprintf('A22=  %g',A22));
disp(sprintf('A24=  %g',A24));
disp(sprintf('A33=  %g',A33));
disp(sprintf('A44=  %g',A44));
disp(sprintf('A42=  %g\n',A42));
disp(sprintf('B22=  %g',B22));
disp(sprintf('B33=  %g',B33));
disp(sprintf('B44=  %g\n',B44));
disp(sprintf('C22=  %g',C22));
disp(sprintf('C33=  %g',C33));
disp(sprintf('C44=  %g',C44));


    F2=1166.38   
  
    F3=19.8    
   
    F4=-(53.569-19.458)*5.95+(-52.569+72.154)*19.13-(52.569-33.75)*1141.3
          
end;
M=[M_22 A_24;A_24 M_44]
B=[B22 B44]
C=[C22 0; 0 C44]
F=[F2 F4]


  end