多维进退-黄金分割法

#include 
#include 
#include 
#define N 2
double objf（double x[]）
{
	double ff;
	ff=60-10*x[0]-4*x[1]+x[0]*x[0]+x[1]*x[1]-x[0]*x[1];
	return（ff）;
}
void jtf（double x0[]double h0double s[]int ndouble a[]double b[]）
{
	int i;
	double *x[3]hf1f2f3;
	for（i=0;i<3;i++）
	   x[i]=（double *）malloc（n*sizeof（double））;
	h=h0;
	for（i=0;i	   *（x[0]+i）=x0[i];
	f1=objf（x[0]）;
	for（i=0;i	   *（x[1]+i）=*（x[0]+i）+h*s[i];
	f2=objf（x[1]）;
	if（f2>=f1）
	{
		h=-h;    
		for（i=0;i		   *（x[2]+i）=*（x[0]+i）;   
		f3=f1;    
		for（i=0;i		{
			*（x[0]+i）=*（x[1]+i）;     
			*（x[1]+i）=*（x[2]+i）;    
		}   
		f1=f2;f2=f3;   
	}   
	for（;;）   
	{
		h=2*h;     
		for（i=0;i		   *（x[2]+i）=*（x[1]+i）+h*