1mesh Th = square(10,10); 2fespace Wh(Th,[P2,P2]); // R2->R2 3Wh [F1,F2]; //def de l'inconnu et des fonct tests 4 5F1 = 1; 6// il y a un bug ici car dans un element fini vectoriel 7// il est impossible defini juste une composante 8// il faut faire: 9 10cout << F1(0.5, 0.6) << endl; 11cout << F2(0.5, 0.6) << endl; 12cout << F1(0.5, 0.6); 13[F1,F2]=[1,2]; 14cout << " -- F1 " << F1(0.5, 0.6) << endl; 15cout << " -- F2 " << F2(0.5, 0.6) << endl;