first.cpp

newton迭代法通用程序（包括算例求解）

//该程序为Newton法解方程f(x)=0的通用程序
# include <iostream.h>
# include <math.h>
double F1 (double x)//定义原函数
{
	return ();//()括号内应写入原函数算法
}
double F2 (double x)//原函数一阶导
{
	return ();//()括号内应写入原函数一阶导算法
}
double Newton(double x0, double e)//Newton迭代法
{
	double x1;
	int i=0;//迭代截断参数
	do{
		i++;
		x1=x0;
		x0=x1-F1(x1)/F2(x1);
	} while (fabs(x0-x1)>e&&i<10000);
	if(i<10000) return x0;
	else return 555;//不收敛函数返回值（该值任取）
}
void main()
{
	double x0,x1;//x0为近似初值
	double e=10E-7;//容许误差
	cout<<"输入一个初值(实数范围)：";
	cin>>x0;
	cout<<"您输入的初值为："<<x0<<endl;
	x1=Newton(x0,e);
	if(x1==555) cout<<"该初值下Newton序列不收敛"<<endl;
	else cout<<"初值为 "<<x0<<" 时，Newton序列收敛于"<<x1<<endl;
}