---------- ---------- ---------- ---------- ---------- ---------- ---------- ----------
# millimetri di pioggia mensile per tre anni, in una data località

mm <- c(44, 59, 75, 80,96,98,97,99,100,102,104,101,
        103,102,108,100,92,71,49,28, 27, 45,66, 81,
         92,99, 93, 85, 75,67,61,62, 64, 68,71, 72)
HF=3.5; BF=3.5; boxmN(-110,110, -110,110)
for(i in 1:4) cerch(0,0, i*25, "grey60")
for(i in 1:12) Direzio(0,0, i*30, 100, "grey60")
ang <- seq(30,30*12*3,30) 
spezzata(xrot(ang)*mm,yrot(ang)*mm,"blue")
A1 <- ang[1:12]; A2 <- ang[13:24]; A3 <- ang[25:36]
mm1 <- mm[1:12]; mm2 <- mm[13:24]; mm3 <- mm[25:36]
puntiO(xrot(A1)*mm1,yrot(A1)*mm1,"black")
puntiD(xrot(A2)*mm2,yrot(A2)*mm2,"black")
puntiV(xrot(A3)*mm3,yrot(A3)*mm3,"black")
text(xrot(ang[1]+5)*115,yrot(ang[1]+5)*115,"Ja")
text(xrot(ang[4]+5)*115,yrot(ang[4]+5)*115,"Ap")
text(xrot(ang[7]+5)*115,yrot(ang[7]+5)*115,"Ju")
text(xrot(ang[10]+5)*115,yrot(ang[10]+5)*115,"Oc")

    

# coordinate geografiche, densità di popolazione di alcune regioni
# es.: 0.117208*10^6 popolaz., area 3262 in km^2   (anno: 1991)
# (ogni disco rappresenta tre dati: la densità di popolazione e le coordinate
#  geografiche del capoluogo)
CERCHIO = function(x,cit) {r <- sqrt(x[3])/25;
  cerch(x[1],x[2],r,"red"); text(x[1],x[2],cit,font=2) }
PIANO(7,14, 41.3,46.5)
aosta <- c(7.31,45.75,0.117208*10^6/3262)
trentino <- c(11.12,46.08,0.934731*10^6/13618)
piemonte <- c(7.68,45.08,4.338*10^6/25399)
lombardia <- c(9.17,45.47,8.941*10^6/23859)
veneto <- c(12.33,45.43,4.453*10^6/18364)
liguria <- c(8.93,44.42,1.702*10^6/5418)
toscana <- c(11.24,43.78,3.599*10^6/22992)
e_romagna <- c(11.34,44.5,3.984*10^6/22123)
umbria <- c(12.39,43.11,0.822972*10^6/8456)
marche <- c(13.5,43.6,1.447*10^6/9693)
lazio <- c(12.5,41.89,5.146*10^6/17203)
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")

# Usando il generatore di numeri casuali runif(n,A,B) che genera n numberi casuali distri-
# buiti uniformemente tra A e B (vedi qui) posso confrontare le densità di popolazione
BF=2.5; HF=2.5
m=-0.5; M=10.5
PIANOss(m,M, m,M)
spezza(c(m,M,M,m,m), c(m,m,M,M,m),"brown")
aosta = 0.117208*10^6/3262; aosta
# 35.93133
liguria = 1.702*10^6/5418; liguria
# 314.1381
for(i in 1:round(aosta) ) Punto(runif(1,0,10),runif(1,0,10),"blue")
sottoassex("val d'aosta")
PIANOss(m,M, m,M)
spezza(c(m,M,M,m,m), c(m,m,M,M,m),"brown")
for(i in 1:round(liguria) ) Punto(runif(1,0,10),runif(1,0,10),"blue")
sottoassex("liguria")