Numerical analysis for a three interacting species model with nonlocal and cross diffusion