Matlab utiliza el método de Newton para resolver ecuaciones

función y=nd(x)

y= 0.036-((x/2090.7).^(1/0.1585))-x/182000

función y=nd0(x)

y=- (1/0.1585)*(x/2090.7).^(1/0.1585-1)*(1/2090.7) - 1/182000

Programa principal

x=0.5; valor inicial de iteración

i=0; número de iteraciones

mientras ilt;=100;

p>

y=x-nd(x)/nd0(x); formato de iteración de Newton

if abs(y-x)gt; juicio

x=y

else break

end

i=i 1

end<; /p>

fprintf('\ns.4f \tsd', 'x =', x, 'i=', i) resultado de salida

x=-100.5000 i=101