We introduce a model order reduction procedure for differential-algebraic equations, which is based on the intrinsic differential equation contained in the starting system and on the remaining algebraic constraints. We implement numerically this procedure and show numerical evidence of its validity
