% Foerst importerer du doedssannsynlighetene.
% Du gjoer det ved aa 
%  1) klikke paa "Import Data" 
%  2) lese inn filen "doedssanns.txt"
%  3) klikke paa "Import Selection" 
% 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('Doedsssannsynligheter per 1000') 

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

plot(alder,qx)
xlabel('Alder')
ylabel('Dødsssannsynligheter') 
set(gca,'YScale','log')

% 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
p = [Fx(1); diff(Fx)]
bar([0:106], p)
terning.BarWidth=1