---------- ---------- ---------- ---------- ---------- ---------- ---------- ---------- # 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")