free-energy code snippets in matlab