///////////////////////////////// // Resolution de l'equation // // de la chaleur dans un carre // ///////////////////////////////// //Definition du maillage mesh Th; //Nombre de mailles suivant x et y int Nbnoeuds=10; Th=square(Nbnoeuds,Nbnoeuds,[x,y]); //Fonction de x et de y func f=1.; //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)) -int2d(Th)(f*vh) +on(1,2,3,4,uh=0) ; //Resolution du probleme chaleur; //On affiche le resultat plot(uh,wait=1);