5
fespace Vh(Th,P3); // P1 FE space
6
Vh uh,vh; // unkown and test function.
7
func f=1; // right hand side function
8
func g=0; // boundary condition function
10
problem laplace(uh,vh) = // Definion of the problem
11
int2d(Th,qft=qf7pT)( dx(uh)*dx(vh) + dy(uh)*dy(vh) ) // bilinear form
12
- int2d(Th,qft=qf7pT)( f*vh ) // linear form
13
+ on(1,2,3,4,uh=g) ; // boundary condition form
15
laplace; // solve the problem plot(uh); // to see the result
16
plot(uh,ps="LaplaceP3.eps",value=true);