fem_neumann

fem_neumann, an Octave code which sets up a time-dependent reaction-diffusion equation in 1D, with Neumann boundary conditions, discretized using the finite element method, by Eugene Cliff.

Licensing:

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

Languages:

fem_neumann is available in a MATLAB version and an Octave version.

Related Data and Programs:

fem_neumann_test

fem1d, an Octave code which applies the finite element method to a 1D linear two point boundary value problem.

fem1d_adaptive, an Octave code which applies the finite element method to a 1D linear two point boundary value problem using adaptive refinement to improve the solution.

fem1d_bvp_linear, an Octave code which applies the finite element method, with piecewise linear elements, to a two point boundary value problem in one spatial dimension.

fem1d_nonlinear, an Octave code which applies the finite element method to a 1D nonlinear two point boundary value problem.

fem1d_pack, an Octave code which contains utilities for 1D finite element calculations.

fem1d_pmethod, an Octave code which applies the p-method version of the finite element method to a 1D linear two point boundary value problem.

Author:

Eugene Cliff

Reference:

1. Jeffrey Borggaard, John Burkardt, John Burns, Eugene Cliff,
   Working Notes on a Reaction Diffusion Model: a Finite Element Formulation.

Source Code:

• basic_hat.m, evaluates the basic hat function on [-1, 1].
• plot_rd.m, makes two plots of the reaction-diffusion solution.
• rd_lin_spline.m, discretizes and solves a 1D reaction diffusion problem.