24小时热门版块排行榜

返回列表

路过081001

铜虫 (正式写手)

应助: 2 (幼儿园)
金币: 42
散金: 332
帖子: 515
在线: 59.1小时
虫号: 1752819
注册: 2012-04-13
性别: GG
专业: 通信理论与系统

[求助] 求高手，，C语言弦截法方程的根,我同谭浩强的视频程序几乎没有差别，为什么求不出呢

#include
#include

float f(float x)
{
return ((x-5)*x+16)*x-80;
}

float xpoint(float x1,float x2)
{
return (x1*f(x2)-x2*f(x1))/(f(x2)-f(x1));

}

float root(float x1,float x2)
{
  float x,y,y1,y2;
  y1=f(x1);y2=f(x2);
  do
  {
x=xpoint(x1,x2);
y=f(x);
if(y*y1>0) {y1=y;x1=x;}
else {x2=x;y2=y;}
  } while (fabs(y)>0.00001);
  return x;
}

void main()
{
float x,x1,x2,y1,y2;

  do
  { printf("input x1,x2:"

;
scanf("%f,%f",&x1,&x2);
y1=f(x1);
y2=f(x2);
} while(y1*y2>0);

x=root(x1,x2);

printf("A root is %f\n",x);

}

回复此楼

» 猜你喜欢

不要光上人人QQ，不然好没出息。

1楼 2012-11-18 14:15:19

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

wpwupingwp

木虫 (正式写手)

应助: 98 (初中生)
金币: 3109.8
散金: 250
红花: 5
帖子: 464
在线: 125.6小时
虫号: 2026298
注册: 2012-09-24
性别: GG
专业: 生物信息学

【答案】应助回帖

感谢参与，应助指数 +1

http://blog.csdn.net/zh458042230/article/details/7748198

赞一下

回复此楼

南无观世音菩萨

2楼2012-11-18 17:29:48

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

chunyuan314

银虫 (正式写手)

应助: 10 (幼儿园)
金币: 1861.4
散金: 123
红花: 2
帖子: 355
在线: 216.4小时
虫号: 2016643
注册: 2012-09-20
性别: GG
专业: 电化学

LZ，问题出在root函数 do-while循环部分。
按弦截法，root函数可以这样：
float root(float x1,float x2)
{
  float x,y,y1,y2;
  y1=f(x1);y2=f(x2);
x=xpoint(x1,x2);
y=f(x);
if(fabs(y)<0.00001){ /* 如果精度达到，Okay */
   return x;
}else if(y*y1>0){  /*即如果新的y与y1同号，则改为在x~x2之间用弦截法*/
   return root(x,x2);
}else{
   return root(x1,x);  /*否则就在在x1~x之间用弦截法*/
}
}

运行root(-10,10)可得到5

函数图形.png

赞一下

回复此楼

Lifeisgood

3楼2013-01-29 20:09:10

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主路过081001 的主题更新

返回列表