SetCask.m

for i =1:COCE.CaskArea
    COCE.Data(COCE.ExternalArea(1)+1:COCE.ExternalArea(1)+ COCE.InternalArea(1),...
      COCE.PopulationSize(2)/2-COCE.CaskArea+i) = COCE.InternalValue - COCE.CaskContrast*(i*i)/COCE.CaskArea^2;
end
for i =1:COCE.CaskArea
   COCE.Data(COCE.ExternalArea(1)+1:COCE.ExternalArea(1)+COCE.InternalArea(1),...
      COCE.PopulationSize(2)/2+i) = ...
         COCE.InternalValue + COCE.CaskContrast*((COCE.CaskArea+1-i)^2)/COCE.CaskArea^2;

%    x(16:48,128+i) = .5 + val*((n+1-i)^2)/n^2;
end
%{
for i =1:COCE.CaskArea
    x(16:48,128-n+i) = .5- val*(i*i)/n^2;
end
for i =1:n
    x(16:48,128+i) = .5 + val*((n+1-i)^2)/n^2;
end
%}