work around a case where Eigen seems to produce an incorrect result (and valgrind flags an uninitialized read)
Improve Natarajan regularization support:
1) Return also B02 (coefficient of input time derivative) and E (feedthrough) matrices
2) Implement iterative regularization for cases where C remains singular after the first iteration