function [p,mu]=quadratic_i(x_i,B,A,J,d)
% hladanie smeru p, x_i pripustne riesenie, J vektor-index aktivnych
% ohraniceni v A (danych rovnicami)

[r,s]=size(A);q=length(J);

d_i=-B*x_i+d;A_it=zeros(q,s);
bb=zeros(q,1);
for i=1:q
A_it(i,:)=A(J(i),:);
end

%[x_r]=reduktion(B,A_it,d_i,b);
[x_r,mu]=reduktion_n(B,A_it,d_i,bb);
p=x_r;