Example 4.  Let  [Graphics:Images/PowerMethodMod_gr_276.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.3 and 0.3, 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_277.gif] or
        [Graphics:Images/PowerMethodMod_gr_278.gif]
Assume that the initial distribution [Graphics:Images/PowerMethodMod_gr_279.gif].
4 (a).  Find the first few terms in the sequence [Graphics:Images/PowerMethodMod_gr_280.gif].  
4 (b).  Verify that  [Graphics:Images/PowerMethodMod_gr_281.gif] is the dominant eigenvector of  A.  
4 (c).  Verify that a corresponding eigenvector is  [Graphics:Images/PowerMethodMod_gr_282.gif].
4 (d).  Conclude that the limiting distribution is  [Graphics:Images/PowerMethodMod_gr_283.gif].  

Solution 4.

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

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



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

[Graphics:../Images/PowerMethodMod_gr_288.gif]
[Graphics:../Images/PowerMethodMod_gr_289.gif]
[Graphics:../Images/PowerMethodMod_gr_290.gif]
[Graphics:../Images/PowerMethodMod_gr_291.gif]
[Graphics:../Images/PowerMethodMod_gr_292.gif]
[Graphics:../Images/PowerMethodMod_gr_293.gif]
[Graphics:../Images/PowerMethodMod_gr_294.gif]
[Graphics:../Images/PowerMethodMod_gr_295.gif]
[Graphics:../Images/PowerMethodMod_gr_296.gif]
[Graphics:../Images/PowerMethodMod_gr_297.gif]
[Graphics:../Images/PowerMethodMod_gr_298.gif]
[Graphics:../Images/PowerMethodMod_gr_299.gif]

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

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


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

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

4 (c).  Verify that  [Graphics:../Images/PowerMethodMod_gr_304.gif] is an eigenvector of  A  and a corresponding eigenvector is  [Graphics:../Images/PowerMethodMod_gr_305.gif].

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



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

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

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

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

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

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

We are done.

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

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


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

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



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


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

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2004