牛顿迭代.cpp

数值计算的源代码

#include<stdio.h>
#include<math.h>
float y(float x0)
{return pow(x0,3)-x0-1;
}


float yd(float x0)
{return 3*pow(x0,2)-1;
}

float x(float x0)
{return x0-y(x0)/yd(x0);
}

void main()
{float x0=1.5,x1;
printf("求x(3)-x-1=0在[1,2]的根,误差10(-3).\n");
x1=x(x0);
for(;1;)
{if(fabs(x1-x0)>=1e-3)
{x0=x1;
x1=x(x0);}
else break;

}

printf("%f",x1);


}