int n=20,nn=n+10;
real[int] xx(nn),yy(nn);
// build square $]-1,1[^2$
mesh Th=square(10,10,[2*x-1,2*y-1]); 
fespace Vh(Th,P1);
Vh u=2-x*x-y*y;
// old code ::
if(1)
{
   savemesh(Th,"mm",[x,y,u*.5]); //  save mm.points and mm.faces file for medit
   // build a mm.bb file 
	{ ofstream file("mm.bb"); 
	file << "2 1 1 "<< u[].n << " 2 \n";
	int j;
	for (j=0;j<u[].n ; j++)  
	  file << u[][j] << endl; 
    }  
    // call ffmedit command 	
    if(!NoGraphicWindow) exec("ffmedit mm");
    // clean files
    exec("rm mm.bb mm.faces mm.points");
}
else 
{
// new code 
load "medit"
    meshS Th3= movemesh23(Th,transfo=[x,y,u*0.5]);
    medit("mm",Th3);// bug un color of u ... FH 
}

// FFCS: testing 3D plots
plot(u);
