from numpy import arange, exp, inf
g=9.8##m/s**2# acceleration due to gravity
m=68.1##kg
c=12.5##kg/sec# drag coefficient
v=[]
j=0
for i in arange(0,12,2):
v.append(g*m*(1-exp(-c*i/m))/c)
print "Time(s) = %d\t"%i,"v(m/s) = %0.2f"%v[j]
j+=1
print "Time(s) = %0.2f\t"%inf,"v(m/s) = %0.2f"%(g*m/c)
from numpy import arange, exp, inf
g=9.8##m/s**2# acceleration due to gravity
m=68.1##kg
c=12.5##kg/sec# drag coefficient
count=0#
v=[]
v.append(0)
print "Time(s) = %d\t"%(0),"v(m/s) = %0.2f"%v[0]
for i in arange(1,13,2):
v.append(v[(count-1)]+(g-c*v[(count)]/m)*(2))
print "Time(s) = %d\t"%(i+1),"v(m/s) = %0.2f"%v[(count+1)]
count=count+1#
print "Time(s) = %0.2f\t"%inf,"v(m/s) = %0.2f\t"%(g*m/c)