function REI = rotace(NI, dI)
%REI = baze(NI, dI)
%vrati rotaci, ktera natoci vektor NI do smeru [1 0 0] a vektor dI tak, ze
%jeho prumet do roviny yz ma smer [0 0 1]
 e1=normc(NI);
 e2=normc(cross(dI, e1));
 e3=cross(e1, e2);
 REI=[e1 e2 e3]';
end
