mm=md=c();

n=20;
out=rep(1000,10)
for(i in 1:length(out))
{
x=rnorm((n-i))
dc=c(x,out[i:1])
mm[i]=mean(dc)
md[i]=median(dc)
}
cbind(mm,md)