We will consider the system of equations obtained by discretizing the two-point boundary value problem:
-y'' + y' + y = 0,
y'(0) = 0,
-y(1) + y'(1) = 1.
We will use second-order finite difference approximations for the differential equation and the boundary conditions. The mesh size will be h = 1/(N-1), resulting in N grid points and N unknowns.