rkdumb.cpp

此程序为VC++常用数值算法这本书附赠的光盘中包含了本书中全部的源代码

void rkdumb(double vstart[5], double nvar, double x1, double x2, double nstep, double xx[201], double y[11][201])
{
	int i,k;
    double v[11], dv[11],x,h;
    for( i = 1; i<=nvar; i++)
	{
        v[i] = vstart[i];
        y[i][ 1] = v[i];
    }
    xx[1] = x1;
    x = x1;
    h = (x2 - x1) / float(nstep);
    for( k = 1; k<=nstep; k++)
	{
        derivs(x, v, dv);
        rk4(v, dv, nvar, x, h, v);
        if( x + h == x)
		{
            cout<<" stepsize not signficant in rkdumb."<<endl;
            _c_exit();
        }
        x = x + h;
        xx[k + 1] = x;
        for( i = 1; i<=nvar; i++)
		{
            y[i][k + 1] = v[i];
        }
    }
    for(i=1; i<=10; i++)
	{
		v[i]=0.0;
		dv[i]=0.0;
	}
}