function h=haria(p,geps)

N=size(p,1);
d=daria(p,[]);
h=6*min(-d*0.2+0.7*ones(N,1),ones(N,1));



