% Foerst importerer du doedssannsynlighetene.
% Du gjoer det ved aa (instruksjonene er testet paa pc):
%  1) Laste ned filen "doedelighet.txt" fra kurssiden
%  2) Klikke paa "Import Data" i MATLAB
%  3) Velge filen "doedelighet.txt" 
%  4) Klikke paa "Import Selection" 
%     (importer dataene som "Column vectors")
% Merk at paa filen er doedssannsynlighetene gitt i promille.
%% 
% Plott doedsssannsynligheter per 1000 menn paa logaritmisk skala:
semilogy(alder,menn)

% Setter navn paa aksene:
xlabel('Alder')
ylabel('Dødsssannsynligheter per 1000') 
%% 

% Bestem de ettaarige doedssannsynlighetene (qx)
qx=menn/1000
%% 

% Bestem overlevelsessannsynlighetene 
% S(x)=P(X>x) for x=0,1,2,...,106:
Sx=cumprod(1-qx)
%% 

% Bestem den kumulative fordelingsfunksjonen F(x)=P(X<=x)=1-S(x)
% for x=0,1,2,...,106:
Fx=1-Sx
%% 

% Plott av den kumulative fordelingsfunksjonen:
stairs(alder,Fx)

% Setter navn paa aksene 
xlabel('Alder')
ylabel('F(x)') 
%% 

% Bestem punktsannsynligheten p(x)=P(X=x)=F(x)-F(x-1) 
% for x=0,1,2,...,106:
px=Fx-[0; Fx(1:106)]
%% 

%Plott av punktsannsynligheten:
bar(alder,px)

% Setter navn paa aksene 
xlabel('Alder')
ylabel('p(x)') 
%% 

% Beregn forventet levealder
EX=sum(alder.*px)
%%