Fish.sci

//Модель хищник-жертва
//a-Коэффициент роста количества травоядных
//b-Коэффициент перенаселённости
//c-Коэффициент смертности хищников
//d-Вероятность съедения жертвы хищником
function fish(a, b, c, d)
//Задаём массив для травоядных и даём начальное значение 0
x=zeros(100,0);
//Задаём массив для хищников и даём начальное значение 0
y=zeros(100,0);
//Начальное количество травоядных
x(1)=100;
//Начальное количество хищников
y(1)=10;
//Задаём цикл
for i=2:100
x(i)=round((a-b*x(i-1))*x(i-1)-d*x(i-1)*y(i-1));
y(i)=round(c*y(i-1)+d*x(i-1)*y(i-1));
end
t=[1:100];
clf;
plot(t, x,'g', t, y,'r');
 set(gca(),"grid",[1 1]);
 xlabel("Дни");
 ylabel("Кол-во рыб");
 legend(['травоядные';'хищники']);
 title(msprintf('Модель Хищник-Жертва a=%f b=%f c=%f d=%f', a, b, c, d));
endfunction