d9r10.java

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

import java.text.*;
public class d9r10
{
	public static void main (String[] args)
	{
		//program d9r10
		//driver for routine medfit
        int i, npt = 100;
        double a, b, siga, sigb, chi2, q, abdev, spread = 0.1;
        a = 0.0; b = 0.0; siga = 0.0; sigb = 0.0; chi2 = 0.0; q = 0.0; abdev = 0.0;
        double x[] = new double[101];
        double y[] = new double[101];
        double sig[] = new double[101];
		int idum = -1984;
        d9r10F g = new d9r10F();
        DecimalFormat form = new DecimalFormat("0.0000");
        for (i = 1; i <= npt; i++)
		{
            x[i] = 0.1 * i;
            y[i] = -2.0 * x[i] + 1.0 + spread * g.gasdev(idum);
            idum = g.gasdev_idum;
            sig[i] = spread;
		}
		int mwt = 1;
		g.fit(x, y, npt, sig, mwt, a, b, siga, sigb, chi2, q);
        a = g.fit_a;
        b = g.fit_b;
        siga = g.fit_siga;
        sigb = g.fit_sigb;
        chi2 = g.fit_chi2;
        q = g.fit_q;
		System.out.println("According to routine fit the result is:");
		System.out.println();
		System.out.print("a = ");
        System.out.print(form.format(a) + "    ");
		System.out.print("Uncertainty: ");
        System.out.println(form.format(siga) + "    ");
		System.out.println();
		System.out.print("b = ");
        System.out.print(form.format(b) + "    ");
		System.out.print("uncertainty:");
        System.out.println(form.format(sigb) + "  ");
		System.out.println();
        System.out.print("chi-squared: " + "   ");
        System.out.print(form.format(chi2) + " for   ");
        System.out.println(form.format(npt) + "points");
		System.out.println();
		System.out.print("Goodness-of-fit: ");
        System.out.println(form.format(q) + "    ");
		System.out.println();
		System.out.println("According to routine medfit the result is:");
		System.out.println();
		g.medfit(x, y, npt, a, b, abdev);
        a = g.medift_a;
        b = g.medift_b;
        abdev = g.medift_abdev;
		System.out.print("a = ");
        System.out.print(form.format(a) + "    ");
		System.out.println();
		System.out.print("b = ");
        System.out.println(form.format(b) + "    ");
		System.out.println();
		System.out.print("Absolute deviation (per data point): ");
        System.out.print(form.format(abdev) + "   ");
		System.out.println();
		System.out.print("(note: Gaussian spread is  ");
        System.out.print(form.format(spread));
		System.out.println(")");
		System.out.println();
		System.out.println();
	}
}