如何用C++实现时滞方程的龙格－库塔法求解？

xsj3917

现有附件中的二维时滞方程，想用四阶龙格库塔法求解，如何用C++去实现？注意该方程的时滞项是变化的
附：一般微分方程的四阶龙格－库塔法C++代码：
for (i=0;i<n;i++)
{
x[i+1]=x+h;
k1=function(x,y);
k2=function(x+h/2,y+h*k1/2);
k3=function(x+h/2,y+h*k2/2);
k4=function(x+h,y+h*k3);
y[i+1]=y+h*(k1+2*k2+2*k3+k4)/6;
}
其中function为方程，例如

double function(double x,double y)
{
return y-2*x/y;
}

[ 本帖最后由 xsj3917 于 2006-11-9 14:41 编辑 ]

风花雪月

这个一般采用隐式Runge-Kutta法进行求解，具体程序参考
http://forum.vibunion.com/forum/viewthread.php?tid=31492

[ 本帖最后由 风花雪月 于 2007-1-3 09:41 编辑 ]

mas3593

<c语言精彩编程百例>上有,中国水利水电出版社