核密度图

> par(mfrow=c(1,1))
> plot(d)
> polygon(d,col="green",border = "blue")
> rug(jitter(mtcars$mpg))