par(mfrow=c(1,2)) ##General code for an lcg lcg <- function(a, c, p, run.length, seed){ x <- rep(0, run.length)