Курсовая работа: Основные межвидовые взаимодействия, эволюция

Mrav=(a2-b21*a1/b1)/b2;

N0=500;

M0=1500;

[t X]=ode45('siste',[0 20],[N0 M0]);

N=X(:,1);

M=X(:,2);

figure

hold on

grid on

plot(t,N,'b')

plot(t,M,'g')

xlabel('Время');

ylabel('Численности популяций');

legend('1ая популяция','2ая популяция')

hold off

figure

hold on

grid on

if(a2>b21*a1/b1)

plot(N,M,'r')

plot(Nrav,Mrav,'b*')

else

plot(N,M,'')

end

plot(N0,M0,'m*')

xlabel('1ая популяция');

ylabel('2ая популяция');

hold off

a2=4.2: