function x=Newton(fname,dfname,x0,e,N)
%用途：Neton迭代法解非線形方程f(x)=0;
%fname和dfname分別表示f(x)及其導函數的M函數句柄或內嵌函數表達式
%x0為迭代初值，e為精度
%x為返回數值解，并顯示計算結果，設置迭代次數的上限
if nargin<5,N=500;end
if nargin<4,e=1e-4;end
x=x0;x0=x+2*e;k=0;
fprintf('It.no=%2d      x[%2d]=%12.9f\n',k,k,x)
while abs(x0-x)>e&k<N,
    k=k+1;
    x0=x;x=x0-feval(fname,x0)/feval(dfname,x0);
    fprintf('It.no=%2d     x[%2d]=%12.9f\n',k,k,x)
end
if k==N,fprintf('已達到迭代次數上限');end