added a 4pl simulation function in Rcpp
simulation function for 4pl incl example