mesh Th=square(10,10); fespace Vh(Th,P1); // P1 FE space Vh uh,vh; // Unkown and test function func f=1; // Right hand side function func g=0; // Boundary condition function // Definition of the problem problem laplace(uh,vh,solver=GMRES,tgv=1e5) = int2d(Th)( dx(uh)*dx(vh) + dy(uh)*dy(vh) ) // bilinear form - int2d(Th)( f*vh ) // linear form + on(1,2,3,4,uh=g) ; // boundary condition form laplace; // Solve the problem plot(uh,value=true); // see the result Save Run Reset

canvas not supported