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