用二分法求方程 x的五次方－3x＋1＝0在(0,1)上的近似解,精确到c=0.001,写出算法和

#include
#include
#include
double calculate(double middle) //将值代入表达式求表达式的结果
{
\x09return (pow(middle,5) - 3 * middle + 1);
}
int main(void)
{
\x09double lower = 0.0;
\x09double high = 1.0;
\x09double middle1 = (lower + high) / 2.0;
\x09double middle2 = 0;
\x09while (abs(middle1 - middle2) >= 0.001) //如果前后两次中值之差 >= 0.001,就执行循环
\x09{
\x09\x09if (calculate(middle1) > 0) //表达式的值大于0,将其中值middle1赋给high
\x09\x09{
\x09\x09\x09high = calculate(middle1);
\x09\x09}
\x09\x09else
\x09\x09\x09lower = middle1;
\x09\x09middle2 = middle1;
\x09\x09middle1 = (lower + high) / 2.0;
\x09}
\x09printf("%lg\n",middle1);
\x09return 0;
}