ass.cpp

来自「工程算法 这是一个很有用的工程数值算法集锦」· C++ 代码 · 共 15 行

double ass(int n, double b)
{
	int i;
    double a[16],temp,db;
    a[1] = exp(-b) / b;
    if( n == 1)  goto yi;
    db = 1.0 / b;
    for (i = 2;i<=n;i++)
	{
        a[i] = a[1] + db * (i - 1) * a[i - 1];
	}
yi:  temp = a[n];
	 return temp;
}