function A=bidiag(a, b, n)
    
A=zeros(n); 
for i=1:n-1
   A(i, i)=a; 
   A(i+1, i)=b;
end
%A(n, n)=a;
end