%%chemostat diferansiyel denklem çözelim%%%
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=3;
alpha2=2;
%%bazı büyüklükler%%
iterasyon=100;
mu=0.1;
%%dif denklem%
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,:))
subplot(3,1,2), plot(durum_v(2,:))
subplot(3,1,3), plot(durum_v(1,:),durum_v(2,:))