| 查看: 1099 | 回复: 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 ] |
回复此楼» 猜你喜欢
建议基金发布提前给出明确的时间点
已经有10人回复
-
科研孤儿太难了
已经有17人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有19人回复
-
看来今天不会放榜了？
已经有11人回复
-
时间戳又变了
已经有19人回复
-
filecode，4个jtjc了
已经有14人回复
-
让我中一个面上吧！
已经有13人回复
-
人气不行了
已经有7人回复
-
感觉是下周放榜了
已经有12人回复
-
我面上完蛋了
已经有8人回复
cxymath
铁杆木虫 (知名作家)
- 应助: 0 (幼儿园)
- 金币: 19886.4
- 红花: 2
- 沙发: 1
- 帖子: 5865
- 在线: 957.5小时
- 虫号: 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
40