Update prob_of_orbits.m