function [H]=limite(I,l)
H=I;
H(H<l)=0;
endfunction;

function [H]=inverse(I)
M=max(max(I));
H=I;
H=abs(H-M);
endfunction;
