function [sol,eval] = coranaMin(sol,options)
numVar=size(sol,2)-1;
eval = corana(sol(1:numVar));
eval = -eval;
%設置變量的邊界

bounds = ones(4,1)*[-10000 10000];
%遺傳算法優(yōu)化
[p,endPop,bestSols,trace]=ga(bounds,'coranaMin');

%性能跟蹤
plot(trace(:,1),trace(:,3),'y-')
hold on
plot(trace(:,1),trace(:,2),'r-')
xlabel('Generation');
ylabel('Fittness');
legend('解的變化','種群平均值的變化');