d9r5.txt

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

import java.text.*;
public class d9r5
{
	public static void main (String[] args)
	{
		//program d9r5
		//driver for routine svdvar
		int i,j,mp = 6;
		int ncvm,ma = 3;
		ncvm = ma;
        double v[] = new double[37];
        double w[] = new double[7];
        double cvm[] = new double[10];
        double tru[][] = new double[4][4];
		w[1] = 0.0;
		w[2] = 1.0;
		w[3] = 2.0;
		w[4] = 3.0;
		w[5] = 4.0;
		w[6] = 5.0;
		tru[1][1] = 1.25; tru[1][2] = 2.5; tru[1][3] = 3.75;
		tru[2][1] = 2.5;  tru[2][2] = 5.0; tru[2][3] = 7.5;
		tru[3][1] = 3.75; tru[3][2] = 7.5; tru[3][3] = 11.25;
        d9r5F g = new d9r5F();
        DecimalFormat form = new DecimalFormat("0.00");
        for (j = 1; j <= mp; j++)
		{
            for (i = 1; i <= mp; i++)
			{
                v[(i - 1) * mp + j] = i;
			}
		}
		System.out.println("Matrix v");
		System.out.println();
        for (j = 1; j <= mp; j++)
		{
            for (i = 1; i <= mp; i++)
			{
                System.out.print("  ");
                System.out.print(v[(j - 1) * mp + i] + "  ");
            }
            System.out.println();
		}
		System.out.println();
		System.out.println("Matrix w");
		System.out.println();
        for (i = 1; i <= mp; i++)
		{
            System.out.print("    ");
            System.out.print(w[i] + "  ");
		}
		g.svdvar(v, ma, mp, w, cvm, ncvm);
		System.out.println();
		System.out.println("Covariance matrix from svdvar");
		System.out.println();
        for (j = 1; j <= ma; j++)
		{
            for (i = 1; i <= ma; i++)
			{
				System.out.print("   ");
                System.out.print(form.format(cvm[(j - 1) * ma + i]) + "  ");
            }
            System.out.println();
		}
		System.out.println();
		System.out.println("Expected covariance matrix");
		System.out.println();
        for (j = 1; j <= ma; j++)
		{
            for (i = 1; i<=ma; i++)
			{
                System.out.print("   ");
                System.out.print(form.format(tru[j][i]) + "  ");
            }
            System.out.println();
		}
	}
}