playRN function

Function to play Royle-Nichols (RN) model