I am trying to plot the temperature distribution of a cake using the heat equation in 3D in MATLAB. The boundary conditions are
.
I am trying to plot the temperature distribution of a cake shaped as a rectangular prism on 2D slices. These slices are placed on the x-axis consecutively. (see below)

How can I visualize these slices? Is it also possible to make an animation of the temperature change over time for one of the slices? Thanks.