这个程序是用二分法求方程X的4次方－3x＋1在区间0.3到0.4的根，要求误差不超过0.2×10的负2次方

#include<stdio.h>
#include<math.h>
#define F(y) (y*y*y*y-3.0*y+1.0)
main()
{
double a,b,a1,b1,x;
int k;
double t=0.2e-2;
a=0.3;
b=0.4;
j=log(200*(b-a))/log(2);
n=(int)j
a1=(F(a)>F(b)?a:b);
b1=(F(a)<F(b)?a:b);/*把大的区间两端函数值放入a1，小的放入b1*/
for(k=1;k<=n;k++)
{
x=(a1+b1)/2.0;/*给X赋初值*/
if(F(x)==0) break;
else if(F(x)<0)
{
b1=x;
if (b1-a1<=t) break;
}
else if(F(x)>0)
{
a1=x;
if (b1-a1<=t) break;
}
}
printf("方程近似跟为%lf.4/n",x);
}