Code from http://doi.org/10.1016/j.jtbi.2012.03.008