%牛頓法極小化函數的通用子函數zuisu.m
%輸入變量為初始的迭代點，輸出變量為極小值點
function x0=newton(x)
 %判斷梯度范數是否滿足計算精度的要求.是,標志變量設為1，輸出結果;否,標志變量設為0，繼續迭代
if sum(abs(tidu(x)).^2)<1e-4
    flag=1;
    x0=x;
else 
    flag=0;
end
%循環求解函數的極小點
while flag==0
      d=-inv(hesse(x))*tidu(x)';
      a=gold(x,d);
      x=x+a*d';
      %判斷梯度范數是否滿足計算精度的要求.是,標志變量設為1，輸出結果;否,標志變量設為0，繼續迭代
      if sum(abs(tidu(x)).^2)<1e-4
         flag=1;
         x0=x;
     else 
         flag=0;
     end
 end