ReturnMatrix dtau_dq(const ColumnVector & q, const ColumnVector & qp, const ColumnVector & qpp);
Matrix