%%chemostat diferansiyel denklem çözelim%%%
clear
%%ilk değer%%
nn(1)=0.1;
cc(1)=0.2;
x_v(:,1)=[nn(1); cc(1)];
%%parametreler%%
alpha1=3;
alpha2=2;
%%bazı büyüklükler%%
iterasyon=50;
mu=0.1;
%%dif denklem%
for i=1:0.1:4
    nn(1)=i;
    cc(1)=i;
    x_v(:,1)=[nn(1); cc(1)];
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