Nonparametric estimation of the copula function using Bernstein polynomials is studied. Convergence in the uniform topology is established. From the nonparametric Bernstein copula, the nonparametric Bernstein copula density is derived. It is shown that the nonparametric Bernstein copula density is closely related to the histogram estimator, but has the smoothing properties of kernel estimators. The optimal order of polynomial under the L2 norm is shown to be closely related to the inverse of the optimal smoothing factor for common nonparametric estimator. In order of magnitude, this estimator has variance equal to the square root of other common nonparametric estimators, e.g. kernel smoothers, but it is biased as a histogram estimator.