d9r5f.java

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

public class d9r5F
{
	void svdvar(double v[], int ma, int np, double w[], double cvm[], int ncvm)
	{
        double wti[] = new double[21];
		int i, j;
		for (i = 1; i <= ma; i++)
		{
			wti[i] = 0.0;
			if (w[i] != 0.0)
			{
				wti[i] = 1.0 / (w[i] * w[i]);
			}
		}
		for (i = 1; i <= ma; i++)
		{
			for (j = 1; j <= i; j++)
			{
				double sum1 = 0.0;
				for (int k = 1; k <= ma; k++)
				{
					sum1 = sum1 + v[(i - 1) * np + k] * v[(j - 1) * np + k] * wti[k];
				}
				cvm[(i - 1) * ma + j] = sum1;
				cvm[(j - 1) * ma + i] = sum1;
			}
		}
	}
}