c语言二分法判断函数的根 二分法求根c++程序( 二 )


#define DEFAULT_E(0.00000001)
#define _MID(x,y)((x+y)/2)
#define _VALUE(x)(2*x*x*x-4*x*x+3*x-6)
double _e;
int getRoot(double lower, double upper, double *result);
main()
{
double root;
printf("Enter a deviation:");
scanf("%lf", _e);
if(_e == 0.0)
_e = DEFAULT_E;
if(getRoot(DEFAULT_LOWER, DEFAULT_UPPER, root))
printf("Root:%2.8lf\n", root);
else
printf("Root:No Solution.\n");
}
int getRoot(double lower, double upper, double *result)
{
*result = _MID(lower,upper);
if(upper - lower = _e)
return 1;
if(_VALUE(lower)*_VALUE(*result) = 0)
return getRoot(lower, *result, result);
else if(_VALUE(*result)*_VALUE(upper) = 0)
return getRoot(*result, upper, result);
else
return 0;
}
很简单的递归 。
不过对于楼主的这道题,还要说一点 , 提示中的二分法求根,只能适用于在给定函数在给定区间中呈单调连续情况下才行得通哦~我也是按照给定的提示来写的上述代码 。
如果连续函数在给定区间不单调 , 很有可能中值*下界值和中值*上界值都大于0 , 那么会跳出认为没有根,而事实上很有可能这个中值点靠近函数极点 。
而真正用二分法求给定区间的思路是:
首先为函数求导,算出导函数的零点,然后再判断零点性质,最后将函数区间分为单调递增和单调递减间隔的形式,对每一段进行二分法求根 。
c语言二分法求方程的根的算法如果连续函数在给定区间不单调,很有可能中值*下界值和中值*上界值都大于0,那么会跳出认为没有根,而事实上很有可能这个中值点靠近函数极点 。
而真正用二分法求给定区间的思路是:
首先为函数求导,算出导函数的零点,然后再判断零点性质,最后将函数区间分为单调递增和单调递减间隔的形式,对每一段进行二分法求根 。
#include stdio.h
#include math.h
#define DEFAULT_UPPER(10)
#define DEFAULT_LOWER(-10)
#define DEFAULT_E(0.00000001)
#define _MID(x,y)((x+y)/2)
#define _VALUE(x)(2*x*x*x-4*x*x+3*x-6)
double _e;
int getRoot(double lower, double upper, double *result);
main()
{
double root;
printf("Enter a deviation:");
scanf("%lf", _e);
if(_e == 0.0)
_e = DEFAULT_E;
if(getRoot(DEFAULT_LOWER, DEFAULT_UPPER, root))
printf("Root:%2.8lf\n", root);
else
printf("Root:No Solution.\n");
}
int getRoot(double lower, double upper, double *result)
{
*result = _MID(lower,upper);
if(upper - lower = _e)
return 1;
if(_VALUE(lower)*_VALUE(*result) = 0)
return getRoot(lower, *result, result);
else if(_VALUE(*result)*_VALUE(upper) = 0)
return getRoot(*result, upper, result);
else
return 0;
}
【c语言二分法判断函数的根 二分法求根c++程序】关于c语言二分法判断函数的根和二分法求根c++程序的介绍到此就结束了,不知道你从中找到你需要的信息了吗 ?如果你还想了解更多这方面的信息,记得收藏关注本站 。

推荐阅读