newton.cpp

牛顿 拉格朗日 分段线性

#include <iostream.h>
#include "math.h"
double f(double x)
{
	return x*exp(x)-1;
}
double ff(double x)
{
	return x*exp(x)+exp(x);
}
void main()
{
	double x0,x1,c; 
	int n,k;
	cin>>x0>>c>>n;
	int flag=0;
	for(k=1;k<=n;k++)
	{
		if(ff(x0)==0)
		{
			cout<<"error!"<<endl;
			break;
		}
		else
		{
		    x1=x0-f(x0)/ff(x0);
            	
		    if(fabs(x1-x0)<c)
			{
			    cout<<"x1="<<x1<<endl;
			    flag=1;
			    break;

			}
			x0=x1;
		}
	}	
	if(k==n+1&&flag==0) cout<<"fault!"<<endl;
}