朗之万方程含时间延迟模拟 c++

dp/dt=f(p，pt)+...

这是主方程，没有时间延迟的已经模拟出来了，现在有时间延迟，请大神提供一个思路，谢谢！

可以分享以下没有时间延迟的吗？