mesh Th=square(10,10); // (細かく)領域$\Omega=]0,1[^2$の三角形分割
plot(Th); //三角形分割を見る
fespace Vh(Th,P2); // 2次要素有限要素空間
Vh u,v; // $u,v\in V_h(\Omega)$

/*弱形式を記述し Poisson と名づける。
$$\int_{\Omega}(\partial_{x}u\partial_{x}v+\partial_{y}u\partial_{y}v) 
dxdy -\int_{\Omega}1\cdot v dxdy = 0$$ */
problem Poisson(u,v) =
    int2d(Th)(  dx(u)*dx(v) + dy(u)*dy(v)) 
  - int2d(Th) ( 1*v )
  + on(1,2,3,4,u=0)  ;

Poisson; // Poisson を解く
//$u$ の等高線を描き、Postscriptファイルで保存
plot(u,ps="hello01.eps");