approx_theta_post_s2 function

Find optimal theta of a CNV segment using forward-backward; uses a 2-state allele HMM