function y=beta_h(x,bl,br,alf,h) 
 
 y = br*heaviside(x,2*h) + bl*(1-heaviside(x,2*h));