A novel three-dimensional wide-angle beam propagation method based on the split-step fast Fourier transform is developed. The formulation is based on the three-dimensional Helmholtz wave equation. Each propagation step is performed by utilizing both the FFT and split-step scheme. The solution of Helmholtz wave equation does not make the slowly varying envelope and one-way propagation approximations. To validate the efficiency and accuracy, numerical results for a propagation beam in a tilted step-index optical waveguide are compared with other beam propagation algorithms.