Recursive calculation of the probability density function characterizing the state estimate of a nonlinear stochastic dynamic system in general cannot be performed exactly, since the type of the density changes with every processing step and the complexity increases. Hence, an approximation of the true density is required. Instead of using a single complicated approximating density, this paper is concerned with bounding the true density from below and from above by means of two simple densities. This provides a kind of guaranteed estimator with respect to the underlying true density, which requires a mechanism for ordering densities. Here, a partial ordering with respect to the cumulative distributions is employed. Based on this partial ordering, a modified Bayesian filter step is proposed, which recursively propagates lower and upper density bounds. A specific implementation for piecewise linear densities with finite support is used for demonstrating the performance of the new approach in simulations.