牛顿迭代法的matlab代码是怎样的?
1个回答
展开全部
主程序:
function [k,x,wuca,yx] = newton(x0,tol)
k=1;
yx1=fun(x0);
yx2=fun1(x0);
x1=x0-yx1/yx2;
while abs(x1-x0)>tol
x0=x1;
yx1=fun(x0);
yx2=fun1(x0);
k=k+1;
x1=x1-yx1/yx2;
end
k;
x=x1;
wuca=abs(x1-x0)/2;
yx=fun(x);
end
分基配程序1:
function y1=fun(x)
y1=sqrt(x^2+1)-tan(x);
end
分程序2:
function y2=fun1(x)
%函数fun(x)的导数
y2=x/(sqrt(x^2+1))-1/((cos(x))^2);
end
结果:
[k,x,wuca,yx] = newton(-1.2,10^-5)
k =8
x =0.9415
wuca =4.5712e-08
yx =-3.1530e-14
[k,x,wuca,yx] = newton(2.0,10^-5)
k 掘大=243
x =NaN
wuca =NaN
yx 搏散指=NaN
本回答被网友采纳
已赞过
已踩过<
评论
收起
你对这个回答的评价是?
推荐律师服务:
若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询