
axis([-1.3 1.3 -1.3 1.3])
daspect([1 1 1])
rybcolormap
caxis([.9 1.1])
title(sprintf('height at time t = %6.3f',t),'fontsize',20);
box on;
