| 查看: 568 | 回复: 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); } |
回复此楼» 猜你喜欢
看来今天不会放榜了？
已经有5人回复
-
我面上完蛋了
已经有8人回复
-
只有每年这种时候来逛逛小木虫
已经有25人回复
-
什么时候开奖？
已经有7人回复
-
时间戳又变了
已经有18人回复
-
时间戳又变了8-15
已经有27人回复
-
今天基金会出结果吗？20260819
已经有16人回复
-
科研孤儿太难了
已经有16人回复
-
应该是下周三26日公布了吧？
已经有4人回复
-
投稿咨询
已经有7人回复
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
40