Dear Bhupesh Dutt Sharma,

VAM in Matlab could be as follows:

function [ Solution, OverallCost ] = VAM(CostsMtx, resources_col, demands_row)

C_start = CostsMtx;

C = C_start;

m = size(C,1);

n = size(C,2);

a = resources_col;

b = demands_row;

X = zeros(m,n);

stop = 0;

while stop == 0

for i = 1:m

for j = 1:n

if a(i,1) == 0

C(i,j) = max(C(:,j));

end

if b(1,j) == 0

C(i,j) = max(C(i,:));

end

end

end

C_sort_col = sort(C,1);

C_sort_row = sort(C,2);

Diff_customer = abs(C_sort_col(1,:) - C_sort_col(2,:)); 

Diff_supplier = abs(C_sort_row(:,1) - C_sort_row(:,2)); 

for i = 1:m

if a(i,1) == 0

Diff_supplier(i,1) = 0;

end

end

for j = 1:n

if b(1,j) == 0

Diff_customer(1,j) = 0;

end

end

Max_Diff_customer = max(Diff_customer);

Max_Diff_supplier = max(Diff_supplier);

Customer_nr = find(Diff_customer==max(Max_Diff_customer,Max_Diff_supplier)); 

Supplier_nr = find(Diff_supplier==max(Max_Diff_customer,Max_Diff_supplier)); 

if isempty(Customer_nr) == 0

Supplier_nr_ = find(C(:,Customer_nr(1)) == min(C(:,Customer_nr(1))));

X(Supplier_nr_(1),Customer_nr(1)) = min(a(Supplier_nr_(1),1),b(1,Customer_nr(1)));

a(Supplier_nr_(1),1) = a(Supplier_nr_(1),1) - X(Supplier_nr_(1),Customer_nr(1));

b(1,Customer_nr(1)) = b(1,Customer_nr(1)) - X(Supplier_nr_(1),Customer_nr(1));

Supplier_nr = [];

end

if isempty(Supplier_nr) == 0

Customer_nr_ = find(C(Supplier_nr(1),:) == min(C(Supplier_nr(1),:)));

X(Supplier_nr(1),Customer_nr_(1)) = min(a(Supplier_nr(1),1),b(1,Customer_nr_(1)));

a(Supplier_nr(1),1) = a(Supplier_nr(1),1) - X(Supplier_nr(1),Customer_nr_(1));

b(1,Customer_nr_(1)) = b(1,Customer_nr_(1)) - X(Supplier_nr(1),Customer_nr_(1));

end

%Stop condition:

a1 = a > 0;

b1 = b > 0;

if sum(a1) == 1

stop = 1;

for j = 1:n

if b(j) > 0;

X(a1 == 1,j) = b(j);

end

end

end

if sum(b1) == 1

stop = 1;

for i = 1:m

if a(i) > 0;

X(i,b1 == 1) = a(i);

end

end

end

end

Solution = X;

OverallCost = sum(sum(C_start .* X));

end

Please, let me know if there are any problems with the code.

Regards,

Remigiusz Iwańkowicz

Hi Remigiusz,

excellent, works like a charm!

Kine regards

Mario Galić

Hello!!!

could you please provide code for Russell's Approximation Method.

thanks

Sorry, I haven't coded this method yet.