Power iteration in Matlab.

How can I implement the power iteration method in matlab? I have a vector, v and a matrix, G, I need to multiply them together using the answer in the previous iterate continuously until it converges to a constant vector. So about 20-30 iterations will need to be done.

   v = 0.2669    0.1413    0.3570    0.3836    0.5105    0.6112
G =

    0.1667    0.1667    0.1667    0.1667    0.1667    0.1667
    0.8750    0.0250    0.0250    0.0250    0.0250    0.0250
    0.8750    0.0250    0.0250    0.0250    0.0250    0.0250
    0.3083    0.5917    0.0250    0.0250    0.0250    0.0250
    0.8750    0.0250    0.0250    0.0250    0.0250    0.0250
    0.3083    0.3083    0.3083    0.0250    0.0250    0.0250

Answer 1

[Walter Roberson]

r = v;
for K = 1 : 30
  r = r * G;
end

[Torsten]

Note that the vector r in the power iteration method must be normalized in each step:

r = r/norm(r);

Best wishes

Torsten.

Answer 2

[Bolivia Correa Tayniz]