function A=condmat(n,kappa,vpisna)

% A=CONDMAT(n,kappa,vpisna) vrne nakljucno matriko n x n
% z obcutljivostjo kappa. Pri tem se generator nakljucnih
% stevil nastavi na stanje vpisna.

if nargin==2
    vpisna=0;
end

rand('state',vpisna)
tmp=rand(n);
[Q1,R]=qr(tmp);
tmp=rand(n);
[Q2,R]=qr(tmp);
D=diag(kappa.^(-(0:n-1)/(n-1)));
A=Q1*D*Q2;