Question

C语言牛顿迭代法求2x^3-4x^2+3x-6=0在1.5附近的根

#include<stdio.h>
#include<math.h>
void main()
{
float b,c,y,y1,root,y2,x;
x=1.5;
do
{
b=pow(x,3);
c=sqrt(x);
y=b*2-4*c+3*x-6;
y1=6*c-8*x+3;
root=x; ①
x=x-y/y1; ②
y2=root-x;
}
while(fabs(y2)>1e-5);
printf("%f\n",root);
}

编译和链接都通过了，就运行的时候显示溢出（一段英文）。。
当我把1式和2式对调位置后,运行有结果。。但貌似结果不对。。
跪求解脱。。大虾们，帮帮忙吧~~~

Answer 1

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

float Newton(float a,float b,float c,float d)
{
float x=1,x0,f,f1;
do
{
x0=x;
f=((a*x0+b)*x0+c)*x0+d;
f1=(3*a*x0+2*b)*x0+c;
x=x0-f/f1;
}
while(fabs(x-x0)>=1e-5);
return x;
}
int main()
{
float a,b,c,d;
scanf("%f%f%f%f\n",&a,&b,&c,&d);
printf("%.2fx^3+%.2fx^2+%.2fx+%.2f=0\n",a,b,c,d);
printf("x=%.5f\n",Newton(a,b,c,d));
return 0;
}
我写了一个通用的，键盘输入系数从高次到低次，
你参考一下。

Answer 2

c = sqrt(x); 这句是错的，sqrt是求根，应该是c = pow(x, 2);