clc
clear all
% ДИСКРЕТНАЯ МОДЕЛЬ ЭКОНОМИЧЕСКИЙ ОСЦИЛЛЯТОР САМУЭЛЬСОНА
%ЗАДАНИЕ ИСХОДНЫХ ДАННЫХ  
%горизонт прогноза
T=50;
%Норма потребления
nk=0.6;
% Акселератор
b=1.57;
%%начальное значение 
for j=1:1:T
 Y(j)=0;
 C(j)=0;
  Inv(j)=0; 
end
Inv(T+1)=0;
Y(T+1)=0;
%Расчет эндогенных переменных модели для t =1
t=1;
Y(t)=20.0;
C(t)=10.0;
%Расчет эндогенных переменных модели для t=2 
%доход
%потребление
%инвестиции
t=2;
C(t)=nk*Y(t-1);
%инвестиции
Inv(t)=b*(C(t)-C(t-1));
%ДОХОД
Y(t)= C(t)+Inv(t);
for  t=3:1:T
 t;
C(t)=nk*Y(t-1);
C(t-1)=nk*Y(t-2);
Inv(t)=b*(C(t)-C(t-1));
Y(t)=C(t)+ Inv(t);
end
Y(T+1)=[];
Inv(T+1)=[];
X=1:1:T;
plot(X,Y, X,C,X,Inv)
title('ОСЦИЛЯТОР САМУЭЛЬСОНА Y(t) C(t)  Inv(t)','FontName','Arial Unicode MS','FontSize',12,'FontWeight','bold')
xlabel('t-ИНТЕРВАЛ ВРЕМЕНИ','FontName','Arial Unicode MS','FontSize',12,'FontWeight','bold' )  
ylabel('ДОХОД  ПОТРЕБЛЕНИЕ  ИНВЕСТИЦИИ','FontName','Arial Unicode MS','FontSize',10,'FontWeight','bold')