function y=glag_TVreg(x)


y=1./sqrt(1+x.^2);