diff --git a/quad12/labmoderna/p3/r/Pic_Am.r b/quad12/labmoderna/p3/r/Pic_Am.r
new file mode 100644
index 0000000..54c0ef4
--- /dev/null
+++ b/quad12/labmoderna/p3/r/Pic_Am.r
@@ -0,0 +1,98 @@
+library(tikzDevice)
+
+nom <- 'Pic_Am'
+scale <- 4
+
+plotbg <- function(table, number) {
+  col1 = ""
+  col2 = ""
+  if (number == 0) {
+    col1 = "#073B4C"
+    col2 = "#C7EDFA"
+  }
+  if (number == 2) {
+    col1 = "#118AB2"
+    col2 = "#B4E7F8"
+  }
+  if (number == 4) {
+    col1 = "#06DBA2"
+    col2 = "#AFFDE8"
+  }
+  if (number == 6) {
+    col1 = "#FFC233"
+    col2 = "#FFF3D6"
+  }
+  if (number == 8) {
+    col1 = "#EF436B"
+    col2 = "#FAC7D3"
+  }
+  if (number == 10) {
+    col1 = "#8377D1"
+    col2 = "#E2E0F5"
+  }
+
+  if (number == 0) {
+    plot(table$Energy, table$Impulses, type="l", col=col2, xlab="Energia (keV)", ylab="Comptes")
+  } else {
+    lines(table$Energy, table$Impulses, type="l", col=col2, xlab="Energia (keV)", ylab="Comptes")
+  }
+}
+
+plotfg <- function(table, number) {
+  col1 = ""
+  col2 = ""
+  if (number == 0) {
+    col1 = "#073B4C"
+    col2 = "#C7EDFA"
+  }
+  if (number == 2) {
+    col1 = "#118AB2"
+    col2 = "#B4E7F8"
+  }
+  if (number == 4) {
+    col1 = "#06DBA2"
+    col2 = "#AFFDE8"
+  }
+  if (number == 6) {
+    col1 = "#FFC233"
+    col2 = "#FFF3D6"
+  }
+  if (number == 8) {
+    col1 = "#EF436B"
+    col2 = "#FAC7D3"
+  }
+  if (number == 10) {
+    col1 = "#8377D1"
+    col2 = "#E2E0F5"
+  }
+
+  lines(table$Energy, table$ImpulsesSmooth, col=col1, type="l")
+}
+
+pdf(paste0('output/', nom, '.pdf'), width = 6, height = 5)
+numbers <- c(0:5)*2
+for (number in numbers) {
+  table <- read.table(paste0('dades/', nom, '_', number, '.txt'), sep="\t", dec='.', header=TRUE)
+  plotbg(table, number)
+}
+for (number in numbers) {
+  table <- read.table(paste0('dades/', nom, '_', number, '.txt'), sep="\t", dec='.', header=TRUE)
+  plotfg(table, number)
+}
+legendColors = c("#073B4C", "#118AB2", "#06DBA2", "#FFC233", "#EF436B", "#8377D1")
+legend("topright", inset = c(0.015, 0.03), col = legendColors, legend = c("Buit", "2 mangueres", "4 mangueres", "6 mangueres", "8 mangueres", "10 mangueres"), lty=c(1, 1, 1, 1, 1, 1))
+dev.off()
+
+tikz(paste0('output/', nom, '.tex'), width = 6, height = 5)
+numbers <- c(0:5)*2
+for (number in numbers) {
+  table <- read.table(paste0('dades/', nom, '_', number, '.txt'), sep="\t", dec='.', header=TRUE)
+  plotbg(table, number)
+}
+for (number in numbers) {
+  table <- read.table(paste0('dades/', nom, '_', number, '.txt'), sep="\t", dec='.', header=TRUE)
+  plotfg(table, number)
+}
+legendColors = c("#073B4C", "#118AB2", "#06DBA2", "#FFC233", "#EF436B", "#8377D1")
+legend("topright", inset = c(0.015, 0.03), col = legendColors, legend = c("Buit", "2 mangueres", "4 mangueres", "6 mangueres", "8 mangueres", "10 mangueres"), lty=c(1, 1, 1, 1, 1, 1))
+dev.off()
diff --git a/quad12/labmoderna/p3/r/Ra_aire.r b/quad12/labmoderna/p3/r/Ra_aire.r
new file mode 100644
index 0000000..a9217aa
--- /dev/null
+++ b/quad12/labmoderna/p3/r/Ra_aire.r
@@ -0,0 +1,34 @@
+library(tikzDevice)
+
+nom <- 'Ra_aire'
+scale <- 4
+
+table <- read.table(paste0('dades/', nom, '.txt'), sep="\t", dec='.', header=TRUE)
+
+# Plot the data
+plotall <- function() {
+  plot(table$Energy, table$Impulses, type="l", col="#EEEEFF", xlab="Energia (keV)", ylab="Comptes", ylim=c(0, 100))
+  lines(table$Energy, table$ImpulsesSmooth, type="l")
+
+  # Add arrows and text labels to the peaks
+  for (j in 1:nrow(peaks)) {
+    i <- peaks[j, 'indexes']
+    label <- peaks[j, 'labels']
+    arrows(table$Energy[i], table$ImpulsesSmooth[i] + scale*(if (j == 1) 5 else 3), table$Energy[i], table$ImpulsesSmooth[i] + scale*0.1, length = 0.05, code = 1, col = "#514E7E")
+    text(table$Energy[i], table$ImpulsesSmooth[i] + scale*(if (j == 1) 6.6 else 4.6), paste0(table$Energy[i], " keV", "\n", label), col = "#514E7E", cex = 0.75)
+  }
+}
+
+# Els valors dins de c(...) són les files del fitxer de dades.
+# Es resta 2 per trobar els índexs (2 columnes són headers).
+indexes <- c(207, 327, 645) - 2
+labels <- c("\\textsuperscript{210}Po + \\textsuperscript{222}Rn", "\\textsuperscript{218}Po", "\\textsuperscript{214}Po")
+peaks <- data.frame(indexes, labels)
+
+pdf(paste0('output/', nom, '.pdf'), width = 6, height = 5)
+plotall()
+dev.off()
+
+tikz(paste0('output/', nom, '.tex'), width = 6, height = 5)
+plotall()
+dev.off()
diff --git a/quad12/labmoderna/p3/r/calibracio_na.r b/quad12/labmoderna/p3/r/calibracio_na.r
new file mode 100644
index 0000000..e86eda1
--- /dev/null
+++ b/quad12/labmoderna/p3/r/calibracio_na.r
@@ -0,0 +1,29 @@
+library(tikzDevice)
+
+table <- read.table('dades/calibracio_ra.txt', sep="\t", dec='.', header=TRUE)
+
+# Plot the data
+plotall <- function() {
+  plot(table$ChannelNumber, table$Impulses, type="l", col="#EEEEFF", xlab="Nombre de canal", ylab="Comptes")
+  lines(table$ChannelNumber, table$ImpulsesSmooth, type="l")
+
+  # Add arrows and text labels to the peaks
+  for (j in 1:nrow(peaks)) {
+    i <- peaks[j, 'indexes']
+    label <- peaks[j, 'labels']
+    arrows(table$ChannelNumber[i], table$ImpulsesSmooth[i]+3, table$ChannelNumber[i], table$ImpulsesSmooth[i]+0.1, length = 0.05, code = 1, col = "#514E7E")
+    text(table$ChannelNumber[i], table$ImpulsesSmooth[i]+4.6, paste0("Canal ", table$ChannelNumber[i], "\n", label), col = "#514E7E", cex = 0.75)
+  }
+}
+
+indexes <- c(1601, 2359, 3373) + 1
+labels <- c("\\textsuperscript{226}Ra", "\\textsuperscript{218}Po", "\\textsuperscript{214}Po")
+peaks <- data.frame(indexes, labels)
+
+pdf('output/calibracio_ra.pdf', width = 6, height = 5)
+plotall()
+dev.off()
+
+tikz('output/calibracio_ra.tex', width = 6, height = 5)
+plotall()
+dev.off()
diff --git a/quad12/labmoderna/p3/r/na_calibrat.r b/quad12/labmoderna/p3/r/na_calibrat.r
new file mode 100644
index 0000000..86cb4fe
--- /dev/null
+++ b/quad12/labmoderna/p3/r/na_calibrat.r
@@ -0,0 +1,34 @@
+library(tikzDevice)
+
+nom <- 'Calibrat_Ra'
+scale <- 4
+
+table <- read.table(paste0('dades/', nom, '.txt'), sep="\t", dec='.', header=TRUE)
+
+# Plot the data
+plotall <- function() {
+  plot(table$Energy, table$Impulses, type="l", col="#EEEEFF", xlab="Energia (keV)", ylab="Comptes", ylim=c(0, 100))
+  lines(table$Energy, table$ImpulsesSmooth, type="l")
+
+  # Add arrows and text labels to the peaks
+  for (j in 1:nrow(peaks)) {
+    i <- peaks[j, 'indexes']
+    label <- peaks[j, 'labels']
+    arrows(table$Energy[i], table$ImpulsesSmooth[i] + scale*3, table$Energy[i], table$ImpulsesSmooth[i] + scale*0.1, length = 0.05, code = 1, col = "#514E7E")
+    text(table$Energy[i], table$ImpulsesSmooth[i] + scale*4.6, paste0(table$Energy[i], " keV", "\n", label), col = "#514E7E", cex = 0.75)
+  }
+}
+
+# Els valors dins de c(...) són les files del fitxer de dades.
+# Es resta 2 per trobar els índexs (2 columnes són headers).
+indexes <- c(396, 492, 581, 830) - 2
+labels <- c("\\textsuperscript{226}Ra", "\\textsuperscript{210}Po + \\textsuperscript{222}Rn", "\\textsuperscript{218}Po", "\\textsuperscript{214}Po")
+peaks <- data.frame(indexes, labels)
+
+pdf(paste0('output/', nom, '.pdf'), width = 6, height = 5)
+plotall()
+dev.off()
+
+tikz(paste0('output/', nom, '.tex'), width = 6, height = 5)
+plotall()
+dev.off()