erfenfa.txt

二分法解非线性方程

(1)二分法求解非线性方程：
#include<stdio.h>
#include <math.h>
#define f(x) ((x*x-1)*x-1)  
void main()
{  float a,b,x,eps;
   int k=0;
   printf("intput eps\n");/* 容许误差 */
   scanf("%f",&eps);
   printf("a,b=\n");
   for(;;)
   {scanf("%f, %f",&a ,&b);
   if(f(a)*f(b)>=0) /* 判断是否符合二分法使用的条件 */
    printf("二分法不可使用,请重新输入:\n"); 
   else break;
   }
   do
   {  x=(a+b)/2;
      k++;
      if(f(a)*f(x)<0) /* 如果f(a)*f(x)<0，则根在区间的左半部分 */
         b=x;
      else if(f(a)*f(x)>0)  /* 否则根在区间的右半部分 */
        a=x;
	  else break;
   }while(fabs(b-a)>eps);/*判断是否达到精度要求,若没有达到,继续循环*/
   x=(a+b)/2; /* 取最后的小区间中点作为根的近似值 */
   printf("\n The root is x=%f, k=%d\n",x,k);
}
运行结果：
intput eps
0.00001
a,b=
2,-5

 The root is x=1.324721, k=20
Press any key to continue