Matlab - Función MILP!!

% Definir la matriz de costos (beneficios)

costos = [3, 1, 5; 2, 6, 4; 7, 1, 8];

% Número de tareas y trabajadores

num_tareas = size(costos, 1);

num_trabajadores = size(costos, 2);

% Definir la función objetivo (maximizar la productividad)

f = -reshape(costos, 1, []);

% Restricciones de asignación (una tarea por trabajador y viceversa)

Aeq = zeros(num_tareas + num_trabajadores, num_tareas * num_trabajadores);

beq = ones(num_tareas + num_trabajadores, 1);

for i = 1:num_tareas

    Aeq(i, (i-1)*num_trabajadores + 1:i*num_trabajadores) = 1;

end

for j = 1:num_trabajadores

    Aeq(num_tareas + j, j:num_trabajadores:num_tareas*num_trabajadores) = 1;

end

% Variables binarias (asignación o no asignación)

intcon = 1:num_tareas * num_trabajadores;

% Resolver el MILP

[x, fval, exitflag] = intlinprog(f, intcon, [], [], Aeq, beq, zeros(num_tareas * num_trabajadores, 1), ones(num_tareas * num_trabajadores, 1));

% Reshape de la solución para obtener la asignación

asignacion = reshape(x, num_tareas, num_trabajadores);

% Mostrar los resultados

disp('Asignación óptima:');

disp(asignacion);

disp(['Productividad total: ' num2str(-fval)]);