三分法

double lm,rm;
while(r-l>eps)//控制精度
{
//cout<<l<<" “<<r<<endl;
lm=l+(r-l)/3;//左三等分点
rm=r-(r-l)/3;//右三等分点
if(f( lm )<=f( rm )) l=lm;//极大值在lm左侧
else r=rm;//反之在rm右侧
}
l = f(l) > f®?l:r;
//printf(”%.5lf",l);
cout.precision(5);
cout << l << endl;*/