#Fare attenzion al percorso nelle istruzioni che seguono

riverflow=scan("po-pontelagoscuro.txt")
duration=seq(1,32850)
duration=duration/90
riverflow1=sort(riverflow,decreasing=T)
plot(duration,riverflow1,type="l")

#Fare attenzione alla soluzione per utilizzare gli apici nelle label degli assi
par(mar=c(5,5,5,5))
plot(duration,riverflow1,type="l", main="Curva di durata delle portate per il Fiume Po a Pontelagoscuro (1920-2009)",lwd=4, col="red",xlab="Durata (giorni)",ylab=expression('Portata fluviale media giornaliera (m'^3*'/s)'))

riverflow2=array(riverflow,dim=c(365,90))
for (i in 1:90)
riverflow2[,i]=sort(riverflow2[,i],decreasing=T)

matplot(riverflow2,type="l")

avg=rep(0,365)
for (i in 1:365)
avg[i]=mean(riverflow2[i,])

plot(duration,riverflow1,type="l",main="Curva di durata delle portate per il Fiume Po a Pontelagoscuro (1920-2009)",lwd=4, col="red",xlab="Durata (giorni)",ylab=expression('Portata fluviale media giornaliera (m'^3*'/s)'))
lines(avg,lwd=4,col="blue")

standev=rep(0,365)
for (i in 1:365)
standev[i]=sd(riverflow2[i,])

matplot(riverflow2,type="l")
lines(avg,lwd=6,col="red")
lines(avg+standev,lwd=6,col="blue")
lines(avg-standev,lwd=6,col="blue")