In this study, a new method for producing initial conditions that are required for the Brownian dynamics simulation of dilute polymer suspension flows is proposed. For this purpose, the equilibrium probability distribution function is employed to generate an equilibrium ensemble of polymers. This approach is programmed and by using it, a polymer suspension in the inception of a simple shear flow is simulated. Also, the results of simulations in a similar flow configuration based on the conventional approach of generating initial conditions are presented. The excellent agreement between the results demonstrates the high accuracy of the proposed method for generating initial conditions. The main advantage of the proposed method is its low computational cost.