A deterministic procedure for optimal approximation of arbitrary probability density functions by means of Dirac mixtures with equal weights is proposed. The optimality of this approximation is guaranteed by minimizing the distance of the approximation from the true density. For this purpose a distance measure is required, which is in general not well defined for Dirac mixtures. Hence, a key contribution is to compare the corresponding cumulative distribution functions. This paper concentrates on the simple and intuitive integral quadratic distance measure. For the special case of a Dirac mixture with equally weighted components, closed-form solutions for special types of densities like uniform and Gaussian densities are obtained. Closed-form solution of the given optimization problem is not possible in general. Hence, another key contribution is an efficient solution procedure for arbitrary true densities based on a homotopy continuation approach. In contrast to standard Monte Carlo techniques like particle filters that are based on random sampling, the proposed approach is deterministic and ensures an optimal approximation with ... mehrrespect to a given distance measure. In addition, the number of required components (particles) can easily be deduced by application of the proposed distance measure. The resulting approximations can be used as basis for recursive nonlinear filtering mechanism alternative to Monte Carlo methods.