diff --git a/xrd/data/.RData b/xrd/data/.RData
new file mode 100644
index 0000000..de34411
Binary files /dev/null and b/xrd/data/.RData differ
diff --git a/xrd/data/.Rhistory b/xrd/data/.Rhistory
new file mode 100644
index 0000000..6b6f8ac
--- /dev/null
+++ b/xrd/data/.Rhistory
@@ -0,0 +1,314 @@
+plastic_count = read.table("plastic.tab")
+length(plastic_count)
+length(plastic_count[,1])
+x = seq(5.00547782
+?Seq
+?seq
+?seq
+x = seq(from = 5.00547782,to = 119.99571965,length.out = 13763)
+x
+length(x)
+length(plastic_count[,1])
+length(x)
+grease_count = table.read("grease.tab")
+grease_count = read.table("grease.tab")
+x = seq(from = 5.00547782,to = 119.99571965,length.out = 13763)
+grease_count = read.table("grease.tab")
+plastic_count = read.table("plastic.tab")
+wood_count = read.table("wood.tab")
+plot(x=x,y=plastic_count[,1],pch='-',ylab=expression(count),xlab=expression(2~theta ~ degree))
+abline(v=42.5)
+axis(3,at=42.5)
+abline(v=14)
+axis(3,at=14)
+abline(v=30.5)
+axis(3,at=30.5)
+dev.print(pdf,"../report/plastic_peaks.pdf",width=4,height=2)
+plot(x=x,y=grease_count[,1],pch='-',ylab=expression(count),xlab=expression(2~theta ~ degree))
+abline(v=12)
+axis(3,at=12)
+abline(v=22)
+axis(3,at=22)
+dev.print(pdf,"../report/grease_peaks.pdf",width=4,height=2)
+plot(x=x,y=wood_count[,1],pch='-',ylab=expression(count),xlab=expression(2~theta ~ degree))
+abline(v=22.5)
+axis(3,at=22.5)
+abline(v=16)
+axis(3,at=16)
+abline(v=34.5)
+axis(3,at=34.5)
+abline(v=46)
+axis(3,at=46)
+dev.print(pdf,"../report/wood_peaks.pdf",width=4,height=2)
+cex.axis=2
+plot(x=x,y=plastic_count[,1],pch='-',ylab=expression(count),xlab=expression(2~theta ~ degree))
+abline(v=42.5)
+axis(3,at=42.5)
+abline(v=14)
+axis(3,at=14)
+abline(v=30.5)
+axis(3,at=30.5)
+dev.print(pdf,"../report/plastic_peaks.pdf",width=4,height=2)
+plot(x=x,y=grease_count[,1],pch='-',ylab=expression(count),xlab=expression(2~theta ~ degree))
+abline(v=12)
+axis(3,at=12)
+abline(v=22)
+axis(3,at=22)
+dev.print(pdf,"../report/grease_peaks.pdf",width=4,height=2)
+plot(x=x,y=wood_count[,1],pch='-',ylab=expression(count),xlab=expression(2~theta ~ degree))
+abline(v=22.5)
+axis(3,at=22.5)
+abline(v=16)
+axis(3,at=16)
+abline(v=34.5)
+axis(3,at=34.5)
+abline(v=46)
+axis(3,at=46)
+dev.print(pdf,"../report/wood_peaks.pdf",width=4,height=2)
+plot(x=x,y=plastic_count[,1],pch='-',ylab=expression(count),xlab=expression(2~theta ~ degree))
+abline(v=42.5)
+axis(3,at=42.5)
+abline(v=14)
+axis(3,at=14)
+abline(v=30.5)
+axis(3,at=30.5)
+dev.print(pdf,"../report/plastic_peaks.pdf",width=4,height=4)
+plot(x=x,y=grease_count[,1],pch='-',ylab=expression(count),xlab=expression(2~theta ~ degree))
+abline(v=12)
+axis(3,at=12)
+abline(v=22)
+axis(3,at=22)
+dev.print(pdf,"../report/grease_peaks.pdf",width=4,height=4)
+plot(x=x,y=wood_count[,1],pch='-',ylab=expression(count),xlab=expression(2~theta ~ degree))
+abline(v=22.5)
+axis(3,at=22.5)
+abline(v=16)
+axis(3,at=16)
+abline(v=34.5)
+axis(3,at=34.5)
+abline(v=46)
+axis(3,at=46)
+dev.print(pdf,"../report/wood_peaks.pdf",width=4,height=4)
+plot(x=x,y=plastic_count[,1],pch='-',ylab=expression(count),xlab=expression(2~theta ~ degree))
+abline(v=42.5)
+axis(3,at=42.5)
+abline(v=14)
+axis(3,at=14)
+abline(v=30.5)
+axis(3,at=30.5)
+dev.print(pdf,"../report/plastic_peaks.pdf",width=2,height=2)
+plot(x=x,y=grease_count[,1],pch='-',ylab=expression(count),xlab=expression(2~theta ~ degree))
+abline(v=12)
+axis(3,at=12)
+abline(v=22)
+axis(3,at=22)
+dev.print(pdf,"../report/grease_peaks.pdf",width=2,height=2)
+plot(x=x,y=wood_count[,1],pch='-',ylab=expression(count),xlab=expression(2~theta ~ degree))
+abline(v=22.5)
+axis(3,at=22.5)
+abline(v=16)
+axis(3,at=16)
+abline(v=34.5)
+axis(3,at=34.5)
+abline(v=46)
+axis(3,at=46)
+dev.print(pdf,"../report/wood_peaks.pdf",width=2,height=2)
+width
+width
+plot(x=x,y=plastic_count[,1],pch='-',ylab=expression(count),xlab=expression(2~theta ~ degree))
+abline(v=42.5)
+axis(3,at=42.5)
+abline(v=14)
+axis(3,at=14)
+abline(v=30.5)
+axis(3,at=30.5)
+dev.print(pdf,"../report/plastic_peaks.pdf",width=3,height=3)
+plot(x=x,y=grease_count[,1],pch='-',ylab=expression(count),xlab=expression(2~theta ~ degree))
+abline(v=12)
+axis(3,at=12)
+abline(v=22)
+axis(3,at=22)
+dev.print(pdf,"../report/grease_peaks.pdf",width=3,height=3)
+plot(x=x,y=wood_count[,1],pch='-',ylab=expression(count),xlab=expression(2~theta ~ degree))
+abline(v=22.5)
+axis(3,at=22.5)
+abline(v=16)
+axis(3,at=16)
+abline(v=34.5)
+axis(3,at=34.5)
+abline(v=46)
+axis(3,at=46)
+dev.print(pdf,"../report/wood_peaks.pdf",width=3,height=3)
+plot(x=x,y=plastic_count[,1],pch='-',log="x",ylab=expression(count),xlab=expression(2~theta ~ degree))
+abline(v=42.5)
+axis(3,at=42.5)
+abline(v=14)
+axis(3,at=14)
+abline(v=30.5)
+axis(3,at=30.5)
+dev.print(pdf,"../report/plastic_peaks.pdf",width=3,height=3)
+plot(x=x,y=grease_count[,1],pch='-',log="x",ylab=expression(count),xlab=expression(2~theta ~ degree))
+abline(v=12)
+axis(3,at=12)
+abline(v=22)
+axis(3,at=22)
+dev.print(pdf,"../report/grease_peaks.pdf",width=3,height=3)
+plot(x=x,y=wood_count[,1],pch='-',log="x",ylab=expression(count),xlab=expression(2~theta ~ degree))
+abline(v=22.5)
+axis(3,at=22.5)
+abline(v=16)
+axis(3,at=16)
+abline(v=34.5)
+axis(3,at=34.5)
+abline(v=46)
+axis(3,at=46)
+dev.print(pdf,"../report/wood_peaks.pdf",width=3,height=3)
+plot(x=x,y=plastic_count[,1],pch='-',log="x",ylab=expression(count),xlab=expression("log(2"~theta~degree~")"))
+abline(v=42.5)
+axis(3,at=42.5)
+abline(v=14)
+axis(3,at=14)
+abline(v=30.5)
+axis(3,at=30.5)
+dev.print(pdf,"../report/plastic_peaks.pdf",width=3,height=3)
+plot(x=x,y=grease_count[,1],pch='-',log="x",ylab=expression(count),xlab=expression("log(2"~theta~degree~")"))
+abline(v=12)
+axis(3,at=12)
+abline(v=22)
+axis(3,at=22)
+dev.print(pdf,"../report/grease_peaks.pdf",width=3,height=3)
+plot(x=x,y=wood_count[,1],pch='-',log="x",ylab=expression(count),xlab=expression("log(2"~theta~degree~")"))
+abline(v=22.5)
+axis(3,at=22.5)
+abline(v=16)
+axis(3,at=16)
+abline(v=34.5)
+axis(3,at=34.5)
+abline(v=46)
+axis(3,at=46)
+dev.print(pdf,"../report/wood_peaks.pdf",width=3,height=3)
+?axis
+?axis
+plot(x=x,y=plastic_count[,1],pch='-',log="x",ylab=expression(count),xlab=expression("log(2"~theta~degree~")"))
+abline(v=42.5)
+axis(3,at=42.5,hadj=1)
+abline(v=14)
+axis(3,at=14)
+abline(v=30.5)
+axis(3,at=30.5)
+dev.print(pdf,"../report/plastic_peaks.pdf",width=3,height=3)
+plot(x=x,y=plastic_count[,1],pch='-',log="x",ylab=expression(count),xlab=expression("log(2"~theta~degree~")"))
+abline(v=42.5)
+axis(3,at=42.5,vadj=1)
+abline(v=14)
+axis(3,at=14)
+abline(v=30.5)
+axis(3,at=30.5)
+dev.print(pdf,"../report/plastic_peaks.pdf",width=3,height=3)
+?axis
+plot(x=x,y=plastic_count[,1],pch='-',log="x",ylab=expression(count),xlab=expression("log(2"~theta~degree~")"))
+abline(v=42.5)
+axis(3,at=42.5,padj=1)
+abline(v=14)
+axis(3,at=14)
+abline(v=30.5)
+axis(3,at=30.5)
+dev.print(pdf,"../report/plastic_peaks.pdf",width=3,height=3)
+plot(x=x,y=plastic_count[,1],pch='-',log="x",ylab=expression(count),xlab=expression("log(2"~theta~degree~")"))
+abline(v=42.5)
+axis(3,at=42.5,padj=-2.4)
+abline(v=30.5)
+axis(3,at=30.5,padj=-1.2)
+abline(v=14)
+axis(3,at=14)
+dev.print(pdf,"../report/plastic_peaks.pdf",width=3,height=3)
+cex.axis=1
+plot(x=x,y=plastic_count[,1],pch='-',log="x",ylab=expression(count),xlab=expression("log(2"~theta~degree~")"))
+abline(v=42.5)
+axis(3,at=42.5,padj=-2.4)
+abline(v=30.5)
+axis(3,at=30.5,padj=-1.2)
+abline(v=14)
+axis(3,at=14)
+dev.print(pdf,"../report/plastic_peaks.pdf",width=3,height=3)
+plot(x=x,y=plastic_count[,1],pch='-',log="x",ylab=expression(count),xlab=expression("log(2"~theta~degree~")"))
+abline(v=42.5)
+axis(3,at=42.5,padj=0)
+abline(v=30.5)
+axis(3,at=30.5,padj=-1.2)
+abline(v=14)
+axis(3,at=14,padj=-2)
+dev.print(pdf,"../report/plastic_peaks.pdf",width=3,height=3)
+plot(x=x,y=wood_count[,1],pch='-',log="x",ylab=expression(count),xlab=expression("log(2"~theta~degree~")"))
+abline(v=16)
+axis(3,at=16,padj=-0.2)
+abline(v=22.5)
+axis(3,at=22.5,padj=-2.5)
+abline(v=34.5)
+axis(3,at=34.5,padj=-1.2)
+abline(v=46)
+axis(3,at=46)
+dev.print(pdf,"../report/wood_peaks.pdf",width=3,height=3)
+plot(x=x,y=wood_count[,1],pch='-',log="x",ylab=expression(count),xlab=expression("log(2"~theta~degree~")"),xlim=c(15,55)
+abline(v=16)
+axis(3,at=16,padj=-0.2)
+abline(v=22.5)
+axis(3,at=22.5,padj=-2.5)
+abline(v=34.5)
+axis(3,at=34.5,padj=-1.2)
+abline(v=46)
+axis(3,at=46)
+dev.print(pdf,"../report/wood_peaks.pdf",width=3,height=3)
+plot(x=x,y=wood_count[,1],pch='-',log="x",ylab=expression(count),xlab=expression("log(2"~theta~degree~")"),xlim=c(15,55))
+abline(v=16)
+axis(3,at=16,padj=-0.2)
+abline(v=22.5)
+axis(3,at=22.5,padj=-2.5)
+abline(v=34.5)
+axis(3,at=34.5,padj=-1.2)
+abline(v=46)
+axis(3,at=46)
+dev.print(pdf,"../report/wood_peaks.pdf",width=3,height=3)
+plot(x=x,y=wood_count[,1],pch='-',log="x",ylab=expression(count),xlab=expression("log(2"~theta~degree~")"),xlim=c(15,55))
+abline(v=16)
+axis(3,at=16)
+abline(v=22.5)
+axis(3,at=22.5)
+abline(v=34.5)
+axis(3,at=34.5)
+abline(v=46)
+axis(3,at=46)
+dev.print(pdf,"../report/wood_peaks.pdf",width=3,height=3)
+plot(x=x,y=plastic_count[,1],pch='-',log="x",ylab=expression(count),xlab=expression("log(2"~theta~degree~")"),xlim=c(12,55))
+abline(v=42.5)
+axis(3,at=42.5)
+abline(v=30.5)
+axis(3,at=30.5)
+abline(v=14)
+axis(3,at=14)
+dev.print(pdf,"../report/plastic_peaks.pdf",width=3,height=3)
+plot(x=x,y=plastic_count[,1],pch='-',log="x",ylab=expression(count),xlab=expression(2~theta~degree),xlim=c(12,55))
+abline(v=42.5)
+axis(3,at=42.5)
+abline(v=30.5)
+axis(3,at=30.5)
+abline(v=14)
+axis(3,at=14)
+dev.print(pdf,"../report/plastic_peaks.pdf",width=3,height=3)
+plot(x=x,y=grease_count[,1],pch='-',log="x",ylab=expression(count),xlab=expression(2~theta~degree))
+abline(v=12)
+axis(3,at=12)
+abline(v=22)
+axis(3,at=22)
+dev.print(pdf,"../report/grease_peaks.pdf",width=3,height=3)
+plot(x=x,y=wood_count[,1],pch='-',log="x",ylab=expression(count),xlab=expression(2~theta~degree),xlim=c(15,55))
+abline(v=16)
+axis(3,at=16)
+abline(v=22.5)
+axis(3,at=22.5)
+abline(v=34.5)
+axis(3,at=34.5)
+abline(v=46)
+axis(3,at=46)
+dev.print(pdf,"../report/wood_peaks.pdf",width=3,height=3)
+quit()