approx_theta_post_s3 function

Find optimal theta of a chromosome using forward-backward; uses a 3-state allele HMM