sim_ev function

Expected values from state probabilities