hft.c

恩威白垩纪我不而不忘记而不回家我回家二万二万而或加我喝酒钨合金我就和vwh而惊闻噩耗我

#include <stdio.h>
#define F(x) (x*x+5*x+5)

main ()
{ double a=-10.0,b=1.0,e=0.0000001,h=0.00001,x1=-1.0,x2,x3;
  x2=x1+h;
  if (F(x1)<F(x2)) {h=-h;x3=x1;x1=x2;x2=x3;x3=x2+h;}
    else h=2*h,x3=x2+h;
  while (F(x2)>=F(x3)) 
      x1=x2,x2=x3,x3=x2+h;	
  a=(x1<x3)? x1:x3,b=(x1>x3)? x1:x3;
  while (b-a>e)	  
  { x2=a+0.618*(b-a);
    x1=b-0.618*(b-a);
    if (F(x1)>F(x2)) a=x1;
      else b=x2;
  }
  a=(a+b)/2;
  printf ("最少值的x为：%f  最少值为：%f",a,F(a));
   
}