n <- 100              
normal.sample <- rnorm(n)               # Generate sample from N(0,1)
weights <- rep(1,n)
density.estimate <- density(normal.sample, weights=weights)
                                        # Compute density estimate
plot(density.estimate)                  # Plot density estimate