title(['heat equation, t=', num2str(tau*n, '%4.2f')])
title(sprintf('heat equation, t=%4.2f', tau*n))