lamename 2015-11-12 09:44 采纳率: 50%
浏览 1659

C语言牛顿迭代法,正常运行输出错误结果,求助

求方程ax^3+bx^2+c^x+d=0的解,其中xn+1=xn-f(xn)/f'(xn)
我的代码这样写的:
#include
#include
int main(void)
{
int a,b,c,d;
printf("Please input four integers:");
scanf("%d %d %d %d",&a,&b,&c,&d);
double x,y;
x=1;
while(1)
{
y=x;
x=(x-(a*pow(x,3)+b*pow(x,2)+pow(c,x)+d)/(3*a*pow(x,2)+2*b*x+pow(c,x)*log(c)));
if((fabs(x)-fabs(y))<1e-5)
{
printf("The solution is:%lf",x);
break;
}
}
return 0;
}

输出的 结果和正确答案不符,求助。。

  • 写回答

1条回答

  • ysuwood 2015-11-12 11:40
    关注
    
    
    

    由于log函数原型有:
    double log(double x);
    float log(float x);
    long double log(long double x);
    没有形参x是int型,所以需要对int 型c进行强制类型转换:
    x=(x-(a*pow(x,3)+b*pow(x,2)+pow(c,x)+d)/(3*a*pow(x,2)+2*b*x+pow(c,x)*log((float)c)));

    
    
    
    评论

报告相同问题?

悬赏问题

  • ¥15 乘性高斯噪声在深度学习网络中的应用
  • ¥15 运筹学排序问题中的在线排序
  • ¥15 关于docker部署flink集成hadoop的yarn,请教个问题 flink启动yarn-session.sh连不上hadoop,这个整了好几天一直不行,求帮忙看一下怎么解决
  • ¥30 求一段fortran代码用IVF编译运行的结果
  • ¥15 深度学习根据CNN网络模型,搭建BP模型并训练MNIST数据集
  • ¥15 C++ 头文件/宏冲突问题解决
  • ¥15 用comsol模拟大气湍流通过底部加热(温度不同)的腔体
  • ¥50 安卓adb backup备份子用户应用数据失败
  • ¥20 有人能用聚类分析帮我分析一下文本内容嘛
  • ¥30 python代码,帮调试,帮帮忙吧