Example 2.  Consider the closed five sector economy  [Graphics:Images/LeontiefModelMod_gr_80.gif] consisting of say: Agriculture, Construction, Energy, Manufacturing, and Services where the input-output matrix is given by

[Graphics:Images/LeontiefModelMod_gr_81.gif] .   Find the production vector  [Graphics:Images/LeontiefModelMod_gr_82.gif].   

Solution 2.

Enter the matrix  [Graphics:../Images/LeontiefModelMod_gr_83.gif].  

[Graphics:../Images/LeontiefModelMod_gr_84.gif]


[Graphics:../Images/LeontiefModelMod_gr_85.gif]

If  [Graphics:../Images/LeontiefModelMod_gr_86.gif] is an eigenvalue of  [Graphics:../Images/LeontiefModelMod_gr_87.gif]  and  [Graphics:../Images/LeontiefModelMod_gr_88.gif]  an eigen-pair, then the solution to  [Graphics:../Images/LeontiefModelMod_gr_89.gif] will be a multiple of the eigenvector  [Graphics:../Images/LeontiefModelMod_gr_90.gif].
Let us verify that
  [Graphics:../Images/LeontiefModelMod_gr_91.gif] is an eigenvalue of  [Graphics:../Images/LeontiefModelMod_gr_92.gif].  

[Graphics:../Images/LeontiefModelMod_gr_93.gif]


[Graphics:../Images/LeontiefModelMod_gr_94.gif]

The eigen-pair   [Graphics:../Images/LeontiefModelMod_gr_95.gif]  is easily computed using Mathematica's subroutines  Eigenvalues and  Eigenvectors.  

[Graphics:../Images/LeontiefModelMod_gr_96.gif]



[Graphics:../Images/LeontiefModelMod_gr_97.gif]

Since  [Graphics:../Images/LeontiefModelMod_gr_98.gif]  is a solution for any constant  c, we are permitted to choose any multiple of  [Graphics:../Images/LeontiefModelMod_gr_99.gif] we desire for the solution.  For illustration:

[Graphics:../Images/LeontiefModelMod_gr_100.gif]



[Graphics:../Images/LeontiefModelMod_gr_101.gif]

Aside.  Any value of  c  is permitted, it is your choice.  When [Graphics:../Images/LeontiefModelMod_gr_102.gif] we get

[Graphics:../Images/LeontiefModelMod_gr_103.gif]



[Graphics:../Images/LeontiefModelMod_gr_104.gif]

We are done.

How did we get that "nice" vector ?
By using exact arithmetic and the following calculations.

Form the homogeneous system [Graphics:../Images/LeontiefModelMod_gr_105.gif].   
and the augmented matrix   [Graphics:../Images/LeontiefModelMod_gr_106.gif]  and find its reduced row echelon form.

[Graphics:../Images/LeontiefModelMod_gr_107.gif]



[Graphics:../Images/LeontiefModelMod_gr_108.gif]

[Graphics:../Images/LeontiefModelMod_gr_109.gif]

 

 

The equations for this augmented matrix are  

    [Graphics:../Images/LeontiefModelMod_gr_110.gif]  

There is one free variables which we choose to be  [Graphics:../Images/LeontiefModelMod_gr_111.gif].   It is used in computing  [Graphics:../Images/LeontiefModelMod_gr_112.gif].  
Use [Graphics:../Images/LeontiefModelMod_gr_113.gif] and solve the previous equations for  [Graphics:../Images/LeontiefModelMod_gr_114.gif]  

    [Graphics:../Images/LeontiefModelMod_gr_115.gif]  

Get    

    [Graphics:../Images/LeontiefModelMod_gr_116.gif]

The solution vector  [Graphics:../Images/LeontiefModelMod_gr_117.gif]  is

 

[Graphics:../Images/LeontiefModelMod_gr_118.gif]

We are done.

Aside.  We can verify that this is the solution by direct multiplication A X.  This is just for fun !

[Graphics:../Images/LeontiefModelMod_gr_119.gif]



[Graphics:../Images/LeontiefModelMod_gr_120.gif]

We are really done.

Aside.  Iteration can be used to solve for the solution of the closed Leontief model.  If the sum of the production levels is known and a starting vector is given  [Graphics:../Images/LeontiefModelMod_gr_121.gif] then the simple iteration  [Graphics:../Images/LeontiefModelMod_gr_122.gif]  will converge to  [Graphics:../Images/LeontiefModelMod_gr_123.gif].  Because the largest eigenvalue is  [Graphics:../Images/LeontiefModelMod_gr_124.gif],  this is a simplified version of the "power method" for  finding the dominant eigen-pair.  For the example above  [Graphics:../Images/LeontiefModelMod_gr_125.gif] and the following iteration will converge.

[Graphics:../Images/LeontiefModelMod_gr_126.gif]



[Graphics:../Images/LeontiefModelMod_gr_127.gif]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2004