% Ref. to A.Goergen, "Dynamiques Economiques. 
% Solutions de Problème avec Maple et Matlab", 
% Economica, 2006.

% v=0.5;mu=0.03;a=0.05;b=0.01;n=0.01;
% y1s=(b+mu)./a; 
% y2s=1-(mu+n).*v;
% ystar=[y1s y2s]
%
% The system must be parametrized according 
% to the Goodwin Model (1967).

a0=[0.9 0.92];
b0=[0.85 0.87];

[t1,X]=ode45(@goodwin,[0 70],a0);
[t2,Y]=ode45(@goodwin,[0 70],b0);
plot(X(:,1),X(:,2),'b',Y(:,1),Y(:,2),'r');
title('Goodwin model (1967), Evolution of Labor and Profit');


%"Complex and Chaotic Nonlinear Dynamics. 
% Advances in Economics and Finance, 
% Mathematics and Statistics" 
% T.Vialar, Springer 2009.
% Copyright(c).