function x=newton(f,df,x0,e,N)
% Newton法解方程f(x)=0根的通用程序
% 格式x=newton(f,df,x0,e,N) f表示f(x),df表示f(x)的導函數,x0表示迭代初值,e為容許誤差,N為迭代次數的上限
x=0;
x=x0; %k=0;
k=0;
disp('迭代結果:')
for k=1:N
    x0=x;
    x=x0-feval(f,x0)/feval(df,x0);
    k=k+1;
    disp(x)
    if abs(x-x0)<e||abs(feval(f,x0))<e break;end
end
if k==N error('已超過迭代上限'); end
disp('迭代次數'),disp(k)