This code is from an Insight Journal publication that can be downloaded here: http://hdl.handle.net/10380/3105
 
