\begin{Schunk}
\begin{Sinput}
 b0 <- 3; b1 <- 0.25; stde <- 0.0
 x <- rnorm(100, m=50, s=10)
 x <- x[order(x)]
 y <- b0 + b1 * x
 rx <- range(x)
 ry <- range(y)
 plot(y ~ x, type="p", xlab="x: the IV", ylab="y: the DV")
 text(.8*rx[2],ry[2],pos=1, paste("Cov[x,y]=",round(cov(x,y),2),"\n r=",round(cor(x,y),2)))
\end{Sinput}
\end{Schunk}
