[en] A computer code is described which calculates the initial fast ion deposition for neutral beam injection into noncircular tokamak plasmas. The beam deposition is calculated by Monte Carlo techniques as a function of the poloidal flux function psi. The flux function and plasma parameters are obtained from the numerical solution of the axisymmetric MHD equilibrium equation. Bounce averages are performed over the initial ion orbits, and the deposition is weighted by the fraction of time the ion spends in each interval of the psi grid. The divergence, focusing, and shape of the neutral beam are modeled. The beam geometry part of this problem was solved by G. G. Lister, D. E. Post, and R. Goldston of the Princeton Plasma Physics Laboratory in a code called FREYA. The code described here is a modification of FREYA. Comparisons of the beam depositions calculated with bounce averaging are made with the ion birth depositions for the PDX and ISX-B tokamaks