An unknown density function $f(x)$, its derivatives, and its characteristic function are estimated by means of Hermite functions $\{h_j\}$. The estimates use the partial sums of series of Hermite functions with coefficients $\hat{a}_{jn} = (1/n) \sum^n_{i=1} h_j(X_i)$ where $X_1\cdots X_n$ represent a sequence of i.i.d. random variables with the unknown density function $f$. The integrated mean square rate of convergence of the $p$th derivative of the estimate is $O(n^{(p/r) + (5/6r)-1})$. The same is true for the Fourier transform of the estimate to the characteristic function. Here the assumption is made that $(x - D)^r f \in L^2$ and $p < r$. Similar results are obtained for other conditions on $f$ and uniform mean square convergence.