function f = fitness(x,PP,b,punish)
    f = 0;
    for p = 1:size(PP,2)
        f = sum(((PP(:,p) * x(p) - b) + abs(PP(:,p) * x(p) - b))*punish)+f;
    end
end