matlab 用牛顿法求解方程

定义函数
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 i10^(-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