This article discusses an alternative method for estimating marginal probability densities of the solution to stochastic differential equations (SDEs). Two algorithms for calculating the numerical–statistical projection estimate for distributions of solutions to SDEs using Legendre polynomials are proposed. The root-mean-square error of this estimate is studied as a function of the projection expansion length, while the step of a numerical method for solving SDE and the sample size for expansion coefficients are fixed. The proposed technique is successfully verified on three one-dimensional SDEs that have stationary solutions with given one-dimensional distributions and exponential correlation functions. A comparative analysis of the proposed method for calculating the numerical–statistical projection estimate and the method for constructing the histogram is carried out.