function y=Besag(x)

y = abs(x).^1.2 ;

