'Augments qdots to the forcing vector above.'

@property
def forcing_full(self):

    if (self.eom is None):
        raise ValueError('Need    to    compute    the    equations    of    motion    first')
    if self.coneqs:
        return self._qdots.col_join(self.forcing).col_join(self._f_cd)
    else:
        return self._qdots.col_join(self.forcing)
