Numerical integral solving method for initial value problem of rigid ordinary differential equation in simulation

The invention discloses a numerical integral solving method suitable for an initial value problem of a rigid ordinary differential equation, provides a d-LM method on the basis of ensuring the solving precision, and further improves the implicit operation stability. On the basis of filling the domestic blank, the operation efficiency is close to that of foreign like products; compared with a foreign commercial system simulation software tool, the numerical method is more novel, a new numerical solution method after 1990s is realized and applied, and the numerical solution process is more stable. 本发明公开了一种适用于刚性常微分方程初值问题的数值积分求解方法，在保证求解精度基础上提出了d-LM方法，进一步提高了隐式运算稳定性。在填补国内空白的基础上，其运算效率接近国外同类产品；与国外商用系统仿真软件工具相比，本发明的数值方法更新颖，对1990年代之后的新的数值求解方法进行了实现与应用，数值求解过程更加稳定。