# coordinate geografiche, densità di popolazione
# es.: 0.1268 pop. in milioni, area 3263 in km^2
dati1 <- read.table("http://macosa.dima.unige.it/om/prg/R/rmacosa/penisola.txt",header=FALSE,sep =",")
dati2 <- read.table("http://macosa.dima.unige.it/om/prg/R/rmacosa/sicilia.txt",header=FALSE,sep =",")
dati3 <- read.table("http://macosa.dima.unige.it/om/prg/R/rmacosa/sardegna.txt",header=FALSE,sep =",")
mo <- min(dati1$V1,dati2$V1,dati3$V1);Mo <- max(dati1$V1,dati2$V1,dati3$V1)
mv <- min(dati1$V2,dati2$V2,dati3$V2);Mv <- max(dati1$V2,dati2$V2,dati3$V2)
par( mai = c(0.5,0.5,0.1,0.1) )
plot(c(mo,Mo),c(mv,Mv),type="n",xlab="", ylab="",asp=1)
abline(h=seq(37,47,1),v=seq(7,18,1),lty=3,col="grey")
lines(dati1,col="black",lty=3)
lines(dati2,col="brown",lty=3)
lines(dati3,col="blue",lty=3)
cerchio <- function(x,cit) {r <- sqrt(x[3]*10^6/x[4])/25
  symbols(x[1],x[2],circles=r,inches=FALSE,add=TRUE,fg="red")
  text(x[1],x[2],cit,font=2) }
aosta <- c(7.31, 45.75, 0.1268, 3263)
trentino <- c(11.12, 46.08, 1.0433, 13607)
friuli <- c(13.23, 46.07, 1.060, 8240)
piemonte <- c(7.68, 45.08, 4.378, 25402)
lombardia <- c(9.17, 45.47, 9.848, 23860)
veneto <- c(12.33, 45.43, 4.893, 18399)
liguria <- c(8.93, 44.42, 1.567, 5420)
toscana <- c(11.24, 43.78, 3.692, 22994)
e_romagna <- c(11.34, 44.5, 4.379, 22447)
umbria <- c(12.39, 43.11, 0.8895, 8456)
marche <- c(13.5, 43.6, 1.542, 9366)
lazio <- c(12.5, 41.89, 5.566, 17236)
abruzzo <- c(13.39, 42.37, 1.326, 10753)
molise <- c(14.67, 41.57, 0.313, 4438) 
campania <- c(14.25, 40.83, 5.765, 13595)
puglia <- c(16.38, 41.06, 4.050, 19358)
basilicata <- c(15.81, 40.64, 0.575, 9992)
calabria <- c(16.60, 38.89, 1.957, 15079)
sicilia <- c(14.14, 37.56, 4.996, 25711)
sardegna <- c(9.08, 40.05, 1.641, 24090)
cerchio(aosta,"ao"); cerchio(trentino,"tr"); cerchio(piemonte,"pi")
cerchio(veneto,"ve"); cerchio(lombardia,"lo"); cerchio(liguria,"li")
cerchio(toscana,"to"); cerchio(e_romagna,"er"); cerchio(umbria,"um")
cerchio(marche,"ma"); cerchio(lazio,"la"); cerchio(friuli,"fr")
cerchio(abruzzo,"ab"); cerchio(molise,"mo"); cerchio(campania,"ca")
cerchio(puglia,"pu"); cerchio(basilicata,"ba"); cerchio(calabria,"ca")
cerchio(sicilia,"si"); cerchio(sardegna,"sa")
text(16,44.5,"Densità di popolazione",font=2)