d11r7.cpp

Visual C++ 常用数值算法集 源代码

#include "iostream.h"
#include "math.h"

double func2(double x[], int n)
{
    double f = 0.0;
    f = f + (x[1] - 1) * (x[1] - 1);
    f = f + (x[2] - 1) * (x[2] - 1);
    f = f + (x[3] - 1) * (x[3] - 1);
    return f;
}

double func(double x)
{
    return f1dim(x);
}

void main()
{
    //program d11r7
    //driver for routine linmin
    int ndim = 3;
    double x,sr2,fret,pio2 = 1.5707963;
    double p[4], xi[4];
    cout<<endl;
    cout<<"minimum of a 3-d quadratic centered";
    cout<<endl;
    cout<<"at (1.0,1.0,1.0). minimum is found";
    cout<<endl;
    cout<<"along a series of radials."<<endl;
    cout<<"  x         y         z        minimum"<<endl;
    int i,j;
	cout.setf(ios::fixed|ios::left);
	cout.precision(6);
	for (i = 0; i<=10; i++)
	{
        x = pio2 * i / 10.0;
        sr2 = sqrt(2.0);
        xi[1] = sr2 * cos(x);
        xi[2] = sr2 * sin(x);
        xi[3] = 1.0;
        p[1] = 0.0;
        p[2] = 0.0;
        p[3] = 0.0;
        linmin(p, xi, ndim, fret);
        for (j = 1; j<=3; j++)
		{
			cout<<"  ";
			cout.width(11);
            cout<<p[j];
        }
		cout<<"  ";
		cout.width(11);
        cout<<fret;
		cout<<endl;
    }
}