#----------------------------------------------------------------------
# Pre-analisi esperimenti simulati
#  -> plot in scala log-log
#
# GdA 21/03/2024
#----------------------------------------------------------------------

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

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

# mettiamoli in comodi 'vettori', 
# sottraendo ai valori la lettura prima dell'immersione
x  <- dati$x  - dati$x[1]
dm1 <- dati$dm1 - dati$dm1[1]
dm2 <- dati$dm2 - dati$dm2[1]
dm3 <- dati$dm3 - dati$dm3[1]
dm4 <- dati$dm4 - dati$dm4[1]
n  <- length(x)

print(dati)

pausa()

plot(x[2:n], dm1[2:n], xlab='affondamento (cm)', ty='b', log='xy', 
     ylab="massa d'acqua spostata (g)", ylim=c(0.05,20),
     pch=1, cex=1.5, col='darkblue',
     main="Massa d'acqua spostata in funzione dell'affondamento")
grid()
points(x[2:n], dm2[2:n], pch=1, ty='b', cex=1.5, col='blue')
points(x[2:n], dm3[2:n], pch=17, ty='b', cex=1.5, col='orange')
points(x[2:n], dm4[2:n], 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')