function ret=Aimfun(x)
% shubert function 

global sizepop  lenchrom  
 
bound=[-10 10;-10 10];   % 函數自變量的范圍

% 將binary數組轉化成十進制數組
n=1;
for i=1:length(lenchrom)
    data(i)=0;
    for j=lenchrom(i)-1:-1:0
        data(i)=data(i)+x(n).*2.^j;
        n=n+1;    
    end
end
data=bound(:,1)'+data./(2.^lenchrom-1).*(bound(:,2)-bound(:,1))';
        
% 計算適應度
temp1=0;
temp2=0;
for j=1:5
    temp1=temp1+j*cos((j+1)*data(1)+j);
    temp2=temp2+j*cos((j+1)*data(2)+j);
end
fitness=temp1*temp2;    % -10<=xi<=10, i=1,2;  

% 對適應度進行變換,使其越大越好
if fitness>=0
    fitness=1;
else fitness=1-0.05*fitness;
end
    
ret=fitness;