function [atil]=contorno(a,n,m)
atil=a;
for i=1:n
  for j=1:m
    if ((i==1 || i==n) || (j==1 || j==n) || (i<n/2 && j<m/3))
          k=(j-1)*n + i;
          atil(k,:)=0;
    end
  end
end