function probs = sofmax(Qs,tau)

probs = exp(Qs/tau)/sum(exp(Qs/tau));