FEM2D_POISSON_LAKE
A 2D Poisson Problem in a Lake

FEM2D_POISSON_LAKE is a MATLAB library which sets up the data and geometry for a 2D Poisson problem in a lake, suitable for solution by fem2d_poisson.

Licensing:

The computer code and data files described and made available on this web page are distributed under the GNU LGPL license.

Related Data and Programs:

FEM2D_POISSON, a MATLAB program which solves Poisson's equation on a triangulated region, using the finite element method. In order to run, it requires user-supplied routines that define problem data.

Source Code:

• dirichlet_condition.m, defines the Dirichlet boundary condition values.
• h_coef.m, defines the coefficient function H(X).
• k_coef.m, defines the coefficient function K(X).
• rhs.m, defines the right hand side function.
• exact_solution.m, defines the exact solution of the problem.

Output files include:

• lake_output.txt, output from a run of the program;
• lake_nodes.txt, a text file containing a list, for each node, of its X and Y coordinates;
• lake_nodes.png, a PNG image of the nodes;
• lake_elements.txt, a text file containing a list, for each element, of the three nodes that compose it;
• lake_elements.png, a PNG image of the element mesh;
• lake_values.txt, a text file containing the solution U at each node (X,Y);
• lake_values.png, a contour plot created by triangulation_order3_contour.

You can go up one level to the MATLAB source code page.