//境界は単位円
//$C=\{(x,y)|x=\cos t, y=\sin t, 0\le t \le 2\pi\}$
border C(t=0,2*pi){x=cos(t); y=sin(t);}
mesh Th = buildmesh (C(50));
plot(Th, wait=1); //三角形分割を見る
// ここからは hello.edp と同じ
plot(Th); //三角形分割を見る
fespace Vh(Th,P1); // 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);