程序1.txt

黄金分割 黄金分割 黄金分割

1.程序：
main()
{
float f(float t);
float absf(float x);
float a=-10.0,b=3.0;
float a1,a2,b1,b2,y1,y2,t;
float p=0.618,e=0.00001;

a1=b-p*(b-a);
a2=a+p*(b-a);
y1=f(a1);
y2=f(a2);

while(b-a>e*absf(b))
{if (y1>y2) {a=a1;
              a1=a2;
              y1=y2;
              a2=a+p*(b-a);
              y2=f(a2);
             }
        else {b=a2;
              a2=a1;
              y2=y1;
              a1=b-p*(b-a);
              y1=f(a1);
	     }
 }

t=(a+b)/2;
printf("%f  %f",t,f(t));
}

float absf(float x)
{if (x<0) return -x;
     else return x;
}

float f(float t)
{return t*t+10*t-10.0;}