Elliptic1D Left Dirichlet and Right Robin Boundary Conditions#
Solves the 1D Poisson equation with left Dirichlet and right Robin boundary conditions.
\[
-\nabla^2 u(x) = \pi^2 \sin(\pi x)
\]
with \(x\in[0,1]\).
The boundary conditions are given by
\[
a_nu + b_n\frac{du}{dx} = g_n
\]
with the left hand side boundary condition (Dirichlet) satisfying
\[
1u(0) + 0\frac{du(0)}{dx} = 10
\]
and the right hand boundary condition (Robin) satisfying
\[
400u(1) + 1\frac{du(1)}{dx} = 15
\]
This corresponds to the call to addScalarBC1D of addScalarBC1D(A,b,k,m,dx,dc,nc,v), where dc, nc, and vc are vectors which hold the coefficients for \(a\), \(b\), and \(g\) in the above system of equations. \(a=[1,400]\), \(b=[0,1]\) and \(g=[10,15]\).
Substituting these values in gives:
\[
u(0) = 10
\]
\[
400u(1) + \frac{du(1)}{dx} = 15
\]
The true solution is:
\[
u(x) = \sin(\pi x) + \frac{\pi - 3985}{401}x + 10
\]
This example is implemented in:
Additional MATLAB/ OCTAVE variants of this example with different boundary conditions: