P = [0    1/4  1/2  0;
     1/2  1/4  0    1/2;
     1/2  1/4  0    0;
     0    1/4  1/2  1/2];

[V, D] = eig(P)

% rescale the eigenvectors to something a human might find:
v1 = V(:, 1);
v1 = v1 / v1(1)

v2 = V(:, 2);
v2 = v2 / v2(end)

v3 = V(:, 3);
v3 = v3 / v3(1)

v4 = V(:, 4);
v4 = v4 / v4(end)

