fem1d_lagrange

fem1d_lagrange, a MATLAB code which sets up the matrices and vectors associated with the finite element method (FEM) solution of a boundary value problem (BVP) -u''+u=f(x), using Lagrange basis polynomials.

Licensing:

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

Languages:

fem1d_lagrange is available in a C version and a C++ version and a FORTRAN90 version and a MATLAB version.

Related Programs:

fem1d, a MATLAB code which applies the finite element method (FEM) to a 1D linear two point boundary value problem (BVP).

fem1d_lagrange_test

fem1d_spectral_numeric, a MATLAB code which applies the spectral finite element method (FEM) to solve the two point boundary value problem (BVP_ u'' = - pi^2 sin(x) over [-1,+1] with zero boundary conditions, using as basis elements the functions x^n*(x-1)*(x+1), and carrying out the integration numerically, using MATLAB's quad() function, by Miro Stoyanov.

lagrange_interp_1d, a MATLAB code which defines and evaluates the Lagrange polynomial p(x) which interpolates a set of data, so that p(x(i)) = y(i).

Source Code:

• fem1d_lagrange_stiffness.m computes the stiffness matrix, mass matrix, and right hand side for a finite element problem that uses Lagrange basis polynomials.
• lagrange_derivative.m, evaluates the derivative of Lagrange basis polynomials.
• lagrange_value.m, evaluates the Lagrange basis polynomials.
• legendre_set.m, sets the points and weights of a Legendre quadrature rule.
• r8vec_print.m, prints an R8VEC.