newton.cpp

牛顿法解决微分方程的实现代码

#include <math.h>
#include <iostream.h>
int main()
{
	float un = 8.0;//initalize
	float du;
	do 
	{
		du =( un * un * 2 + 5 * un - 4)/ (4 * un + 5);
		un = un - du;
		cout<<un<<endl;		
	} while(fabs(du) > 5e-5);
	return 0;
}