Thanks for the contribution. There is also an example code by Chris Jackson in this tutorial (example 2.1):