| 查看: 595 | 回复: 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); } |
回复此楼» 猜你喜欢
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有6人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有5人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有7人回复
-
浙江师范大学是怎么坑我的
已经有7人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有5人回复
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
