bessy.cpp

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

double bessy(int n, double x)
{
	double tox,by,bym,byp,temp;
	int j;
    if (n < 2)
	{
        cout<<"bad argument n in bessy";
        _c_exit();
	}
    tox = 2.0 / x;
    by = bessy1(x);
    bym = bessy0(x);
    for( j = 1;j<=n-1;j++)
	{
        byp = j * tox * by - bym;
        bym = by;
        by = byp;
    }
    temp = by;
	return temp;
}