求方程根的近似值,
分别用matlab和c编译
例 exp(x)+10x-2
matlab:建立m文件(牛顿迭代)
function x=newton(fname,dfname,x0,e)
%fname为函数名dfname为导函数x0为迭代初值N为最大迭次数
N=100;
x=x0;
x0=x+2*e;
k=0;
while abs(x0-x)>e&k<N
k=k+1
x0=x;
x=x0-feval(fname,x0)/feval(dfname,x0);
disp(x)
end
if k==N warning(‘error’);
end
命令窗口输入:
fun=inline(exp(x)+10*x-2);
dfun=inline(exp(x)+10);
x=newton(fun,dfun,1.5,1e-6)
用c语言编译如下
#include<stdio.h>
#include<math.h>
float(float x){
float y;
y=x-exp(-x);
return y;
}
float df(float x){
float y;
y=1+exp(-x);
return y;
}
main(){
int N,k;
float x0,eps,x1;
printf("input x0,eps,N:");
scanf("%f,%f,%d",&x0,&eps,&N);
for(k=1;k<=N;k++){
if(x0=0){printf("\n结果为%f\n",&x0);break;}
x1=x0-f(x0)/df(x0);
printf("\nx=%d=%f\n",k,x1);
if(fabs(x1-x0)<eps){printf("\n结果为%f\n k=%d",x1,k);break;}
x0=x1;
}
if(k>N)printf("error!!!");
}