function(y, npoints, k)
{
	#Esta funcion estima una densidad univariada usando el metodo k-nn
	#inputs:
	#y es la muestra tomada
	# npoints: es el numero de puntos donde se va a estimar la funcion de densidad
	# k : es el numero de vecinos mas cercanos
	#Output
	# la densidad estimada en los puntos npoints y su plot
	x <- seq(min(y) - 1, max(y) + 1, length = npoints)
	d1 <- abs(outer(x, y, "-"))
	m1 <- t(apply(d1, 1, sort))
	d2 <- 2 * m1[, k]
	fest <- k/(length(y) * d2)
	plot(x, fest, type = "l")
	fest
}