h=hist(state(to+1:n),[1,2,3]); p=h/(n-t0) p= 0.1852 0.3369 0.4779h=hist(state(t0+1:n),[1,2,3]); p=h/(n-t0) p = 0.1852 0.3369 0.4779