import numpy
v0 = 2
a = 0.2
dt = 0.1  # Increment
n = int(round(2/dt)) + 1  # No of t values

t_values = numpy.linspace(0, 2, n+1)
s_values = v0*t_values + 0.5*a*t_values**2

# Make nicely formatted table
for t, s in zip(t_values, s_values):
    print('%.2f  %.4f' % (t, s))
