function res = residual(x,xold)

res=norm(x-xold,1);
