NIT=100; [q, parm] = hmm_reest(parm, x, istart, nsamp, NIT);The algorithm prints the total log likelihood at each iteration. At the end, it prints the final state transition matrix and initial probabilities. These should be close to the correct ones.