close all

Ta = 30;    % Umgebungstemp in K
g_r = 13.6; % Kupfersprung

R = 10.557; % Pk/P0
u = 1.3;    % 
v = 0.8;    % 
k11 = 0.5;

dT_Or = 46;     % K

K = [190 1
     365 0.6
     500 1.5
     710 0.3
     735 2.1
     750 0];  
anzahl = size(K); % Ermittlung der Zeilen von K
step_init = 0 % Initalwert des Sprungs; in der For-Schleife, wird Initalwert auf den Vergangenheitswert gesetzt
for i=1:anzahl(1,1) % Durchlauf durch Matrix bis zur letzten Zeile --> in anzahl (max_Zeilen, max_Spalten)
    
    if (t<K(i,1)) %während t kleiner als t2
        step_final = K(i,2);
        
        
       
    end
    if (i>1)
        step_init = K(i-1,2); % setzten von Initalwert, wenn i>1
    end
end
sim_time = K(anzahl(1,1),1);   
%[t,x,y] = sim('Loesung_mit_DGL_06',450); 
[t,x,y] = sim('test_03',sim_time);
%simulation bis zum großten t von Matrix K
plot(t,y(:,2));


hold on