编程小白求助，计算牛顿迭代法，为什么编出来的程序输出为 -nan

Earl_Flo

#include <stdio.h>
#include <math.h>

int main(){
double newton(double a, double b, double c,double x_0);
double x_n = newton(1,1,0.25,-1);
printf("x_n is:%lf",x_n);
return 0;
}


double newton(double a, double b, double c, double x_0) {

double x_n = x_0;
double f,fd;
f = a*x_n*x_n+b*x_n + c;
  while(fabs(f)>0.001){
   f = a*x_n*x_n+b*x_n+c;
   fd = 2*a*x_n +b;
   x_n = x_n + f/fd;
   }
  return x_n;
}

jackz007

       你的 "牛顿迭代法" 是用来算什么的？

Earl_Flo

#include <stdio.h>
#include <math.h>




int main(){
double newton(double a, double b, double c,double x_0);
double x_n = newton(1,1,0.25,-1);
printf("x_n is:%lf",x_n);
return 0;
}


double newton(double a, double b, double c, double x_0) {

  double x_n=x_0;
  double f,fd;
  do{
   x_0=x_n;
   f = a*x_0*x_0+b*x_0+c;
   fd = 2*a*x_0 + b;
   x_n = x_0-f/fd;
   }while(fabs(f)>0.001);
    return x_n;

}


这个是正确的程序 但是不明白为什么要x_0=x_n

Earl_Flo

#include <stdio.h>
#include <math.h>




int main(){
double newton(double a, double b, double c,double x_0);
double x_n = newton(1,1,0.25,-1);
printf("x_n is:%lf",x_n);
return 0;
}


double newton(double a, double b, double c, double x_0) {

  double x_n=x_0;
  double f,fd;
  while(fabs(f)>0.001){
   x_0=x_n;
   f = a*x_0*x_0+b*x_0+c;
   fd = 2*a*x_0 + b;
   x_n = x_0-f/fd;
   }
    return x_n;

}


想先判断 再执行 把do while换成了while 输出的值从正确的 -0.51.... 变成了-1.0000

jackz007

Earl_Flo 发表于 2022-11-7 21:17
#include
#include

1. #include <stdio.h>
   
2. 
   
3. #define abs(n) ((n) >= 0 ? (n) : -(n))
   
4. 
   
5. double newton(double a , double b , double c , double x0)
   
6. {
   
7.          double f , fd , x1 = x0                        ;
   
8.          do {
   
9.                  x0 = x1                                ;
   
10.                  f = a * x0 * x0 + b * x0 + c           ;
    
11.                  fd = 2 * a * x0 + b                    ;
    
12.                  x1 = x0 - f / fd                       ;
    
13.         } while(abs(x1 - x0) > 1e-15)                   ;
    
14.         return x1                                       ;
    
15. }
    
16. 
    
17. int main(void)
    
18. {
    
19.         printf("x = %lf\n" , newton(1 , 1 , 0.25 , -1)) ;
    
20. }

复制代码

        编译、运行实况：

1. D:\[00.Exerciese.2022]\C>g++ -o x x.c
   
2. 
   
3. D:\[00.Exerciese.2022]\C>x
   
4. x = -0.500000
   
5. 
   
6. D:\[00.Exerciese.2022]\C>

复制代码

Earl_Flo

jackz007 发表于 2022-11-7 21:54
编译、运行实况：

谢谢回答 但是回答和我提的问题不太一样

jackz007

Earl_Flo 发表于 2022-11-7 22:41
谢谢回答 但是回答和我提的问题不太一样

         有什么不一样？

Earl_Flo

jackz007 发表于 2022-11-7 22:41
有什么不一样？

正确代码我已经的出来了 我只是想知道我的另外两个和正确的那个有什么不一样 第二个代码是个正确的  我的题目是当f<=0.001时就结束程序了，而且也不需要abs n。 我就是想弄明白为什么我写的不对 我想使用while语句 你写的这个程序之前在网上看见过了  谢谢

jackz007

Earl_Flo 发表于 2022-11-7 22:56
正确代码我已经的出来了 我只是想知道我的另外两个和正确的那个有什么不一样 第二个代码是个正确的  我的 ...

       你只是笼统说了一个 "牛顿迭代法" ，至于计算什么也不知道，我辛辛苦苦在网上查资料，才猜出你是想求一元二次方程的根，然后，才参考别人的算法原理（网上没有找到计算一元二次方程根的代码），比照你的参数写出了这个代码，我才原创的东西，你却说在网上见过？？？

Earl_Flo

jackz007 发表于 2022-11-7 23:04
你只是笼统说了一个 "牛顿迭代法" ，至于计算什么也不知道，我辛辛苦苦在网上查资料，才猜出你是 ...

我也是第一次发帖 可能这个题比较基础 我再csdn上看到很多极度相似的 如有冒犯 很抱歉