I see.  The matrix I am building allows enforcement of neumann boundary condition in a linear system with '0' as the boundary derivative value, but here I have 'lambda * u(0)' as the derivative value.  How to fix?