The aim of this project is to implement a boundary element method program in order to solve the soil mechanics problem of a two-layered soil subject to a contact on his top surface. For this, I use Python programming language and follow the fictitious element method numerical procedures described in the literature. In the first section of this report, a brief explanation of the context and method is given. Then the methodology and model implementation are described gradually. Finally, the validation with several soil mechanics problems and the results are presented.