/////////////////////////////////
// Resolution de l'equation    //
// de la chaleur dans un carre //
/////////////////////////////////

//Definition du maillage
mesh Th;

//Nombre de mailles suivant x et y
int Nbnoeuds=50;

Th=square(Nbnoeuds,Nbnoeuds,[x,y]);



//Fonction de x et de y
func f=x*y;

//Definition de l'espace des elements finis P1 associe
//au maillage Th
fespace Vh(Th,P1);

//uh et vh sont des elements de Vh
Vh uh,vh;


//Definition du probleme variationnel
problem chaleur(uh,vh,solver=LU)=
	int2d(Th)(dx(uh)*dx(vh)+dy(uh)*dy(vh)+ uh*vh)
	-int2d(Th)(f*vh)
;


//Resolution du probleme
chaleur;

//On affiche le resultat
plot(uh,wait=1);