Reindex the interior points, we get a low bandwidth matrix represents the nonzero relations between interior points

Then the only remaining thing is to solve a sparse linear system

Then the only remaining thing is to solve a sparse linear system

Graphs are supposed to fit the screen width

Refresh if graphs are not displayed correctly