R code for RC_KP model