HLsim=function(Nsim=100,m=15,n=15,avec=c(.2,.4,.6,.8,1,1.2),sig2=2){
k=length(avec)
out=matrix(rep(0,k*Nsim),ncol=k)
sig=sqrt(sig2)
for(i in 1:Nsim){
x=rnorm(m,0,sig)
y=rnorm(n,0,sig)
HL=median(outer(y,x,"-"))
pHL=abs(HL)<=avec
out[i,]=pHL
}
x=apply(out,2,mean)
x
}