%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% MAO 2002/2003
%
%   =================================================
%   Resolution numérique d'EDP par différences finies
%   =================================================
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%   Resolution d'un probleme de Poisson
%   sur un maillage régulier composé de n points
%   internes ; les conditions aux limites sont supposées
%   de type Dirichlet
%
%   - en dimension 1 si d=1
%   - en dimension 2 si d=2
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%==============================
% Données du problème
%==============================
d=2; %dimension du problème
n=15;%nombre de points de discrétisation 
     %par direction du domaine
     
%-----------------------------------------
%Choix des conditions aux limites
% lab_bord=1 : Dirichlet homogène
% lab_bord=2 : x (si d=1)  x*y (si d=2)
%-----------------------------------------
lab_bord=2;
%==============================
% Construction de la matrice
%==============================

A=matlapla(d,n);

%===========================================
% Définition des conditions aux limites
%===========================================

G=cond_lim(d,n,lab_bord);

%=================================
% Construction du second membre F
% 
%=================================

F=sec_mem(G,d,n);

% Prise en compte des conditions aux limites

%==============================
% Résolution du système
%==============================

X=A\F;

%==============================
% Représentation de la solution
%==============================

 repres_sol