We present a robust and accurate numerical algorithm for calculating energy-minimizing wavelengths of equilibrium states for diblock copolymers. The phase-field model for diblock copolymers is based on the nonlocal Cahn-Hilliard equation. The model consists of local and nonlocal terms associated with short- and long-range interactions, respectively. To solve the phase-field model efficiently and accurately, we use a linearly stabilized splitting-type scheme with a semi-implicit Fourier spectral method. To find energy-minimizing wavelengths of equilibrium states, we take two approaches. One is to obtain an equilibrium state from a long time simulation of the time-dependent partial differential equation with varying periodicity and choosing the energy-minimizing wavelength. The other is to directly solve the ordinary differential equation for the steady state. The results from these two methods are identical, which confirms the accuracy of the proposed algorithm. We also propose a simple and powerful formula: h = L*/m, where h is the space grid size, L* is the energy-minimizing wavelength, and m is the number of the numerical grid steps in one period of a wave. Two- and three-dimensional numerical results are presented validating the usefulness of the formula without trial and error or ad hoc processes. (C) 2014 Elsevier B.V. All rights reserved.