
用二分法求方程的根(C语言)
想问下大家,如何求2*x*x*x-4*x*x+3*x-6=0在(-10,10)内的全部的根,要是全部的根,不能只求一个,大家帮帮忙,帮忙分析下,并给出源代码,谢谢大家了啊...
想问下大家,如何求2*x*x*x-4*x*x+3*x-6=0在(-10,10)内的全部的根,要是全部的根,不能只求 一个,大家帮帮忙,帮忙分析下,并给出源代码 ,谢谢大家了啊!
展开
展开全部
刚才那个只能找一个根,现在改了下,可以找所有根:
刚才有个小BUG,改了下::
=================================================
#include <math.h>
void FindRoot(void)
{
double dbLeftX = -10.0;
double dbRightX = 10.0;
double dbResultY = 0.0;
double dbTmp = 0.0;
double dbMidleX =0.0;
//当二分点Y的值不为0不是根(即,二分点X不满足方程)
//当区间无限小时完成计算(下面的语句表示当区间不是无限小时继续计算)
while ( abs(dbRightX-dbLeftX) > 0.1e-6 )
{
//二分区间
dbMidleX = ( dbLeftX + dbRightX ) / 2.0;
//算二分点Y的值
dbResultY = 2 * pow(dbMidleX,3) - 4 * pow(dbMidleX,2) + 3 * dbMidleX - 6;
//如果计算的Y值为0时,找到一个根
if ( abs(dbResultY) < 0.1e-6 )
{
printf("Finde a root in [-10,10] , value = [%lf].\n",dbMidleX);
}
//判断方程的根落于刚才二分区间的哪个区间,并生成新的区间.
if ( dbResultY < 0.0 )
{
dbTmp = 2 * pow(dbLeftX,3) - 4 * pow(dbLeftX,2) + 3 * dbLeftX - 6;
if ( dbTmp < 0.0 )
{
dbLeftX = dbMidleX;
}
else
{
dbRightX = dbMidleX;
}
}
else
{
dbTmp = 2 * pow(dbLeftX,3) - 4 * pow(dbLeftX,2) + 3 * dbLeftX - 6;
if ( dbTmp > 0.0 )
{
dbLeftX = dbMidleX;
}
else
{
dbRightX = dbMidleX;
}
}
}
//区间无限小,如果区间两端计算的Y值一个为正一个为负则区间中点也是一个根
if( 2 * pow(dbLeftX,3) - 4 * pow(dbLeftX,2) + 3 * dbLeftX - 6 < 0 &&
2 * pow(dbRightX,3) - 4 * pow(dbRightX,2) + 3 * dbRightX - 6 >0 )
{
printf("Finde a root in [-10,10] , value = [%lf].\n",dbMidleX);
}
}
刚才有个小BUG,改了下::
=================================================
#include <math.h>
void FindRoot(void)
{
double dbLeftX = -10.0;
double dbRightX = 10.0;
double dbResultY = 0.0;
double dbTmp = 0.0;
double dbMidleX =0.0;
//当二分点Y的值不为0不是根(即,二分点X不满足方程)
//当区间无限小时完成计算(下面的语句表示当区间不是无限小时继续计算)
while ( abs(dbRightX-dbLeftX) > 0.1e-6 )
{
//二分区间
dbMidleX = ( dbLeftX + dbRightX ) / 2.0;
//算二分点Y的值
dbResultY = 2 * pow(dbMidleX,3) - 4 * pow(dbMidleX,2) + 3 * dbMidleX - 6;
//如果计算的Y值为0时,找到一个根
if ( abs(dbResultY) < 0.1e-6 )
{
printf("Finde a root in [-10,10] , value = [%lf].\n",dbMidleX);
}
//判断方程的根落于刚才二分区间的哪个区间,并生成新的区间.
if ( dbResultY < 0.0 )
{
dbTmp = 2 * pow(dbLeftX,3) - 4 * pow(dbLeftX,2) + 3 * dbLeftX - 6;
if ( dbTmp < 0.0 )
{
dbLeftX = dbMidleX;
}
else
{
dbRightX = dbMidleX;
}
}
else
{
dbTmp = 2 * pow(dbLeftX,3) - 4 * pow(dbLeftX,2) + 3 * dbLeftX - 6;
if ( dbTmp > 0.0 )
{
dbLeftX = dbMidleX;
}
else
{
dbRightX = dbMidleX;
}
}
}
//区间无限小,如果区间两端计算的Y值一个为正一个为负则区间中点也是一个根
if( 2 * pow(dbLeftX,3) - 4 * pow(dbLeftX,2) + 3 * dbLeftX - 6 < 0 &&
2 * pow(dbRightX,3) - 4 * pow(dbRightX,2) + 3 * dbRightX - 6 >0 )
{
printf("Finde a root in [-10,10] , value = [%lf].\n",dbMidleX);
}
}
本回答被提问者和网友采纳
已赞过
已踩过<
评论
收起
你对这个回答的评价是?
展开全部
^是位运算符
不是多少次方
改成q1=pow(10,-3),q2=pow(10,-6);
不是多少次方
改成q1=pow(10,-3),q2=pow(10,-6);
已赞过
已踩过<
评论
收起
你对这个回答的评价是?
推荐律师服务:
若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询