An electric dipole is a system of two opposite point charges when their separation goes to zero and their charge goes to infinity in a way that the product of the charge and the separation remains finite.
How can we have a continuous volume charge distribution from such a collection of point charges?
In the article 'Electric dipole' at Knowino it is said that
The charge distribution is written in terms of Dirac delta functions: $$\rho (\mathbf{r})=q_1 \delta (\mathbf{r}-\mathbf{r}_1)+q_2 \delta (\mathbf{r}-\mathbf{r}_2)$$
Here $\mathbf{r}_1$ and $\mathbf{r}_2$ are the position vectors of $q_1$ and $q_2$ and $\mathbf{r}=\mathbf{r}_1-\mathbf{r}_2$. Please explain why do we need Dirac delta in describing dipolar charge distribution?