Example from Luke Tierney (2003)'s introductory article on namespaces
mydnorm <- function(z) 1/sqrt(2 * pi) * exp(- z^2 / 2) mydnorm(-2:2)
## [1] 0.05399097 0.24197072 0.39894228 0.24197072 0.05399097
x <- seq(-5,5, length=100) all.equal(dnorm(x), mydnorm(x))
## [1] TRUE