%Homotopy_Optimization_Method
clear;
global x00 lambda;
m=20;deta=1/m;
format long;
%%%%%%%%%%%%%%%%%%
%%%%%%%%%%(1)
n=4;
for i=1:n
    x1(i,1)=rand(1);
end
lambda0(1)=0;x00=x1;

xtemp1=x1;
for i=2:m+1
    lambda0(i)=lambda0(i-1)+deta;lambda=lambda0(i);
    
    x1=mulNewton_H(x1);BB='牛頓法(基本方法)';
    %x1=mulDNewton_H(x00);BB='牛頓下山法';
    % x1=mulGSND(x0);BB='高斯牛頓法';
   %  x1=mulDamp(x0);BB='阻尼最小二乘法';
    %x1=mulConj(x0);BB='共軛梯度法???????????????????';
    %x1=mulFastDown(x0);BB='最速下降法??????????';
    %x1=mulNumYT(x0);BB='數值延拓法!???/????????';
   %x1=newtoniteration6(x0);BB='牛頓迭代法改進《機床與液壓》';
    x1=x1;f1=myfun(x1);
  end
%%%%%%%%%%
format short;
x1',f1'