rmvbinary_EP function

Simulating correlated binary variables using the algorithm by Emrich and Piedmonte (1991)