function [b,H] = is_indep(P,m,n)

%[b,A,C,dd] = is_markov(P,m,n);
[hh,Ht,H] = gethhn(P,m,n);
b = (H<=1+1e-10);

return