ksone.cpp

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

void ksone(double data[], int n, double& d, double& prob)
{
    sort(n, data);
    int en = n;
    d = 0.0;
    double fq,ff,dt,fo = 0.0;
    for (int j = 1; j<=n; j++)
	{
        fq = (double)j /(double)en;
        ff = func(data[j]);
        dt = fabs(fo - ff);
        if (fabs(fq - ff) > dt)
		{
			dt = fabs(fq - ff);
		}
        if (dt > d)
		{
			d = dt;
		}
        fo = fq;
    }
    prob = probks(sqrt(en) * d);
}