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