牛頓迭代法matlab實現如下:
function [x_star,index,it] = Newton(fun,x,ep,it_max)
%求解非線性方程的牛頓法
%第一個分量是函式值,第二個分量是導數值
% x為初始點
% ep為精度,當 | x(k)-x(k-1) |=it_max)
% it為迭代次數
if nargin