Perform smoothing using the forward-backward algorithm

hmm_forward_backward(ys, x0, transition_matrix, observation)

Arguments

ys	The observed values
x0	the distribution of the initial state
transition_matrix	the transition matrix
observation	the observation distribution