table <- read.table('dades/dades.txt', sep="\t", dec='.', header=TRUE)
table$LogIntensitatReal <- log(table$IntensitatReal)

fit <- lm(LogIntensitatReal~GruixMassic, data = table)
summary(fit)

svg(file='output/p1a.svg', width=6, height=5)
plot(table$GruixMassic, table$LogIntensitatReal, xlab="Gruix màssic (mg/cm^2)", ylab="log(Intensitat · s)")
abline(fit, col = "mediumorchid2")
dev.off()