%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 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 F=sec_mem(G,d,n)

if d==1
   h=1/(n+1);
   F=ones(n,1)+G/h^2;
end

if d==2
   h=1/(n+1);
   F=ones(n*n,1)+G/h^2;
end
