请教时滞微分方程的4阶龙格库塔算法求解

sch_LNQ

对于延迟微分方程用dde23可以求解，但是我觉得控制起来不太方便，所以想利用4阶龙格库塔算法写个，但是我不不知道怎么处理延迟项，比如在写龙格库塔算法时，变量会叠加0.5*h*k_x（其中h为步长，k_x为变量x的叠加量），那么x(t-tau)是否必须也叠加这个量呢（tau为延迟时间）。
希望得到指导。
真切希望得到无水大哥和其他高手的帮忙。大家给我个求解程序也行，谢谢了

无水1324

呵呵，没有做过这样子的。不过你可以尝试一下，把方程离散化，然后求解。