newton.cpp

牛顿法——数值分析,程序代码

// diedai.cpp : Defines the entry point for the console application.
//
//第104页24题：用牛顿法求方程的根
#include "stdafx.h"
#include "iostream.h"
#include "math.h" 

float fun(float x)
{
	float m;
	m=x-(x*x*x+2*x*x+10*x-20)/(2*x*x+4*x+10);
	return m;
}

int main(int argc, char* argv[])
{
	int k=0,N;
	float x0,x1;
	x0=1.0;
	cout<<"求方程x*x*x+2*x*x+10*x-20=0的根"<<endl;
	cout<<"输入需要迭代的次数：";
	cin>>N;
	for(;k<=N;k++)
	{
		x1=fun(x0);	
		if(fabs(x1-x0)<0.000001)
			break;
		x0=x1;
	}
	cout<<"求出的根为x="<<x1<<endl;
	return 0;
}