\begin{Schunk}
\begin{Sinput}
 x <- rnorm(100, m=50,s=10)
 y <- 2 + 3 * x + 80 * rnorm(100)
 plot(x,y, type="n",xlab="Indep. Var.", ylab="Dep. Var",ylim=c(0,350))
 for(i in 1:100){
    y <- 2 + 3 * x + rnorm(100, m=0, s=80)
    mod <- lm (y ~ x)
    abline(mod,col=gray(.80))
  }
 points(x,y)
\end{Sinput}
\end{Schunk}
