finished Example 4 and tested vs. Chris' implementation. Seems to work well. Also Boltzmann potentials implemented