Dependent mass matrix with RadauIIA5

No, don’t invert it (that won’t be possible with DAEs anyways). Just write the version that’s the equivalent equation with constant mass matrices through substitution.