%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%
%% FILE: esim_yyr.m
%%
%% Luvun 2 Yrjo Yhta-Reipas esimerkin ratkaisu funktiolla pmatriisi.
%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%% Palkkiomatriisi
R = [ 1   1   0   0;
      0   0   1   0;
      0   1   0   1];
%% Skenaarioiden todennakoisyydet (pystyvektori)
p = [0.05  0.45  0.25  0.25]';
      
%% Paatosongelma rakennemuuttujaan po. Kentat psaanto, pmat ja tn tarvitaan
po.pmat = R;
po.psaanto = "oa";
po.tn = p;
pr = pmatriisi(po)

%% Tiedoston ajo ja tuloste
%%>> esim_yyr
%%pr =
%%
%%  scalar structure containing the fields:
%%
%%    parvot =
%%
%%       0.50000
%%       0.25000
%%       0.70000
%%
%%    poptimi =  3 