| 查看: 806 | 回复: 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人回复
-
一志愿南昌大学324求调剂
已经有7人回复
-
材料专硕322分
已经有5人回复
-
材料与化工（0856）304求B区调剂
已经有12人回复
-
材料调剂
已经有10人回复
-
省双一流重点一本大学招收调剂
已经有6人回复
-
321求调剂
已经有6人回复
-
288求调剂 一志愿哈工大 材料与化工
已经有26人回复
-
材料调剂
已经有5人回复
-
一志愿北交大材料工程总分358
已经有4人回复
cxymath
实习版主
- 应助: 0 (幼儿园)
- 金币: 19193.4
- 红花: 2
- 沙发: 1
- 帖子: 5861
- 在线: 955.9小时
- 虫号: 787220
- 注册: 2009-06-04
- 专业: 教育统计学
2楼2009-10-11 12:29:21
swc0339
无虫
天天向上
- 应助: 0 (幼儿园)
- 金币: 506.3
- 红花: 2
- 帖子: 289
- 在线: 6.8小时
- 虫号: 858237
- 注册: 2009-09-27
- 性别: GG
- 专业: 工业工程与管理
3楼2009-10-11 14:45:18
4楼2009-10-11 19:39:36
swc0339
新虫
天天向上
- 应助: 0 (幼儿园)
- 金币: 506.3
- 红花: 2
- 帖子: 289
- 在线: 6.8小时
- 虫号: 858237
- 注册: 2009-09-27
- 性别: GG
- 专业: 工业工程与管理
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
