Calculate a as a function of v, from this the new v and the new s in the time interval DeltaT. Now first set c=0 (no air resistance). Then create a copy of the simulation and change the start value to c>0.