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

Содержимое work5.m

global a1 a2 b11 b22 b12 b21;

a1=10;

a2=0.24;

b11=0.004;

b22=0.0035;

b12=0.042;

b21=0.033;

Mrav=(-a2*b11+b21*a1)/(b22*b11+b21*b12);

Nrav=(a1*b22+b12*a2)/(b22*b11+b21*b12);

N0=900;

M0=800;

[t X]=ode45('sistc',[0 15],[N0 M0]);

N=X(:,1);

M=X(:,2);

figure

hold on

grid on

plot(t,N,'b')

plot(t,M,'g')

xlabel('Время');

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

legend('Жертва','Хищник')

hold off

figure

hold on

grid on

if(b21*a1>a2*b11)

plot(N,M,'r')