def ds1(v): d=0.005*v**2+v/3.6 return d