subroutine mnnd(x,y,n,dminbig,dminav)
implicit double precision(a-h,o-z)
dimension x(n), y(n)

dminav = 0.d0
do i = 1,n {
	dmin = dminbig
	do j = 1,n {
		if(i!=j) {
			d = (x(i)-x(j))**2 + (y(i)-y(j))**2
			if(d < dmin) dmin = d
		}
	}
	dminav = dminav + sqrt(dmin)
}

dminav = dminav/n

return
end