//2次要素の方が良い近似になる。
mesh Th1=square(10,10); //領域$\Omega=]0,1[^2$の三角形分割
mesh Th2=square(20,20); //細かく
fespace Vh2(Th1,P2); // 粗いメッシュで2次要素
fespace Vh1(Th2,P1); // 細かいメッシュで1次要素
Vh1 u1,v1; // $u_1,v_1\in V_{h}1(\Omega)$
Vh2 u2,v2; // $u_2,v_2\in V_{h}2(\Omega)$

// 1次要素
solve Poisson1(u1,v1) =
    int2d(Th2)(  dx(u1)*dx(v1) + dy(u1)*dy(v1)) 
  - int2d(Th2) ( 1*v1 )
  + on(1,2,3,4,u1=0)  ;
// 2次要素
solve Poisson2(u2,v2) =
    int2d(Th1)(  dx(u2)*dx(v2) + dy(u2)*dy(v2)) 
  - int2d(Th1) ( 1*v2 )
  + on(1,2,3,4,u2=0)  ;

//$u1$ と $u2$ を重ねる
plot(u1,u2);