Thrust gradients and quadrotor dircol example

I think inline differentiation (taylorvar) will always return zero until you add sparse(...)*q(1). Someone also used that trick in the TimeSteppingRigidBodyManipulator. It's not something your test would catch because your dforce is also always zero, but I guess it's right for the wrong reasons?