%%chemostat diferansiyel denklem çözelim%%%
%%farklı parametre değerleri için%%%
clear
%%ilk değer%%
nn(1)=0.1;%denge noktası 0
cc(1)=0.2;%denge noktası 2
x_v(:,1)=[nn(1); cc(1)];
%%parametreler%%
      alpha1=2;
      alpha2=1;
%%bazı büyüklükler%%
iterasyon=100;
mu=0.1;
%%dif denklem%
for i=1:0.1:3
    alpha1=i+1;
    alpha2=i;
for k=1:iterasyon
    nn(k+1)=nn(k)+mu*(alpha1*(cc(k)/1+cc(k))*nn(k)-nn(k));
    cc(k+1)=cc(k)+mu*(alpha2-(cc(k)/1+cc(k))*nn(k)-cc(k));
    durum_v(:,k)=[nn(k); cc(k)];
end
subplot(3,1,1), plot(durum_v(1,:)),hold on
subplot(3,1,2), plot(durum_v(2,:)),hold on 
subplot(3,1,3), plot(durum_v(1,:),durum_v(2,:)),hold on
end