function [a]=faia(n,gamma)
w=[2,-1];
w(n)=0;
L=toeplitz(w',w);
a=eye(n)+gamma*L;
end