3楼: Originally posted by langka1573 at 2015-12-11 08:10:46
原来学过计算方法这门课程，里面有详细内容，包括原理和代码例程。可以找书看看。

function f=fun(x) f=1000*exp(x)+(435/x)*(exp(x)-1)-1564; end

function y=df(x) y=1000*exp(x)+435*(exp(x)/x)-435*exp(x)/(x*x)+435/(x*x); end

clear all; clc; x0=1; r=1.0e-6; tic; i=1; x1=x0-fun(x0)/df(x0); while (abs(x1-x0)>r)&&(i<=60)     x0=x1;    x1=x0-fun(x0)/df(x0);   i=i+1; end   toc;  i  x1=vpa(x0-fun(x0)/df(x0),10) e=vpa(abs(x1-x0),10)  fm=vpa(fun(x1),10)