Fixed bug: program no longer dependant on dt size. The same maximum velocity was found with dt = .01 and .001.

Added *rho parameter to Fdrag, which is passed through from velocity to acceleration which allows simulation 2 to be run without using a different function. Just pass rho = 1.