function ShowBlockMarkov()
% Suppose B (pxp) and C(qxq) are stochastic and Bu = u and Cv = v,
% E.g. p = 3 and q = 4...
B = RandomStoch(3);
u = StationaryVec(B);
C = RandomStoch(4);
v = StationaryVec(C);
% This matrix is stochastic..
A = [B(1,1)*C B(1,2)*C B(1,3)*C;...
B(2,1)*C B(2,2)*C B(2,3)*C;...
B(3,1)*C B(3,2)*C B(3,3)*C]; % A = kron(B,C);
% This method for computing A's stationary vector is O((pq)^3)...
w = StationaryVec(A);
% On the otherhand, it can be shown that w is structured
wA = [u(1)*v; u(2)*v; u(3)*v];
clc
format short
disp(' w wA')
disp([w wA])