This is an IPython notebook showing how to reproduce the solitary waves called stegotons, discovered by LeVeque & Yong in 2003.