Example 5.  Let  [Graphics:Images/PowerMethodMod_gr_319.gif]  record the number of people in a certain city who use brands X, Y, and Z, respectively.  
Each month people decide to keep using the same brand or switch brands.  
The probability that a user of brand
X will switch to brand Y or Z is 0.4 and 0.2, respectively.
The probability that a user of brand
Y will switch to brand X or Z is 0.3 and 0.2, respectively.
The probability that a user of brand
Z will switch to brand X or Y is 0.1 and 0.3, respectively.
The transition matrix for this process is  
[Graphics:Images/PowerMethodMod_gr_320.gif] or
        [Graphics:Images/PowerMethodMod_gr_321.gif]
Assume that the initial distribution [Graphics:Images/PowerMethodMod_gr_322.gif].
5 (a).  Find the first few terms in the sequence [Graphics:Images/PowerMethodMod_gr_323.gif].  
5 (b).  Verify that  [Graphics:Images/PowerMethodMod_gr_324.gif] is the dominant eigenvector of  A.  
5 (c).  Verify that a corresponding eigenvector is  [Graphics:Images/PowerMethodMod_gr_325.gif].
5 (d).  Conclude that the limiting distribution is  [Graphics:Images/PowerMethodMod_gr_326.gif].  

Solution 5.

5 (a).  Enter the matrix A and vector [Graphics:../Images/PowerMethodMod_gr_327.gif] and use the subroutine Markov to find the first few terms in the sequence [Graphics:../Images/PowerMethodMod_gr_328.gif].  

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



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

[Graphics:../Images/PowerMethodMod_gr_331.gif]
[Graphics:../Images/PowerMethodMod_gr_332.gif]
[Graphics:../Images/PowerMethodMod_gr_333.gif]
[Graphics:../Images/PowerMethodMod_gr_334.gif]
[Graphics:../Images/PowerMethodMod_gr_335.gif]
[Graphics:../Images/PowerMethodMod_gr_336.gif]
[Graphics:../Images/PowerMethodMod_gr_337.gif]
[Graphics:../Images/PowerMethodMod_gr_338.gif]
[Graphics:../Images/PowerMethodMod_gr_339.gif]
[Graphics:../Images/PowerMethodMod_gr_340.gif]
[Graphics:../Images/PowerMethodMod_gr_341.gif]
[Graphics:../Images/PowerMethodMod_gr_342.gif]
[Graphics:../Images/PowerMethodMod_gr_343.gif]

5 (b).  Verify that  [Graphics:../Images/PowerMethodMod_gr_344.gif] is an eigenvector of  A.  

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


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

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

5 (c).  Verify that  [Graphics:../Images/PowerMethodMod_gr_348.gif] is an eigenvector of  A  and a corresponding eigenvector is  [Graphics:../Images/PowerMethodMod_gr_349.gif].

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



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

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

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

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

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

5 (d).  The iteration in part (a) appears to be converging to  [Graphics:../Images/PowerMethodMod_gr_356.gif].

We are done.

Aside.  We can graph the situation, this is just for fun !

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


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

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



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


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

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2004