d9r1.java

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

import java.text.*;
public class d9r1
{
	public static void main (String[] args)
	{
		//program d9r1
		//driver for routine fit
		int i, mwt, npt = 100;
		double spread = 0.5;
        double x[] = new double[101];
        double y[] = new double[101];
        double sig[] = new double[101];
		double a, b, siga, sigb, chi2, q;
		int idum = -117;
        d9r1F g = new d9r1F();
		DecimalFormat form = new DecimalFormat("0.000000");
        a = 0; b = 0; siga = 0; sigb = 0; chi2 = 0; q = 0;
		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;
		}
		for (mwt = 0; mwt <= 1; mwt++)
		{
            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;  
            if (mwt == 0)
			{
                System.out.println("Ignoring standard deviation");
			}
            else
			{
                System.out.println("Including standard deviation");
            }
            System.out.print("a = " + form.format(a) + "       ");
            System.out.println("Uncertainty: " + form.format(siga));
            System.out.print("b = "+form.format(b) + "       ");
            System.out.println("Uncertainty: " + form.format(sigb));
            System.out.println("Chi-squared: " + form.format(chi2));
            System.out.println("Goodness-of-fit: " + form.format(q));            
		}
	}
}