function [ A ] = genera_qbd_circ( k,m )
% Genera le matrici A di dimensione k*m che definiscono un QBD.
% k=dimensione a blocchi, m=dimensione dei blocchi
% A(:,:,3) e' 'circolante' a blocchi
% A(:,:,1) e' A3 trasposta + correzione random piccola
% A(:,:,2) e' random piccola
% Le it fun, anche la veloce con X0=I, arrancano
% 
 
A(:,:,3)=zeros(m*k);
 for i=1:k-1,A(i*m+1:(i+1)*m,(i-1)*m+1:i*m,3)=rand(m);
 end;
 A(1:m,(k-1)*m+1:k*m,3)=rand(m);
A(:,:,1)=A(:,:,3)'+1.e-10*rand(m*k);
A(:,:,2)=1.e-10*A(:,:,3);
v=sum(sum(A,3),2);
v=1./v;
for l=1:3
    A(:,:,l)=diag(v)*A(:,:,l);
end

end

