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

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

pausa()

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

plot(x, y1, xlab='x', ty='b', ylab='y',
     pch=1, cex=1.5, col='darkblue',
     main="Variazione di letture sulla bilancia 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, y1, xlab='x (cm)', ty='b', ylab='variazione di lettura (g)',
     pch=1, cex=1.5, col='darkblue',
     main="Variazione di letture sulla bilancia 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, 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')