%% 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