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;