Obtaining solutions to large-scale linear ordinary differential equations (ODEs) inevitably relies on numerical integration methods that have developed in various contexts to ultimately approximate the differentiation of continuous-time values using differences in discrete values. This approximation typically induces integration error, which we hypothesize as arising from the deviation of the characteristic roots from those of the original ODEs. Herein, this hypothesis was proved by numerical experiments on integr…