#----------------------------------------------------------------------
# Esperimento di affondamenti in acqua (FisAI 1/3/2023)
# 
#----------------------------------------------------------------------

pausa <- function() { cat ("\n >> press Enter to continue\n"); scan() }

# leggi i dati dal file
dati <- read.table("misure_affondamenti.dat", header=TRUE)

# mettiamoli in comodi 'vettori', 
# sottraendo ai valori la lettura prima dell'immersione
x  <- dati$x  - dati$x[1]
y1 <- dati$y1 - dati$y1[1]
y2 <- dati$y2 - dati$y2[1]
y3 <- dati$y3 - dati$y3[1]
y4 <- dati$y4 - dati$y4[1]
n  <- length(x)

print(dati)
plot(dati, col='blue', pch=19)

pausa()

plot(x, y1, xlab='affondamento (cm)', ty='b',
     ylab="massa d'acqua spostata (g)",
     pch=1, cex=1.5, col='darkblue',
     main="Massa d'acqua spostata in funzione dell'affondamento")
grid()
points(x, y2, pch=1, ty='b', cex=1.5, col='blue')
points(x, y3, pch=17, ty='b', cex=1.5, col='orange')
points(x, y4, pch=19, ty='b', cex=1.5, col='darkgray')
legend('topleft', legend=c("Cilindro 1", "Cilindro 2", "Prisma", "Cono"),
       lty=1, cex=1.5,
       col=c("darkblue", "blue","orange", "darkgray"),
       text.col=c("darkblue", "blue","orange", "darkgray"),
       box.col='gray')

pausa()
plot(x[2:n], y1[2:n], xlab='affondamento (cm)', ty='b',
     ylab="massa d'acqua spostata (g)", ylim=c(0.04,20),
     pch=1, cex=1.5, col='darkblue', log='xy',  
     main="Massa d'acqua spostata in funzione dell'affondamento -- scale logaritmiche")
grid()
points(x[2:n], y2[2:n], pch=1, ty='b', cex=1.5, col='blue')
points(x[2:n], y3[2:n], pch=17, ty='b', cex=1.5, col='orange')
points(x[2:n], y4[2:n], pch=19, ty='b', cex=1.5, col='darkgray')
legend('bottomright', legend=c("Cilindro 1", "Cilindro 2", "Prisma", "Cono"),
       lty=1, cex=1.5,
       col=c("darkblue", "blue","orange", "darkgray"),
       text.col=c("darkblue", "blue","orange", "darkgray"),
       box.col='gray')