src/test/TEST R/covariance.R

   1 covariance= function(p,a){
   2   A = matrix(a, p,p)
   3   for{i in 1:p}{
   4     for{k  in 1:p}{
   5       A[i,]= A[i,]^abs(i-k)
   6     }
   7   }
   8   return(A=A)
   9 }