a <- 10 b <- 5 x <- 8 * pi * runif(1000) delta <- rnorm(1000, sd = 2) y <- a * sin(x) + b + delta plot(x, y)