function A=qbd(n,delta)
% QBD tale che A2=A3=r*ones(n)-r*eye(n)
% A1=A2+delta*eye(n)
% r= (1-delta)/(3*(n-1))
% Vale: drift=-delta; tanto piu' delta' e' piccolo,
% tanto piu' il problema e' "difficile"
A=zeros(n,n,3);
r=(1-delta)/(3*(n-1));
A(:,:,2)=r*ones(n)-r*eye(n);
A(:,:,3)=A(:,:,2);
A(:,:,1)=A(:,:,2)+delta*eye(n);
