clear;              % clear memory
close all;          % close windows


P = [.55  .40  .05
     .35  .55  .10
     .20  .20  .60];

r = 0.04;
 
beta = 1/(1+r);


% cash flow

cf = [20
      10 
      -5];
 
v = [0
     0
     0];                           % arbitrarne zvolena pocatecni hodnotova funkce


 % iterations
convcrit = 0.0001;                  % konvergencni kriterium
diff = 1;                           % pocatecni hodnota rozdilu (jakakoliv hodnota vetsi nez convcrit)
iter = 0;                           % citac iteraci (jen pro zajimavost, kdy nam to zkonverguje)

while diff > convcrit
    Tv = cf + beta*P*v;
    diff = abs(v-Tv);               % rozdil mezi vektorama v pro dve iterace 
    v = Tv;                         % 
    iter = iter + 1;                % zvyseni hodnoty citace
end


disp('Pocet iteraci')        % pocet iteraci
disp(iter)

disp('Hodnota firmy [g n b]')   % hodnotova funkce
disp(v);

figure
plot(cf,v)
title('Hodnota firmy jako funkce dnesniho stavu (cash-flow)')
xlabel('cash flow v obdobi t')
ylabel('hodnota firmy v obdobi t')

%% to stejne pro jinou urokovou miru