function [H,C]=init_hyperstate(regulator_parameters, apriori) n_grid=regulator_parameters.n_grid; H=zeros(regulator_parameters.dim,n_grid); rozsah=apriori.Hmax-apriori.Hmin; a=0; for i=1:size(regulator_parameters.H_tilde,1) for j=1:size(regulator_parameters.H_tilde,2) if (regulator_parameters.H_tilde(i,j)) a=a+1; H(a,:)=ones(1,n_grid)*apriori.Hmin(i,j)+(0:n_grid-1)*rozsah(i,j)/(n_grid-1); end end end for i=1: end