function [atil b]=contorno(nv,a,n,m)
atil=a;
b=zeros(nv,1);
for i=1:n
  for j=1:m
    if ((i==1 || i==n) || (j==1 || j==n))
          k=(j-1)*n + i;
          atil(k,:)=0;
          atil(k,k)=1;
          b(k)=0;
          if (i==n) b(k)=1.; end
    end
  end
end