24小时热门版块排行榜

>论坛更新日志 (4546)
>文献求助 (475)
>导师招生 (389)
>虫友互识 (350)
>考博 (167)
>论文投稿 (159)
>博后之家 (139)
>休闲灌水 (131)
>硕博家园 (116)
>考研 (96)
>公派出国 (92)
>基金申请 (79)
>论文道贺祈福 (71)
>教师之家 (71)
>找工作 (54)
>招聘信息布告栏 (45)

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

路过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的回帖

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