第1个回答 推荐于2017-09-30
定义函数
function y=nd(x)
y= 0.036-((x/2090.7).^(1/0.1585))-x/182000
function y=nd0(x)
y=- (1/0.1585)*(x/2090.7).^(1/0.1585-1)*(1/2090.7) - 1/182000
主程序
x=0.5; %迭代初值
i=0; %迭代次数计数
while i<=100;
y=x-nd(x)/nd0(x); %牛顿迭代格式
if abs(y-x)>10^(-5); %收敛判断
x=y;
else break
end
i=i+1;
end
fprintf('\n%s%.4f \t%s%d','x=',x,'i=',i) %输出结果
x=-100.5000 i=101本回答被提问者采纳