> n <- 1e+5
> y <- rnorm(n,mean=1)
> x <- rexp(n)-2
> mean(y)/sqrt(var(y) + mean(y)^2 + mean(x)^2*var(y)/var(x))
[1] 0.5788539
> cor(x,x*y)
[1] 0.5751899