%% DGLHL ist die Funktion, die die rechte Seite der DDE angibt 
%  DGLHLL muss von der Form (t,y,Z) sein.

function X = DelayF31 (t,y,Z)

ylag1 = Z(:,1);   
ylag2 = Z(:,2);
ylag3 = Z(:,3);

%% Parameter fuer Halbleiterlaser 
   
P=0.136;                       % Dimensionlose Pumpstrom              
T=550;                         % Photonverlustrate
k=4.55e-5;                     % Rueckkopllungsstaerke 
a=3.5;                         % Linienbreite Verbesserungsfaktor
Cp=31.49;                      % Rueckkopplungsphase
                 
taus=1.1;
%Gn=1.1*10^-19;
Gn=1.1*10^-21;


%% Die 3 gekoppelte Differentialgleichungen E(t), P(t) und N(t)

X=[ (y(3)*y(1)+k*(ylag1(1))*(cos((ylag2(2))-y(2)-Cp)));...              % E(t)
    (a*y(3)+k*((ylag1(1))/y(1))*sin((ylag2(2))-y(2)-Cp));...            % P(t)
    ((1/T)*(P-y(3)-(1+2*(y(3)))*(abs(y(1))^2)))];                       % N(t)              
           
end