function ObjVal = objharv(Chrom,sswitch);
% 目標函數的維數
   Dim = 20;
% 參數設置
   a = 1.1;
   x0 = 100;
   xend = x0;
   XENDWEIGHT = 0.4/(Dim^0.6);
   
% 種群的大小
   [Nind,Nvar] = size(Chrom);

% 種群檢查
   if Nind == 0
      %返回輸出圖形的標題文字
      if sswitch == 2
         ObjVal = ['HARVEST PROBLEM-' int2str(Dim)];
      %返回全局最小值
      elseif sswitch == 3
         ObjVal = -sqrt(x0*(a^Dim-1)^2/(a^(Dim-1)*(a-1)));
      %定義邊界矩陣大小和值
      else   
         % 變量的邊界        
         ObjVal1 = [0; 10*Dim];
         ObjVal = rep(ObjVal1,[1 Dim]);
      end
   elseif Nvar == Dim
      ObjVal = zeros(Nind,1);
      X = rep(x0,[Nind 1]);
      for irun = 1:Nvar,
         X = a*X - Chrom(:,irun);
      End
X;
ObjVal = -(sum(sqrt(Chrom)')' - XENDWEIGHT * abs(X-x0));
%否則，染色體格式錯誤
else
error('size of matrix Chrom is not correct for function evaluation');
   end