clear all  

% Bei 50kg - 100kg gilt 200N/m
D = 180; % 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.78; % 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 = 15;
t = ta:td:tn; 
a =9.81;

s1 = 'Values ??for work';
for i = 1 : 3
    if i >= 2 
        m = m *1.1; 
    end
    d = 100; % d=sqrt(m*D);
    sim('Bungee_Jump_Sim2_1')
    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