clear all  

% Bei 50kg - 100kg gilt 200N/m
D = 200; % N/m Federkonstante 
s0 = 50; % m Höhe  
L = 20; % m Seillänge
rho = 1.29; % kg/m Dichte der Luft 
A = 0.4; % m² Querschnittsfläche des Springers 
cw = 0.8; % Wiederstandsbeiwert 
v0 = 0; % m/s Anfangsgeschwindigkeit
m = 80; % kg Masse

g = 9.81; % Beschleunigung m/s² 

c=0.5*cw*rho*A/m;
cf = D/m; 
ta = 0; td = 0.1; tn = 120;
t = ta:td:tn;

% Federdämpfung fehlt 

s1 = 'Werte für die Arbeit';
for i = 1 : 3
    if i >= 2 
        m = m *1.1;
    end 
    sim('Bungee_Jump_Sim3')
    s2 = [' für die Masse: ', num2str(m)];
    s3 = [s1, s2]; 
    figure( 'name', s3 )
    subplot(2,1,1); hold on;
    plot(t,Wkin,t,Wsp,t,Wges,t,Whub);grid on;legend('Wkin', 'Wspann', 'Wges', 'Whub') 
    subplot(2,1,2); hold on;
    plot( t, v ); hold on; 
    plot( t, s, 'm'); hold on; grid on; 
end


figure( 'name', 'Energiewerte' )
plot(t,Wkin,t,Wsp,t,Wges,t,Whub);grid on;legend('Wkin', 'Wspann', 'Wges', 'Whub')
figure( 'name', 'Geschwindigkeit und Weg nach der Zeit' )

plot( t, v ); hold on; 
plot( t, s, 'm'); hold on; grid on;