二分法.txt

来自「数值计算方法所用的上机实验C++源程序,个人所编.」· 文本 代码 · 共 24 行

#include <math.h>
#define f(x) ((x*x+10)*x-20)
#define eps 0.0001  /* 容许误差 */
main()
{  float a,b,y,x;
   int k=0;
   printf("a,b=");
   scanf("%f, %f ",&a ,&b);
   if(f(a)*f(b)>=0) /* 判断是否符合二分法使用的条件 */
   {  printf("Not bisect of bisect");   return;   }
   do
   {  x=(a+b)/2;
      k++;
      if(f(a)*f(x)<0) /* 如果f(a)*f(x)<0，则根在区间的左半部分 */
         b=x;
      else  /* 否则根在区间的右半部分 */
        a=x;
   }while(fabs(b-a)>eps);/*判断是否达到精度要求,若没有达到,继续循环*/
   x=(a+b)/2; /* 取最后的小区间中点作为根的近似值 */
   printf("\n The root is x=%f, k=%d\n",x,k);
}
     a,b=1↙
     The root is x=1.594574, k=14