Data at t = 0:
Data at t = 0.2:
Solution at t = 0.5:
Solution in x-y plane at t = 0.5, z = 0.425 using view(zSlice): (Compare to Figure 11.2(a) in [LeV96])