diff --git a/sed/plots.R b/sed/plots.R
new file mode 100644
index 0000000..8562276
--- /dev/null
+++ b/sed/plots.R
@@ -0,0 +1,76 @@
+data = read.table("mehdipour2013.tab")
+plot (data,log="xy",type="l")
+
+data=read.table("magdziarz1997.tab")
+plot (data,log="xy",type="l")
+
+
+
+
+data = read.table("mehdipour2013.tab")
+plot (data,log="xy",type="l",xlim=c(.001,300),ylim=c(10e11,10e13))
+# mehdipour axes
+abline(v=0.001)
+abline(h=40e11)
+
+data = read.table("mehdipour2013.tab")
+plot (data,log="xy",type="l")
+# mehdipour axes
+abline(v=0.001)
+abline(h=40e11)
+
+
+
+data=read.table("magdziarz1997.tab")
+plot (data,log="xy",type="l",xlim=c(.001,300),ylim=c(1e-3,1e-1))
+#magdziarz boundaries
+abline(v=.0004)
+abline(v=.03)
+abline(v=.4)
+abline(v=3)
+abline(v=50)
+abline(v=500)
+
+#low uv peak
+abline(h=0.04)
+abline(v=0.008)
+
+#high uv peak
+abline(h=0.024)
+abline(v=0.012)
+
+#xray peak
+abline(h=0.03)
+abline(v=60)
+
+coords=read.table("magdziarz1997_samples.tab")
+points(coords)
+
+
+
+
+
+hc=9.113e2
+
+inci = read.table("cont_combined_inci",header=TRUE)
+diffuse = read.table("cont_combined_diffuse",header=TRUE)
+
+
+# To plot inci with error bars for bin size
+plot(Enr,nFn,log="x",xlim=c(20,25),type="l")
+arrows(Enr-d.anu.,nFn,Enr+d.anu.,nFn,length=0.05,angle=90,code=3)
+arrows(Enr-d.anu.,nFn+5,Enr+d.anu.,nFn-5,length=0.05,angle=90,code=3)
+arrows(Enr-d.anu.,nFn+5,Enr+d.anu.,nFn+5,length=0.05,angle=90,code=3)
+arrows(Enr-d.anu.,nFn+50,Enr+d.anu.,nFn+50,length=0.05,angle=90,code=3)
+arrows(Enr-d.anu.,nFn+1e9,Enr+d.anu.,nFn+1e9,length=0.05,angle=90,code=3)
+
+
+
+plot(energy.Ryd,Total,log="x",type="l",xlim=c(1e-2,5))
+arrows(energy.Ryd-d.anu.,Total,energy.Ryd+d.anu.,Total,length=0.05,angle=90,code=3)
+plot(energy.Ryd,ConEmitLocal,log="x",type="l",xlim=c(1e-2,5))
+
+lambda=hc/energy.Ryd
+width_lambda=lambda*d.anu./energy.Ryd
+plot(lambda,Total,log="x",type="l",ylim=c(0,0.2),xlim=c(2000,9400))
+arrows(lambda-width_lambda/2,Total,lambda+width_lambda/2,Total,length=0.05,angle=90,code=3)