wolffd@0: function p = is_stochastic(T)
wolffd@0: % IS_STOCHASTIC Is the argument a stochastic matrix, i.e., the sum over the last dimension is 1.
wolffd@0: % p = is_stochastic(T)
wolffd@0: 
wolffd@0: p = approxeq(T, mk_stochastic(T));