ita <- read.table("http://macosa.dima.unige.it/om/prg/R/rmacosa/penisola.txt",header=FALSE,sep =",")
# il grafico
plot(ita,asp=1)
lines(ita)
abline(h=seq(37,47,1),v=seq(7,18,1),lty=3,col="grey")
x <- ita$V1; y <- ita$V2
# l'area usando come unità di misura i "quadri" delimitati da
# longitudini e latitudini
n <- length(x); area <- (y[n]+y[1])*(x[n]-x[1])
for (i in 1:(n-1)) area <- area + (y[i]+y[i+1])*(x[i]-x[i+1]); area <- area/2
area
# la sua trasformazione in km^2 tenendo conto che 
# 1 grado di latitudine è 111 km (valore arrotondato)
# 1 grado di longitudine in Italia è mediamente circa 79 km
area*111*79; round(area*111*79,-3)
# l'area della penisola italiana è effettivamente 252000 km^2