function x=nanewton(fname,dfname,x0,e,N)
%用途：牛頓法解非線性方程f(x)=0
%格式：x=nanewton(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