#include <stdio.h>
#include <math.h>
int main()
{
double x0, x1; // 定义迭代初值x0和迭代结果x1
scanf("%lf", &x0); // 从键盘输入x0
do
{
x1 = x0 - (2 * pow(x0, 3) - 4 * pow(x0, 2) + 3 * x0 - 6) / (6 * pow(x0, 2) - 8 * x0 + 3); // 根据牛顿迭代法计算x1
if (fabs(x1 - x0) < 1e-6) // 如果前后两次求出的x的差的绝对值小于10-6,跳出循环
{
break;
}
x0 = x1; // 否则,将x1赋值给x0,继续迭代
} while (1);
printf("%.6lf\n", x1); // 输出方程的根,保留6位小数,末尾换行
return 0;
}
Problem Description
用牛顿迭代法求下面方程在输入初值点附近的根:
2x3-4x2+3x-6=0
要求前后两次求出的x的差的绝对值小于10-6
牛顿迭代法公式如下:
将给定给定方程写成f(x)=0的形式,在给定初值x0的情况下,按如下公式迭代计算:
xn+1=xn-f(x)/f'(x)
提示:C语言数学库中有求指数an的函数pow(a, n)以及求x绝对值的函数fabs(x)
浮点型数据请定义为双精度double类型。
Input Description
双精度浮点型初值x0
Output Description
x0附件方程的根,小数点后保留6位小数,末尾换行。
Sample Input
3
Sample Output
2.000000