黄金分割法（）.txt

#include <stdio.h>
#include <math.h>
double fun(double x)
{
return (4*x*x*x-6*x+2);//函数定义
}

main()
{
double y[5],x[5];
double c=0,a=0,b=0,h=1;//定义c,h(为x初值，h为歩长）
double R=0.618,g=0.0001;        //定义缩短率,精度
y[1]=fun(c);x[1]=c;
x[2]=(c+h);
y[2]=fun(x[2]);
if(y[2]>y[1])
{
h=-h;x[3]=x[1];y[3]=y[1]; 
x[1]=x[2];y[1]=y[2];
x[2]=x[3];y[2]=y[3];   
}
x[3]=(x[2]+h);y[3]=fun(x[3]);
while(y[3]<y[2])
{
h=2*h;
x[1]=x[2]; y[1]=y[2];
x[2]=x[3]; y[2]=y[3];
x[3]=x[2]+h; y[3]=fun(x[3]);
}
if(x[3]<x[1]){a=x[3];b=x[1];}