xgfc.cpp

来自「在数值计算中利用弦割发解方程组的c++程序」· C++ 代码 · 共 36 行

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

int main()
{
	double X0,X1,X2;
	double eps;
	double f0,f1,f2;

	cout<<"请输入X0,X1的值：";
	cin>>X0>>X1;

	cout<<"请输入控制精度：";
	cin>>eps;

	do
	{
		f0=pow(X0,3)-3*X0-1;
		f1=pow(X1,3)-3*X1-1;

		X2=X1-f1*(X1-X0)/(f1-f0);
		X0=X1;
		X1=X2;
	}
    while(fabs(X1-X0)>=eps);
	
	cout<<"方程F(X)的近似解:X="<<X2<<endl;

	return 0;
}