| 查看: 489 | 回复: 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); } |
回复此楼» 猜你喜欢
压汞仪和BET测气凝胶孔隙率
已经有4人回复
-
博士申请都是内定的吗？
已经有14人回复
-
博士申请
已经有3人回复
-
投稿返修后收到这样的回复，还有希望吗
已经有5人回复
-
谈谈两天一夜的“延安行”
已经有13人回复
-
氨基封端PDMS和HDI反应快速固化
已经有11人回复
-
之前让一硕士生水了7个发明专利，现在这7个获批发明专利的维护费可从哪儿支出哈？
已经有11人回复
-
论文投稿求助
已经有4人回复
-
Applied Surface Science 这个期刊。有哪位虫友投过的能把word模板发给我参考一下嘛
已经有3人回复
-
投稿精细化工
已经有6人回复
wpwupingwp
木虫 (正式写手)
- 应助: 98 (初中生)
- 金币: 3109.8
- 散金: 250
- 红花: 5
- 帖子: 464
- 在线: 125.6小时
- 虫号: 2026298
- 注册: 2012-09-24
- 性别: GG
- 专业: 生物信息学
2楼2012-11-18 17:29:48
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 |
3楼2013-01-29 20:09:10