MATLAB kommandoer til notatet om livsforsikring

 

 

% Vi leser f?rst inn d?delighetstabellen (jf. det f?rste obligatoriske oppgavesettet)

 

dod = textread('dod.txt', '', 'headerlines', 7);

alder = dod(:,1);

kvinner = dod(:,3);

qkvinner=kvinner./1000;

 

 

% Vi beregner s? punktsannsynligheten for gjenst?ende levetid og plotter den.

%(Disse beregningene er tilsvarende som beregningene i det f?rste obligatoriske

% oppgavesettet, men vi starter p? alder 30 ?r i stedet for 0 ?r.)

 

gjenst=alder(31:100)-30;

levesanns=cumprod(1-qkvinner(31:100));

levesanns=[1;levesanns(1:69)];

px=levesanns.*qkvinner(31:100);

scatter(gjenst,px)

 

 

% Vi beregner forventningsverdien gitt i formel (2) i notatet:

 

EhX=1000000*sum(px(1:35)./(1.03.^gjenst(1:35)))

 

 

% Vi beregner b?ken i formel (4) i notatet:

sum1=sum((1/1.03).^(gjenst(1:34)+1).*px(1:34));

EgXdivK=(1-sum1-(1/1.03)^(35)*sum(px(35:70)))/(1-1/1.03)

 

 

% Vi beregner til slutt den ?rlige premien:

 

premie=EhX/EgXdivK