C语言用二分法求方程的近似解的方法

C语言用二分法求方程的近似解的方法

基础知识以下地址：http://jingyan.baidu.com/article/597a06438def54312a524376.html



我的应用如下：求x^5+x^3+7=0,求取x的值是多少？

编译环境：Window8 +C Free
源程序如下：

#include "stdio.h"

#include "math.h"

#include "stdlib.h"

double func (double x)//这个函数主要实现求y=x^3+x^5+7 的值

{

    double y;

    y=pow(x,3)+pow(x,5)+7;       //函数表达式可自己定义

    return y;

}
double CalRoot(double m,double n)//这里用的迭代法和二分法的结合吧！

{//这个自定义函数就是不断的取y(m)和y(n)的值，和他们的中间值y(x)的比较

//不断取中间值y=func(x),当y接近于0时，那么结束循环！此时x的值便是我们求到

//的近似值！

     double x,y,y1,y2,a,b;

     a=m;  b=n;

     while (1)

     {

           x=(a+b)/2;

           y=func(x);

          if(fabs(y)<=1e-5) //相当于|y|<0.00001,当y为浮点型数据时，用这样的方式来判断与0的大小！          

    {

     return x;//这个if语句也是判断循环结束的标志！

    }

          y1=func(a);

          y2=func(b);

          if(y1*y<0)  b=x;//如果y1和y为异号的话就将a,b中间值x传递给b

          if(y2*y<0)  a=x;//如果y1和y为异号的话就将a,b中间值x传递给a

     }

}

void main ( )//原来程序中过多的部分我已经全部删除掉了，这版本比较简洁！

{

     double a,b,a1,b1,x;

     a=-3;  b=-1;       //a，b的值是搜索根的区间范围，可以适当调整,范围不可过大这里！

     a1=a;  b1=b;//还要注意func(a)和func(b)那么他们的值必须为异号才行！这是二分法使用的前提！

     x=CalRoot(a1,b1);//也就是说范围a,b的取值要注意很多这里我取得是a=-3;b=-1

     printf("X=%f\n",x);
 }