function y=BZ(x)

y = min( abs(x).^2, 0*x+1);
