Solving Linear System Eigen library - implicit smoothing C++

Could somebody point me into right direction how to solve this equation in C++ using Eigen library?

I need to make implicit smoothing for the mesh using this equation:

And I do not know how to create this matrix in Eigen:

Could somebody point me to right resource, I believe this kind of implicit smoothing using Eigen library was already done many times...