#include <iostream>
#include <cmath>
using namespace std;

int main()
{
	double X0,X;
	double eps;
	double FX,dFX;

	cout<<"請輸入X0的值: ";
	cin>>X0;

	cout<<"請輸入控制精度：";
	cin>>eps;

	X=X0;

	do
	{	
		FX=pow(X,3)-3*X-1;
		dFX=3*pow(X,2)-3;
		X0=X;
		X=X-FX/dFX;
	}
	while(fabs(X-X0)>=eps);

	cout<<"方程F(X)的近似解:X="<<X<<endl;

	return 0;
}