figure();

%% PARAMETER fuer Halbleiterlaser 
Gn=2.142e-5;        %Gewinnfunktion G(N)=Gn(N(t)-Ntr), Gn Differentieller optischer Gewinn in [ns] 
Ntr=1.54e8;         %Trasparenzladungstr?gerzahl 
ggamma=0.357e3;     %Photonenverlustrate in [ns] 
kgamma=0.909;       %Ladungstraegerverlustrate in [ns] 

%% Festen Pumpstromwert 
J=2.32e8;           %Pumpstrom in [ns] 


%Loesung der Dgl mit ode45 
 
  [t,y]=ode45(@(t,y)[((Gn*(y(2)-Ntr)-ggamma)*y(1));...      % E(t)
     J-kgamma*y(2)-(Gn*(y(2)-Ntr))*y(1)],...                % N(t)
     [0;4],[1e0;1e0]);

%% Bild 
 subplot(1,2,1)  
 plot(t,y(:,1),'blue')
 title('P(t)')
 xlabel('Zeit t in [ns]');
% 
 subplot(1,2,2)
 plot(t,y(:,2),'magenta')
 title('N(t)')
 xlabel('Zeit t in [ns]');
 
 
 sol.x 
 sol.y