+#! /usr/bin/env python
+
+import matplotlib
+#matplotlib.use('GTK')
+matplotlib.use('Agg')
+import matplotlib.pyplot as plt
+import matplotlib.cm as cm
+import numpy as np
+import sys
+
+x,y= np.loadtxt('saxs.data',usecols=(0,1),unpack=True)
+x1,y1= np.loadtxt('MODEL1_saxs.data',usecols=(0,2),unpack=True)
+x2,y2= np.loadtxt('MODEL2_saxs.data',usecols=(0,2),unpack=True)
+x3,y3= np.loadtxt('MODEL3_saxs.data',usecols=(0,2),unpack=True)
+x4,y4= np.loadtxt('MODEL4_saxs.data',usecols=(0,2),unpack=True)
+x5,y5= np.loadtxt('MODEL5_saxs.data',usecols=(0,2),unpack=True)
+
+
+plt.clf()
+plt.xlabel('r')
+plt.ylabel('P(r)')
+plt.xlim(x[0], x[-1])
+plt.plot(x,y,'-',color='black')
+plt.plot(x1,y1,'-',color='red')
+plt.plot(x2,y2,'-',color='lime')
+plt.plot(x3,y3,'-',color='blue')
+plt.plot(x4,y4,'-',color='magenta')
+plt.plot(x5,y5,'-',color='cyan')
+plt.legend(['reference','MODEL1','MODEL2','MODEL3','MODEL4','MODEL5'])
+plt.savefig('saxs.png')