Hi team,

After a brief chat with @fonnesbeck in Southampton, UK a few months back, I thought I would try using the `pm.Potential`

function to help teach MCMC in physics and it seems to work excellently.

I’d be keen to write up a notebook for the documentation if you think it would fit?

Here’s a notebook I produced for our group link. It would need adapting a little but I think an example of the Heisenberg model would be a nice example of how pm.Potential can be used.

I think the important concepts for those looking to use PyMC3 for physics problems would be:

- Showing how energy and the probability distribution link
- Explaining the necessity of the prior distributions
- Plotting the energy and probability surfaces

Any other suggestions?