用C语言写一段用牛顿迭代法求根的程序

#include <stdio.h>
int main(){
double x0, x1;
double fx0, fx1;
double epsilon;
printf("请输入牛顿迭代法的初值x0:");
scanf("%lf", &x0);
printf("请输入牛顿迭代法的精度epsilon:");
scanf("%lf", &epsilon);
do {
fx0 = x0*x0*x0-x0-1;
fx1 = 3*x0*x0-1;
x1 = x0 - fx0/fx1;
x0 = x1;
}while(fabs(x1-x0) >= epsilon);
printf("最终求得的根为:%lf\n", x1);
return 0;
}