A = sparse(6,6);
A(1,1) = 1;
A(2,[1,3]) = 1;
A(3,5) = 1;
A(4,[2,3,5]) = 1;
A(5,6) = 1; A(6,5) = 1;
n = 6;

d = sum(A,2);
P = diag(1./d)*A;

% show normout












%% Backup
Pt = P';

v = ones(n,1)/n;

alpha = 0.85;

x = (speye(6) - alpha*Pt)\((1-alpha)*v)

sum(x)