Comment by edp

This was code from another article, linked in the introduction.