## Copyright Ton Dieker
## Centre of Mathematics and Computer Science (CWI) Amsterdam
## April 2002

## ton@cwi.nl
	
cov<-function(N, left, right, H) {
dyn.load("rmdcov.so")
result<-rep(0,2^(2*N))
output<-.C("rmdcov", as.integer(N),
		     as.integer(left),
		     as.integer(right),
		     as.double(H),
		     result = as.double(result))
matrix(output$result,2^N,2^N)
}
