% Rydd opp
clc
clear
close
format shortG

% Stokastisk matrise for endringer i Stortingsvalg 2013 -> 2017.
% Tatt fra SSB, men har jukset for å få kolonnesummene til å bli 1 ved
% å legge til 3, 1, 1, 7, 2, 3, 2, 3, 2, 8, 8, 13 på diagonalen.

P = (1/100)*[80+3 3 7 3 23 30 1 23 18 27;     % AP
              1 66+1 7 0 3 1 1 0 0 20;        % FrP
              4 22 77+1 6 8 4 39 9 3 22;      % H
              1 3 3 79+7 5 1 2 0 2 2;         % KrF
              1 1 0 0 43+2 0 0 1 2 0;         % SP
              4 1 0 0 2 47+3 1 7 20 2;        % SV
              5 1 3 1 8 8 53+2 23 2 0;        % V
              1 1 0 0 4 5 1 28+8 7 4;         % MDG
              0 0 1 0 1 1 0 0 35+8 0;         % R         
              0 1 1 4 1 0 0 1 3 10+13         % Andre
]

pause

% Sannsynlighetsvektor for Stortingsvalg 2017.
x = (1/98)*[27;	
            15;	
            25;
            4;	
            10;	
            6;	
            4;	
            3;	
            2;	
            2]
pause
       
% Valgresultat 2021 med endringer
Valgresultat = P*x
Endringer = P*x-x
pause

% Egenverdier og egenvektorer
[Egenvektorer Egenverdier] = eig(P)
pause
E1 = Egenvektorer(:,1,:);           % Ta ut første egenvektor.
Likevekt = 1/sum(E1) * E1
pause

% Regulær?
P2 = P*P