On bias of transfer function of numerical ordinary differential equation solvers in the frequency domain

The bias of transfer function in numerical solution of ordinary differential equations is studied in frequency domain. Explicit formulas are given for two particular cases widely used in reactor physics. Comparing the frequency characteristics of numerical method with the real transfer function one can reduce the inaccuracy with a suitable selection of step-size