% this code make a movie of a the surface as it changes

function make_movie(H,delta_t)

if nargin==1
    delta_t=.1;
end

m=min(min(H));
M=max(max(H));

figure(100)

for i=1:size(H,1)
    plot(H(i,:))
    axis([1,size(H,2),m,M])
    drawnow
%    pause(delta_t)
end

