In [6]:
#Plotting the data and then fitting it using polyfit()
plt.plot(x,y,'D')

model = np.polyfit(x,y,1)
model
Out[6]:
array([  4.68585196, 633.9931736 ])