Карибский кризис - Matlab - Ответ 7380155

function du = Richard(t,u)
x = u(1);
y = u(2);
 
%эти коэффициенты описывают положительную обратную связь
%чем большим количеством располагает y, 
%тем больше захочет приобрести x в ответ на воспринимаемую им угрозу
k = 2;
m = 0.8;
 
%эти коэффициенты описывают отрицательную обратную связь
%государство x вынуждено решать социальные задачи и не
%может всю свою экономику свести к производству.
%поэтому, чем большим количеством располагает x,
%тем меньше дополнительных оно сможет
%приобрести из-за существующего бремени расходов
a = 0.02;
b=2;
 
%константы g и h – положительные или отрицательные в
%зависимости от того, насколько в целом враждебно или
%дружественно настроены x и y по отношению друг к другу.
g = 0.0002;
h=3;
 
du = zeros(2,1);
du(1) = k*y - a*x + g;
du(2) = m*x - b*y + h;

%модель Ричардсона
Tk = [1, 5];
u0 = [500, 500]; 
[T, M] = ode45('Richard',Tk,u0); 
X = M(:,1); 
Y = M(:,2); 
plotyy(T,X, T,Y), grid on
legend('Государство А','Государство Б')
ylabel('Военные расходы')
xlabel('Время, года')