24小时热门版块排行榜    

查看: 520  |  回复: 2
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

路过081001

铜虫 (正式写手)

[求助] 求高手,,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,不然好没出息。
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chunyuan314

银虫 (正式写手)

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的回帖
查看全部 3 个回答

wpwupingwp

木虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
南无观世音菩萨
2楼2012-11-18 17:29:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料与化工 322求调剂 +5 然11 2026-03-19 5/250 2026-03-25 16:26 by 功夫疯狂
[考研] 333求调剂 +6 87639 2026-03-21 11/550 2026-03-25 16:17 by 87639
[考研] 274求调剂 +7 顾九笙要谦虚 2026-03-24 7/350 2026-03-25 15:18 by baoball
[考研] 284求调剂 +15 Zhao anqi 2026-03-22 15/750 2026-03-25 12:51 by wht0531
[考研] 282求调剂 +3 wcq131415 2026-03-24 3/150 2026-03-25 12:16 by userper
[考研] 求调剂,一志愿:南京航空航天大学大学 ,080500材料科学与工程学硕,总分289分 +6 @taotao 2026-03-19 6/300 2026-03-25 08:37 by 木托莫露露
[考研] 上海电力大学材料防护与新材料重点实验室招收调剂研究生(材料、化学、电化学,环境) +4 我爱学电池 2026-03-23 4/200 2026-03-25 00:59 by 1027_324
[考研] 求调剂 +5 林之夕 2026-03-24 5/250 2026-03-24 17:16 by dick_runner
[考博] 申博26年 +4 八6八68 2026-03-19 4/200 2026-03-24 15:49 by 小Ben呵呵
[考研] 307求调剂 +3 余意卿 2026-03-21 6/300 2026-03-24 15:03 by 余意卿
[考研] 305分求调剂(食品工程) +5 Sxy112 2026-03-21 7/350 2026-03-24 12:27 by 544594351
[考研] 化学308分求调剂 +3 你好明天你好 2026-03-23 3/150 2026-03-23 20:11 by macy2011
[论文投稿] 急发核心期刊论文 +3 贤达问津 2026-03-23 5/250 2026-03-23 17:13 by 妹子不好惹
[考研] 328求调剂 +4 LHHL66 2026-03-23 4/200 2026-03-23 14:55 by lbsjt
[考研] 初试 317 +7 半拉月丙 2026-03-20 7/350 2026-03-21 22:26 by peike
[考研] 0703化学调剂 +4 妮妮ninicgb 2026-03-21 4/200 2026-03-21 18:39 by 学员8dgXkO
[基金申请] 学校已经提交到NSFC,还能修改吗? 40+4 babangida 2026-03-19 9/450 2026-03-21 16:12 by babangida
[考研] 一志愿 南京航空航天大学大学 ,080500材料科学与工程学硕 +5 @taotao 2026-03-20 5/250 2026-03-20 20:16 by JourneyLucky
[考研] 261求B区调剂,科研经历丰富 +3 牛奶很忙 2026-03-20 4/200 2026-03-20 19:34 by JourneyLucky
[考研] 本科郑州大学物理学院,一志愿华科070200学硕,346求调剂 +4 我不是一根葱 2026-03-18 4/200 2026-03-19 09:11 by 浮云166
信息提示
请填处理意见