Add bell function

[zx8754]

https://github.com/oncogenetics/bell/blob/master/bell.R