应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 计算数学 » 【求助】数值分析计算机应用c语言问题
91/1返回列表
查看: 749  |  回复: 8
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

adjoinff

银虫 (小有名气)

[交流] 【求助】数值分析计算机应用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 ]
回复此楼

» 猜你喜欢
1楼2009-10-11 11:34:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cxymath

铁杆木虫 (知名作家)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
要找计算数学方面的人了
赞一下(1人) 回复此楼
2楼2009-10-11 12:29:21
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

swc0339

捐助贵宾 (小有名气)

天天向上
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
点点风(金币+2,VIP+0):感谢您的积极应助 10-11 21:10
引用回帖:
Originally posted by adjoinff at 2009-10-11 11:34:
#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;
   ...

if(f2+1.0=1.0);
应该改为
if(f2+1.0==1.0);

你写的是一条赋值语句，所以每次都输出END
赞一下(11人) 回复此楼
锻压青春 焊接理想 处理未来 铸造辉煌
3楼2009-10-11 14:45:18
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

adjoinff

银虫 (小有名气)
点点风(金币+0,VIP+0):问题解决了么？ 10-11 21:10
if(f2+1.0==1.0);
我就是这么写的，上面是少打了一个
赞一下(10人) 回复此楼
4楼2009-10-11 19:39:36
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

swc0339

捐助贵宾 (小有名气)

天天向上
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
引用回帖:
Originally posted by adjoinff at 2009-10-11 19:39:
if(f2+1.0==1.0);
我就是这么写的，上面是少打了一个

那你还少打了什么。。。。。。
赞一下(1人) 回复此楼
锻压青春 焊接理想 处理未来 铸造辉煌
5楼2009-10-11 19:55:50
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

formleaf

木虫 (正式写手)
★ ★ ★ ★ ★ ★
小木虫(金币+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)后面多打了一个；
赞一下(11人) 回复此楼
6楼2009-10-15 14:19:00
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

adjoinff

银虫 (小有名气)
多谢这位仁兄
回复此楼
7楼2009-10-15 15:02:32
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yaoyusui

木虫 (正式写手)
那要问计算数学专业的人了
赞一下 回复此楼
8楼2009-12-01 21:04:43
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

403968649

新虫 (初入文坛)
我还以为是程序问题，看啦半天，
赞一下 回复此楼
9楼2009-12-03 22:30:35
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 adjoinff 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭