d9r3f.java

解具有导数的多元函数的无约束极小值点或具有导数的非线性方程组的近似解

public class d9r3F
{	
	void covsrt(double covar[], int ncvm, int ma, int lista[], int mfit)
	{
		int i, j;
		double swap;
		for (j = 1; j <= ma - 1; j++)
		{
			for (i = j + 1; i <= ma; i++)
			{
                covar[(i - 1) * ma + j] = 0.0;
			}
        }
        for (i = 1; i <= mfit - 1; i++)
		{
            for (j = i + 1; j <= mfit; j++)
			{
                if (lista[j] > lista[i])
				{
                    covar[(lista[j] - 1) * ma + lista[i]] = covar[(i - 1) * ma + j];
				}
                else
				{
                    covar[(lista[i] - 1) * ma + lista[j]] = covar[(i - 1) * ma + j];
                }
            }
        }
        swap = covar[1];;
        for (j = 1; j <= ma; j++)
		{
            covar[j] = covar[(j - 1) * ma + j];
            covar[(j - 1) * ma + j] = 0.0;
        }
        covar[(lista[1] - 1) * ma + lista[1]] = swap;
        for (j = 2; j <= mfit; j++)
		{
            covar[(lista[j] - 1) * ma + lista[j]] = covar[j];
        }
        for (j = 2; j <= ma; j++)
		{
            for (i = 1; i <= j - 1; i++)
			{
                covar[(i - 1) * ma + j] = covar[(j - 1) * ma + i];
			}
		}
	}
}