c - 牛顿-拉夫逊算法

Question

我要编译并运行此程序以使用 Newton-Raphson 算法求解方程，但 do-while 循环中的迭代次数始终为 1，事实并非如此。

代码

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

float fonc(float x){
   float result;
   result=2*(pow(x,3))-x-2;
   return result;
 }

float foncprime(float x){
   float result;
   result= 6*(pow(x,2))-1;
   return result;
}

int main(void)
{
   long double x1,x0;
   float eps;
   eps=0.5*pow(10,-4);
   x0=1;
   x1=1;
   int i=0;
   do
   {
      x0=x1;
      x1=x0-(fonc(x0)/foncprime(x0));
      i++;

   }
   while(x1-x0>eps);
   printf("%d",i);
}

Answer 1

您可能应该检查 while (fabs(x1-x0) > eps)。