newton下山法

课本27页例2———8

#include  
#include  
using namespace std;  
int main()  
{  
    cout<<"请输入初始值X："<>x0;  
    F0=x0*x0*x0-1-x0;  
    F1=3*x0*x0-1;  
      
    int i,j;  
	double  t=1.0;
    for( i=0;i<50;i++)  
    {  
        if(F1==0)  
            break; 
		for(t=1;;)
       {
		   x1=x0-1.0/t*(F0/F1);  
		   
		   if(abs(x1*x1*x1-x1-1)