jubilee.optimal_tb3ms function

Calculate the optimal TB3MS