> x<-rnorm(100)
> wal_whileloop<-function(x)
+ {
+   n<-length(x)
+   i<-1
+   j<-1
+   k<-1
+   vec<-NULL
+   while(j<=n)
+   {
+     while(i<=j)
+     {
+       vec[k]<-(x[i]+x[j])/2
+       k<-k+1
+       i<-i+1
+     }
+     j<-j+1
+     i<-1
+   }
+   return(median(vec))
+ }
> wal_forloop<-function (x)
+ {
+   n <- length(x)
+   w <- vector(n * (n + 1)/2, mode = "numeric")
+   ind <- 0
+   for (i in 1:n) {
+     for (j in i:n) {
+       ind <- ind + 1
+       w[ind] <- 0.5 * (x[i] + x[j])
+     }
+   }
+   return(median(w))
+ }
> wal_whileloop(x)
[1] 0.03776875
> wal_forloop(x)
[1] 0.03776875