| 查看: 762 | 回复: 8 | |||
| 当前主题已经存档。 | |||
[交流]
【求助】数值分析计算机应用c语言问题
|
|||
|
#include "math.h" main() {float x1,x2,f1,f2; x1=1.9; f1=pow(x1,7)-28*pow(x1,4)+14; f2=7*pow(x1,6)-4*28*pow(x1,3); x2=x1-f1/f2; do {x1=x2; f1=pow(x1,7)-28*pow(x1,4)+14; f2=7*pow(x1,6)-4*28*pow(x1,3); x2=x1-f1/f2; if(f2+1.0=1.0); {printf("**end**"\n);} } while(fabs(x1-x2)>1e-5); printf("x=%f\n",x2); } 题目是用newton法求方程，x（7）-28*x（4）+14=0 （）内为x的指数 在（0.1,1,9）内有近似根，用6次迭代，使误差小于0.00001 我的程序在上面，用TC运行，x1初值为0.1时，输出为 **end** 共30行 x=3.030577 x1初值为1.9时，输出为 **end** 共5行 x=0.845497 与我的设想不同，end为什么会有这么多行， 谁能帮忙改一下，谢谢各位大大 [ Last edited by wuguocheng on 2009-10-11 at 11:51 ] |
回复此楼» 猜你喜欢
所感
已经有3人回复
-
要不要辞职读博？
已经有7人回复
-
不自信的我
已经有11人回复
-
北核录用
已经有3人回复
-
实验室接单子
已经有3人回复
-
磺酰氟产物，毕不了业了！
已经有8人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有10人回复
-
26申博（荧光探针方向，有机合成）
已经有4人回复
-
论文终于录用啦！满足毕业条件了
已经有26人回复
-
2026年机械制造与材料应用国际会议 （ICMMMA 2026）
已经有4人回复
cxymath
铁杆木虫 (知名作家)
- 应助: 0 (幼儿园)
- 金币: 18997.9
- 红花: 2
- 沙发: 1
- 帖子: 5845
- 在线: 951.6小时
- 虫号: 787220
- 注册: 2009-06-04
- 专业: 教育统计学
2楼2009-10-11 12:29:21
3楼2009-10-11 14:45:18
4楼2009-10-11 19:39:36
5楼2009-10-11 19:55:50
formleaf
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 贵宾: 1.097
- 金币: 3383.1
- 散金: 2780
- 红花: 3
- 帖子: 991
- 在线: 69.6小时
- 虫号: 698652
- 注册: 2009-02-09
- 性别: GG
- 专业: 组合数学
★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
wuguocheng(金币+5,VIP+0): 感谢参与, 以后多支持数学版块 10-15 19:50
小木虫(金币+0.5):给个红包，谢谢回帖交流
wuguocheng(金币+5,VIP+0): 感谢参与, 以后多支持数学版块 10-15 19:50
|
修改后的程序 #include "stdio.h" #include "math.h" main() {double x1,x2,f1,f2; x1=1.9; f1=pow(x1,7)-28*pow(x1,4)+14; f2=7*pow(x1,6)-4*28*pow(x1,3); x2=x1-f1/f2; do {x1=x2; f1=pow(x1,7)-28*pow(x1,4)+14; f2=7*pow(x1,6)-4*28*pow(x1,3); x2=x1-f1/f2; if(f2+1.0==1.0)printf("**end**\n"  ;} while(fabs(x1-x2)>1e-5); printf("x=%f\n",x2); } 这样就正常了，老兄， if(f2+1.0==1.0)后面多打了一个； |
6楼2009-10-15 14:19:00
7楼2009-10-15 15:02:32
8楼2009-12-01 21:04:43
9楼2009-12-03 22:30:35