牛顿迭代法c++代码

数值分析实验报告(二) 一.实验名称:

牛顿法求非线性方程的根。

二.实验目的:

运用牛顿法求方程根

三.题目:

用牛顿法求解人口方程

156.4=100e l+43.5

l

(e l-1)

要求误差小于4

10-.

四.程序:

#include

#include

double f(double x)

{

return 100*exp(x)+43.5/x*(exp(x)-1)-156.4;

}

double g(double y)

{

return 100*exp(y)-43.5/(y*y)*(exp(y)-1)+43.5/y*exp(y);

}

double newton(x0,eps)

double x0,eps;

{

double x;

while(1)

{

x=x0-f(x0)/g(x0);

if(fabs(x-x0)

else x0=x;

}

}

main()

{

clrscr();

printf("the root is %10.6f",newton(0.2,1e-4));

getchar();

}

五.运行结果:

The root is 0.1010

六.体会:

牛顿法在单根附近有较快的收敛速度.通过这次学习,我又学会了一种方法,在以后会更好的掌握和应用.

信计0 5 0 2 13 孟松见