%%% MMO : Cviceni 01
clear;

cilove_procento=5;


kurzy=load('akcie_2007.dat');

[roky,akcie]=size(kurzy);

vynosy=zeros(roky-1,akcie);

for i=1:roky-1
    vynosy(i,:)=(kurzy(i+1,:)-kurzy(i,:))./kurzy(i,:)*100;
end

prumerny_vynos=mean(vynosy);
riziko=cov(vynosy);

%% Objektivni funkce
H=riziko;
f=zeros(akcie,1);

%% Nerovnosti
A=-prumerny_vynos;
b=-cilove_procento;

%% Rovnosti
Aeq=ones(size(A));
beq=1.;

%% Omezeni
lb=zeros(akcie,1);
ub=ones(akcie,1);

[x,fval] = quadprog(H,f,A,b,Aeq,beq,lb,ub)

kurzy_skutecne=load('akcie_skutecne_2008.dat');
vynosy_skutecne=(kurzy_skutecne-kurzy(end,:))./kurzy(end,:)*100;

vynosy_skutecne*x