cost_knapsack(X , p , W , c , alpha)

function S = cost_knapsack(X , p , W , c , alpha)

[n , N] = size(X);

ON      = ones(1 , N);

S       = alpha*sum(min((c(: , ON) - W*X) , 0))  + sum(p(: , ON).*X);