坐标轮换法C语言相关程序

#include<stdio.h>
#include<math.h>
 float f(float x[])
{ float z;
  z=4*(x[0]-5)*(x[0]-5)+(x[1]-6)*(x[1]-6);
  return z;
}
float y,y1,y2,y3;
float x1[2], x2[2], x3[2],s[2],a[2],b[2];
int i,n=2;
huangjin() /*黄金分割法确定每轮最优值*/
{
  float e=0.01,h0=0.03,r=0.618,h=h0;  /*外推法确定区间*/
   {y1=f(x1);
    for(i=0;i<n;i++) s[i]=s[i]*h;
    for(i=0;i<n;i++) x2[i]=x1[i]+s[i];
    y2=f(x2);
   if(y2>y1)
    {for(i=0;i<n;i++) {s[i]=-s[i]; x3[i]=x1[i];y3=y1;}
     for(i=0;i<n;i++) {x1[i]=x2[i];y1=y2;  x2[i]=x3[i];  y2=y3;}
    }
   for(i=0;i<n;i++) {x3[i]=x2[i]+s[i]; y3=f(x3);}
   while(y3<y2)
     {for(i=0;i<n;i++) s[i]=s[i]*2;
      for(i=0;i<n;i++) {x1[i]=x2[i]; y1=y2; x2[i]=x3[i]; y2=y3; }
      for(i=0;i<n;i++) {x3[i]=x2[i]+s[i];y3=f(x3); }
     }
   for(i=0;i<n;i++)   {a[i]=x1[i];b[i]=x3[i];}
  }