| 查看: 808 | 回复: 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 ] |
回复此楼» 猜你喜欢
291求调剂
已经有25人回复
-
一志愿北京科技大学085601材料工程英一数二初试总分335求调剂
已经有6人回复
-
一志愿北京科技，085601总分305求调剂
已经有10人回复
-
化学工程专硕324分，一志愿中国矿业大学求调剂
已经有6人回复
-
279求调剂
已经有3人回复
-
26年申博
已经有4人回复
-
材料工程322分
已经有5人回复
-
318求调剂，计算材料方向
已经有7人回复
-
292分，材料与化工，申请调剂
已经有20人回复
-
085600 295分求调剂
已经有18人回复
cxymath
铁杆木虫 (知名作家)
- 应助: 0 (幼儿园)
- 金币: 19193.4
- 红花: 2
- 沙发: 1
- 帖子: 5861
- 在线: 955.9小时
- 虫号: 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
