function fordelning=psimulering(startp,P,n);
% fordelning=psimulering(startp,P,n);
% Funktion för att beräkna hur fördelningen utvecklas med tiden.
% Indata: startfordelning (startp)
%         overgangsmatris (P)
%         tidsintervall (n)
% Utdata: matris där varje kolumn innehåller sannolikheten för tillståndet
%         vid de olika tidpunkterna. Startvardet ar tidpunkten 1.

% Anders Svensson 941107
% Finn Lindgren 971113,981110,981125,991124

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

[a,b]=size(P);                                      % Antal tillstand

fordelning=zeros(n+1,a);                            % Utdata-matrisen
fordelning(1,:)=startp;

for j=2:n+1
  fordelning(j,:)=fordelning(j-1,:)*P;              % Uppdaterad slh
end;