Add `HeatEq`

As a playground to test different numerical flux methods for parabolic equations. Mainly to make sure that the fluxes we use in the AV method already work with parabolic equations.