%% Dgl definieren 
%wird benutzt um die funktionswerte zu berechen 
function X = Delaydgl(t,y,Z)
ylag1 = Z(:,1);
X = zeros(2,1); %Vektor mit 2 Zeilen und 1 Spalte, 2 ist wegen 2 dgl 


% Physikalische Konstanten 
a= 3.5;              
g= 2.12*10^-4;
V= 1.6;
n0= 4*10^6;
epsilon=3.375*10^-7;
gp= 500;      %  Photonenlebensdauer in sekunden
k= 20;         % k R?ckkopplungsst?rke 
tau= 4.3;               % tau verz?gerungszeit
phi = 0;    % phi die Phase 
I=1.2;
J= I/exp(1);               % J Stromdichte 
Tn= 1.8;           %Tn Lebensdauer der Ladungstr?ger


 
X(1) = (1+1i.*a)./2.* ((g.*(y(2)./V-n0)./(1+epsilon.*abs(y(1))^2))-gp).*y(1)+k.*ylag1(1).*exp(-1i.*phi);
X(2) = J-y(2)./Tn - (g.*(y(2)./V-n0)/(1+epsilon.*abs(y(1))^2)).*abs(y(1))^2;

end