axis([-60 40 -.1 1.1])
title(['density at time t = ',num2str(t)])