function [atil b]=Atilde(nv,a,n1,n2,i)
atil = a;       b = zeros(nv,1);
atil(n1,:) = 0; atil(n1,n1) = 1;
b(n2)=i;
return 
end