C语言程序设计 牛顿迭代法

/* 牛顿迭代法求解方程 在2.0附近的一个根 */

#include<stdio.h>
#include<math.h>
void main()
{
  double x1,x,f,f1;
  scanf("%lf",&x1);
  do
  {
    x=x1;
    f=x*x*x*x-3*pow(x,3)+1.5*x*x-4;
    f1=4*pow(x,3)-9*x*x+3*x;
    x1=x-f/f1;
  }while(fabs(x1-x)>1e-6);
  printf("x=%lf\n",x);
  x=2.648937;
  printf("%lf",x*x*x*x-3*pow(x,3)+1.5*x*x-4);
}