Nous allons définir la fonction f et construire une triangulation du cercle unité avec en écrivant dans un fichier de nom (par exemple) "trace.pde" le texte suivant:
/* fichier trace.pde */
n:=50;
border(1,0,2*pi,n) begin
x := cos(t);
y := sin(t);
end;
buildmesh(1000);
savemesh('trace.msh');
f = x * y;
plot(f);
save('trace.dta',f);
Ce texte est écrit dans la syntaxe Gfem. Il définit une frontière
qui porte le numéro 1, dont le paramètre t dans la description paramétrique
varie de 0 à 
et qui sera découpé en segments par n sommets.
Le mailleur automatique s'arrêtera si le nombre de sommets générés par
l'algorithme Delaunay-Voronoi dépasse 1000
(cette incongruité facilite l'optimisation
de la mémoire allouée à l'intérieru de freefem).
Enfin le résultat sera stocké
dans le fichier "trace.msh" selon un format qui sera détaillé dans le paragraphe
suivant.
Le fichier des valeurs de f aux sommets de la triangulation s'appelle "trace.dta". Le format est
- Le nombre de sommets
- les valeurs de f sur chaque sommets sur une ligne par sommet
Certaine version de freefem ayant prévue le cas f a valeur complexe, les
valeurs de f aux sommets sont 2 nombes, la partie reel et la partie imaginaire.
Verifiez quel version vous avez.
