电磁暂态程序求助，编程中遇到了个新问题

lvyinchenxi

近来，想完成这样一个事情：用自己编的程序实现故障前后即稳态和暂态的过程，程序是用C语言编写的。现在有这么个问题：在稳态和突然发生故障的交界处，网络结构发生变化，即导纳矩阵发生变化，造成了以后每一步运算都与EMTP的仿真有个误差（500KV系统误差在10KV左右），使得总的仿真结果不正确。模型中线路采用贝瑞容模型、双端三相电源、系统阻抗。我的qq404713548，有对这个问题感兴趣的或者以前有做个这个方面工作的，请帮个忙。