We show, using quasi-exact numerical simulations, that Anderson localization in a disordered one-dimensional potential survives in the presence of attractive interaction between particles. The localization length of the particles' center of mass—computed analytically for weak disorder—is in good agreement with the quasi-exact numerical observations using the time evolving block decimation algorithm. Our approach allows for simulation of the entire experiment including the final measurement of all atom positions.