dx^dy^dz == volume(3)
p <- 1
for(i in 1:7){p <- p ^ as.kform(i)}
p
p == volume(7) # should be TRUE
o <- volume(5)
M <- matrix(runif(25),5,5)
det(M) - as.function(o)(M) # should be zero
is.volume(d(1) ^ d(2) ^ d(3) ^ d(4))
is.volume(d(1:9))
Run the code above in your browser using DataLab