%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% MAO 2002/2003
%
%   =================================================
%   Resolution numérique d'EDP par différences finies
%   =================================================
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%   Prise en compte des 
%   conditions aux limites de type Dirichlet
%
%   - en dimension 1 si d=1 (fonction bord1)
%   - en dimension 2 si d=2 (fonction bord2)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function G=cond_lim(d,n,lab_bord)

if d==1
   G=zeros(n,1);
   G(1)=bord1(0,lab_bord);
   G(n)=bord1(1,lab_bord);
end

if d==2
   h=1/(n+1);
   G=zeros(n*n,1);
   G(1)=bord2(h,0,lab_bord)+bord2(0,h,lab_bord);
   for i=2:n
      G(i)=bord2(i*h,0,lab_bord);
   end
   G(n)=bord2(n*h,0,lab_bord)+bord2(0,n*h,lab_bord);

   for j=2:n-1
      G(1+(j-1)*n)=bord2(0,(j-1)*h,lab_bord);
      G(n+(j-1)*n)=bord2(1,(j-1)*h,lab_bord);
   end
   G(1+n*(n-1))=bord2(h,1,lab_bord)+...
      bord2(0,1-h,lab_bord);
   
   for i=2:n-1
      G(i+n*n-n)=bord2(i*h,1,lab_bord);
   end
   G(n*n)=bord2(1,1-h,lab_bord)+...
      bord2(1-h,1,lab_bord);

   
end
