data fit hw

import numpy
import pylab

data = numpy.loadtxt('data_fitting.dat')
t = data[:,0]
v = data[:,1]
s = data[:,2]

poly = numpy.polyfit(t,v,3)
vals = numpy.polyval(poly,t)
pylab.figure()
pylab.title('Velocity vs. Time (with Error)')
pylab.errorbar(t, vals, yerr=s, fmt='-')
pylab.show()