d9r9f.java

来自「解具有导数的多元函数的无约束极小值点或具有导数的非线性方程组的近似解」· Java 代码 · 共 20 行

public class d9r9F
{
	double fgauss_y; 
	void fgauss(double x, double a[], double y, double dyda[], int na)
	{
		y = 0.0;
        for (int i = 1; i <= na - 1; i = i + 3)
		{
			double arg = (x - a[i + 1]) / a[i + 2];
			double ex = Math.exp(-(arg * arg));
			double fac = a[i] * ex * 2.0 * arg;
			y = y + a[i] * ex;
			dyda[i] = ex;
			dyda[i + 1] = fac / a[i + 2];
			dyda[i + 2] = fac * arg / a[i + 2];
		}
        fgauss_y = y;
	}
}