#include?
#include?
int?const?n=2;?
double?fun（double?x[n]）;?
void?compare（double?x[n+1][n]double?LGH[3][n]double?LGH_tap[3]）;//在調用之前LGH先給x[n+1]賦值?
double?H（double?x[n+1][n]double?LGH[3][n]double?c）;//判別準則?
void?zx（double?x[n+1][n]double?LGH[3][n]double?LGH_tap[3]double?xn_1[n]double?xn_2[n]）;?
void?main（）?
{?
	double?x[n+1][n]={0};?
	double?c=0.001;?
	double?a=2;?
	double?b=0.75;?
	double?LGH[3][n]={0};?
	double?LGH_tap[3]={0};?
	double?xn_1[n]={0};?
	double?xn_2[n]={0};?
	double?xn_3[n]={0};?
	double?xn_4[n]={0};?
	int?ijtap=0;?
	//給出一個初始點X0X1X2?
	x[0][0]=8;x[0][1]=9;?
	x[1][0]=10;x[1][1]=9;?
	x[2][0]=8;x[2][1]=11;?
	cout<<“初始點：“<	for（i=0;i	{?
		cout<<“x[“<		for（j=0;j		{?
			cout<			if（j!=n-1）?
				cout<<““;?
		}?
		cout<<“）“<	}?
	compare（xLGHLGH_tap）;//返回XlXhXg及它們在x[n+1][n]中的下標?
	zx（xLGHLGH_tapxn_1xn_2）;//返回中心X（n+1）和反射點X（n+2）?
	while（H（xLGHc）==0）?
	{?
		if（fun（xn_2）		{?
			for（i=0;i			{?
				xn_3[i]=xn_1[i]+a*（xn_2[i]-xn_1[i]）;//擴張?
			}?
			if（fun（xn_3）			{?
				for（i=0;i				{?
					x[int（LGH_tap[2]）][i]=LGH[2][i]=xn_3[i];?
				}?
				//goto?step?2;?
			}?
			else?
			{?
				for（i=0;i				{?
					x[int（LGH_tap[2]）][i]=LGH[2][i]=xn_2[i];?
				}?
				//goto?step?2;?
			}?
		}?
		else??
		{?
			if（fun（xn_2）			{?
				for（i=0;i				{?
					xn_4[i]=xn_1[i]+b*（xn_2[i]-xn_1[i]）;//壓縮?
				}?
				//goto?step?2;?
			}?
			else??
			{?
				for（i=0;i				{?
					xn_4[i]=xn_1[i]-b*（xn_2[i]-xn_1[i]）;//壓縮?
				}?
			}?
			if（fun（xn_4）			{?
				for（i=0;i				{				?
					x[int（LGH_tap[2]）][i]=LGH[2][i]=xn_4[i];?
				}?
			}?
			else//縮邊?
			{?
				cout<<“進行縮邊。“<				for（i=0;i				{?
					if（i==LGH_tap[0]）?
						continue;?
					for（j=0;j					{?
						x[i][j]=（LGH[0][j]+x[i][j]）/2;?
					}?
				}?
			}		?
		}?
		compare（xLGH