step-2: Added trilayer heterostructure as example.

Still computes a few quantities by hand, needs to use more of the now built-in Gubbins functionality.