Earth Mover's distance

henon <- function(a,b,x){
	c(x[2]+1-a*x[1]^2,b*x[1])
}

henon.series <- function(a,b,x,t){
	ret <- matrix(0,t+1,2)
	ret[1,] <- x
	for(i in 1:t){
		ret[i+1,] <- henon(a,b,ret[i,])
	}
	ret
}

a <- 1.4
b <- 0.3
#x <- c((sqrt(609)-7)/28,3*(sqrt(609)-7)/280)
x <- c(0,0)
t <- 100000

h <- henon.series(a,b,x,t)
plot(h,cex=0.001)