Fixed a few things, added some comments. Now the density calculation is
correct and the terminal velocity consistent with analytically found
values.