'Return the jacobian matrix of a transformation.'

def jacobian(self, to_sys, coords):

    with_dummies = self.coord_tuple_transform_to(to_sys, self._dummies).jacobian(self._dummies)
    return with_dummies.subs(list(zip(self._dummies, coords)))
