# carica l'algoritmo isto3D presente qui e quindi esegui quanto segue.
dati <- read.delim("http://macosa.dima.unige.it/om/voci/bivar/bersaglio2.txt",header=FALSE,skip=10,sep =",")
# ho caricato un po' di uscite di "bersagl2.htm"
n <- length(dati$V1); n
#  5990
xy <- c(dati$V1,dati$V1); min(xy); max(xy)
#  -0.9851912  0.9957661
# traccio un istogramma 10*10; i dati vanno da -1 ad 1
nq <- 10; min <- -1; max <- 1; d <- (max-min)/nq
# il dato x cade nella classe floor((x-min)/d)+1
dati2 <- array(0,dim=c(nq,nq))
for (i in 1:n) {j <- dati$V1[i]; k <- dati$V2[i];
j <- floor((j-min)/d)+1; k <- floor((k-min)/d)+1; dati2[j,k] <- dati2[j,k]+1}
par( mai=c(0,0,0.4,0) )
isto3D(dati2,45,-5,1.2)