function seep(startp,P,n)
% seep(startp,P,n) visar utvecklingen av fördelningen för
% en Markovkedja med övergångsmatris P och
% startfördelning startp, för n tidssteg.

% Finn Lindgren 1998-11-11, 1998-11-25, 1999-11-24

timedelay=0.6;

[ok,msg]=pchk(3,startp,P);
if ~ok
  error(msg)
end

f=psimulering(startp,P,n);
figure(gcf)
bar(f(1,:))
set(gcf,'doublebuffer','on')
a=axis;
axis([0 length(startp)+1 0 a(4)*1.1])
title(['Tidssteg = ' int2str(0)])
xlabel('Tillstånd')
ylabel('Sannolikhet')
a=[0 length(startp)+1 0 max(max(f(2:n+1,:)))*1.1];
drawnow
%Före ca. 5.3: tic;while toc<timedelay,drawnow;end
pause(timedelay)
for k=1:n
  bar(f(k+1,:))
  axis(a)
  title(['Tidssteg = ' int2str(k)])
  xlabel('Tillstånd')
  ylabel('Sannolikhet')
  drawnow
  %Före ca. 5.3: tic;while toc<timedelay,drawnow;end
  pause(timedelay)
end