run.cpp

DFP变尺度法

//run.cpp
#include "dfp.h"

double fun(CMatrix& X);

void main()
{
	ofstream o(".\\result.txt");   //用于清除 result.txt中的记录
	CMatrix X0(2,1);
	X0(1,1)=-2,X0(2,1)=2;
	CDfp dfp;
	CMatrix result=dfp(fun,X0,0.01);
	return;
}

 //a=1
double fun(CMatrix& X)
{
	double x1=X(1,1),x2=X(2,1);
	return pow(x1,4)-2*x1*x1*x2-2*x1*x2+x1*x1+2*x2*x2+4.5*x1-4*x2+5;
}


/* //a=0
double fun(CMatrix& X)
{
	double x1=X(1,1),x2=X(2,1);
	return pow(x1-5,2)*4+pow(x2-6,2);
}
*/

/*
double fun(CMatrix& X)
{
	double x1=X(1,1),x2=X(2,1);
	return x1*x1+x2*x2-x1*x2-10*x1-4*x2+60;
}
*/