Although the study of surface segregation has a great technological importance, the work done in the field was for a long time largely restricted to experimental studies and the theoretical work was neglected. However, recent improvements in both first principles and semi-empirical methods are opening a new era for surface scientists. A method developed by Bozzolo, Ferrante, and Smith (BFS) is particularly suitable for complex systems and several aspects of the computational modeling of surfaces and segregation, including alloy surface segregation, structure and composition of alloy surfaces and the formation of surface alloys. In the following work I introduce the BFS method and apply it to model the Ni-Al alloy through a Monte-Carlo simulation. A comparison between my results and those results published by the group mentioned above was my goal. This thesis also includes a detailed explanation of the application of the BFS method to surfaces of multi-component metallic systems, beyond binary alloys.