%% Halbleiterlaser mit zeitverzoegerte optische Rueckkopplung 

tau=121;                         % Verzoegerungszeit
lags=[tau,tau,tau];              % Vektor, der die Zeitverzoegerung darstellt                 
tspan=[0,1000];                  % Integrationsintervall als zwei-elementiger Vektor
tau_p=2;                         % Zeit t ist skaliert auf tau_p -> laenger integrieren

%% Ausgabe um die Funktionen mit dde23 numerisch zu berechen 
sol = dde23(@DelayF31,lags,@delayhist31,tspan);
 
%% Grafen fuer die Darstellung der 3 Differentialgleichungen 

% Grafen fuer E(t)

subplot(1,3,1) 
plot((((sol.x)/tau_p)),sol.y(1,:),'blue');
title('E(t)');
xlabel('t [ps]');
 
% Grafen fuer P(t) 
subplot(1,3,2);
plot(((sol.x)/tau_p),sol.y(2,:),'magenta');
title('P(t)');
xlabel('t [ts]');
 
% Grafen fuer N(t) 
subplot(1,3,3);
plot(((sol.x)/tau_p),sol.y(3,:),'cyan');
title('N(t)');
xlabel('t [ps]');

% Ausgabe der Ergebnisse 
sol.x
sol.y