function y=Asigma(x)
global A df n a delta 
y = ([-A  (a*a')/delta;speye(n) -A]-df*speye(2*n))\x;