function x = newton(fname,dfname,x0,e,N)
% 用途：Newton迭代法解非線性方程f(x) = 0
% 格式：x = newton(fname,dfname,x0,e,N)  fname和dfname分別為表示f(x)
%      及其導函數的M函數句柄或內嵌函數，x0為迭代初值，e為精度要求（默認
%      1e-4），x返回數值解，并顯示計算過程設置迭代次數上限N以防止發散（默認500）
if nargin<5,N = 500;end
if nargin<4,e = 1e-4;end
x = x0;x0 = x+2*e;k = 0;
while abs(x0-x)>e&k<N,
    k = k+1;
    x0 = x;x = x0-feval(fname,x0)/feval(dfname,x0);
    disp(x)
end
if k == N,warning('已達迭代次數上限');end