java视频教程马士兵:一道c语言题

程序流程分析：
① 赋值x0＝1.5,即迭代初值；
② 用初值x0代入方程中计算此时的f(x0)及f’(x0)，程序中用变量f描述方程的值，用fd描述方程求导之后的值；
③ 计算增量d=f/fd；
④ 计算下一个x，x=x0-d;
⑤ 把新产生的x替换x0，为下一次迭代做好准备;
⑥ 若d绝对值大于1e-3，则重复②③④⑤步。

源程序代码：

#include <math.h>
main()
{
float x,x0,d,f,fd;
x0=0;
do {
f=2*x0*x0*x0-4*x0*x0+3*x0-6;
fd=6*x0*x0-8*x0+3;
d=f/fd;
x=x0-d;
x0=x;
}while(fabs(d)>1e-3);
printf("x=%f\n",x);
}

这是一个具体的例子

/*----------用牛顿迭代法求解方程的根-----------*/
#include<stdio.h>
#include<math.h>
#include<stdlib.h>
#define LOWLIMIT 1.E-6
#define MAX 100
#define infinity 1.E300
double f(double x)
{
return x*x+3*x-3;/*依你的需要自定义修改函数*/
}
double diff(double x)
{
return 2*x+3; /*定义函数的导数*/
}
finite(double x)
{
if(x > infinity)return 0; /*此函数的功能是防止上溢*/
else return 1;
}
/*-------牛顿法迭代算法--------*/
void newton()
{
int k;
double x;
double delta_x;
printf("Enter x0,the initial guess:");
scanf("%lf",&x);
for(k=0;k<MAX;k++){
delta_x=-f(x)/diff(x);
if(!finite(delta_x))break;
x+=delta_x;
if(fabs(delta_x) < LOWLIMIT)break;
}
/*分析结果*/
if(k==MAX)
printf("\t No convergence in %d iterations.\n",MAX); /*与X轴无交点*/
else if(!finite(delta_x))
printf("\tOverflow on delta_x.\n");/*溢出*/
else
printf("\tProcess converged.\n"); /*与X轴有交点，成功求出根*/
printf("\tx=%e,delta_x=%e.\n",x,delta_x);
printf("f(x)=%e,f'(x)=%e.\n",f(x),diff(x));

}

void main()
{
int do_again;
puts("\n\n-------Solve f(x)=x^2+3x-3 by Newton's method.--------");
do{newton();
printf("\n\nEnter 1 to continue ,0 to exit:");
scanf("%d",&do_again);
}while(do_again);
}

建议看看计算方法一书。