#include "MpUtil.h"


int Func1(const Vector<double> &x, double &F, Vector<double> &CF)
{
	F = Sqr(x(0) - 0.5) + Sqr(x(1) - 0.5);
	//CF(0) = -x(1) + One / x(0);
	CF(0) = x(1) - Sqr(x(0));
	return 0;
}

int main()
{
	Print("Begin");

	int n = 2;

	Matrix<double> A(n,0);
	Vector<double> b(0);

	MinNewton m(Func1,A,b,n,0,0,1,0);
	m.x(0) = 0.1;
	m.x(1) = 0.1;
	m.Solve();
	Print("m.x",m.x);
	PauseExecution();

	MinPenalty m2(Func1,n,1,0);
	m2.x(0) = 0.1;
	m2.x(1) = 0.1;
	m2.Solve();
	Print("m2.x",m2.x);
	PauseExecution();

	MinPenalty m3(Func1,n,1,0);
	m3.x(0) = 0.1;
	m3.x(1) = 0.1;
	m3.PenaltyInitial = 10000000000.0;
	m3.Solve();
	Print("m2.x",m2.x);
	PauseExecution();

	Print("End");
	return 0;
}