function [v,m] = rekurse(P,R,n)
% funkce pocita maximalni stredni hodnotu celkoveho vynosu 
% a vektor optimalnich strategii po n krocích
% pro h ruznych strategii
% syntaxe: [v,m] = rekurse(P,R,n)
% vstupni parametry:
% P je tvaru [P1;P2;...Ph], kde Pi jsou matice prechodu
% R je tvaru [R1;R2;...Rh], kde Ri jsou matice vynosu
% n ... pocet kroku
% vystupni parametry:
% v ... vektor optimalnich strategii
% m ... vektor maximalnich strednich hodnot celkovych vynosu
% autor: Veronika Hajnova

%počet matic Pi

[a,b]=size(P);

c=b\a;

%vektory qi

for i=1:c
q(:,i)=diag(P(1+(i-1)*b:i*b,1:b)*R(1+(i-1)*b:i*b,1:b)');
end

%výpočet v1

[v(1,:),m(1,:)]=max(q');

%výpočet vi

for i=2:n
pom=q;
for p=1:b 
pom=pom+reshape(P(:,p),b,c).*v(i-1,p);
end
[v(i,:),m(i,:)]=max(pom');
end

v=v';
m=m';

end

