Scheme for solving Laplace's (Poisson's) equation
If we isolate on the left-hand side the unknown quantities
u11,
u12,
u21 and
u22, that is
the inner points not constrained by the boundary conditions, we can
rewrite the above equations as a matrix
A times an unknown vector
x, that is
Ax=b,
or in more detail
[4−1−10−140−1−104−10−1−14][u11u12u21u22]=[u01+u10−˜ρ11u13+u02−˜ρ12u31+u20−˜ρ21u32+u23−˜ρ22].