function [pi] = statProb(P)
    Ptmp = P'-eye(size(P));
    Ptmp(end,:) = 1;
    pi = Ptmp\[zeros(200,1);1];
end