function C=maxJ_L(A,B);

%
% Produto max-J_L
%

[m,n]=size(A);

if size(B,1)~=n
	disp('Inner matrix dimensions must agree.');
	return
end

C=zeros(m,size(B,2));
for i=1:m
	for j=1:n
		C(i,:)=max(C(i,:),A(i,j)-B(j,:));
	end
end