三分算法

double deal()
{
double l=0,r=PI,mid,mmid;
while(fabs(r-l)>EPS){
mid=(r+l)/2;
mmid=(mid+r)/2;
cal(mid)<cal(mmid)? r=mmid:l=mid;
// 这个地方要注意，并不是所有的地方都是这么用的
// 诀窍就是，如果mid更靠近极值点，那么r=mmid
}
return cal(l)*cal(l);
}