TASK 3: Initial version
TASK 3: Very basic version works

Can solve thermal conductivity w/o sources within system with rectangular
boundaries.

NB! Local-1D method here.