#include
#include“math.h“
using?namespace?std;
#define	MIU???0.1
#define?SIGMA?0.5
void?wolfe（）
{
	double?x1[2]={00};
	double?x2[2]={00};
	int????pk[2]={10};
	double?a=0b=10000000000alpha=1j=0;
	x2[0]=x1[0]+alpha*pk[0];
	x2[1]=x1[1]+alpha*pk[1];
	double?fg2[2]g1[2];
	f=100*（x2[1]-x2[0]*x2[0]）*（x2[1]-x2[0]*x2[0]）+（1-x2[0]）*（1-x2[0]）;
	g1[0]=-400*（x1[1]-x1[0]*x1[0]）*x1[0]-2*（1-x1[0]）;
	g1[1]=200*（x1[1]-x1[0]*x1[0]）;
	g2[0]=-400*（x2[1]-x2[0]*x2[0]）*x2[0]-2*（1-x2[0]）;
	g2[1]=200*（x2[1]-x2[0]*x2[0]）;
	double?g2pk=g2[0]*pk[0]+g2[1]*pk[1];
	double?g1pk=g1[0]*pk[0]+g1[1]*pk[1];
	bool?conditionA=100*（x1[1]-x1[0]*x1[0]）*（x1[1]-x1[0]*x1[0]）+（1-x1[0]）*（1-x1[0]）-f>=-MIU*alpha*g2pk?true:false;
	bool?conditionB=g2pk>=SIGMA*g1pk?true:false;
???cout<<‘j‘<<‘\t‘<<“alpha“<<‘\t‘<<“Xk+1“<<‘t‘<<‘f‘<<‘\n‘;
	while（!conditionA||!conditionB）
	{
		if（!conditionA）
		{
			j=j+1;
			b=alpha;
			alpha=（a+alpha）/2;
			x1[0]=0x1[1]=0;
????????????x2[0]=x1[0]+alpha*pk[0];
	????????x2[1]=x1[1]+alpha*pk[1];
			f=100*（x2[1]-x2[0]*x2[0]）*（x2[1]-x2[0]*x2[0]）+（1-x2[0]）*（1-x2[0]）;
			g1[0]=-400*（x1[1]-x1[0]*x1[0]）*x1[0]-2*（1-x1[0]）;
	????????g1[1]=20