void dTdqi(Matrix & dRot, ColumnVector & dp, const int i); ReturnMatrix dTdqi(const int i);
Matrix or None (in this case dRot and dp are modified on output)