covariance = function(p,a)
{
	A = matrix(a, p,p)
	for (i in 1:p){
		A[i,] = A[i,]^abs(i-(1:p))
	}
	return (A)
}