The detection of flaws in steel pipes using Magnetic Flux Leakage (MFL) consists in detecting magnetic flux leaks outside the pipe, either with a magnetic sensor or with an induction coil, while the pipe is rotating. In the Vallourec group, many NDT units use MFL for testing ferromagnetic pipes. In order to improve the performances of flaw detection, CEA LIST and the Vallourec Research Aulnoye (VRA) group are collaborating on MFL modelling. The aim is to be able to perform parametric studies thanks to a fast 3D numerical model dedicated to MFL systems. A simplified 2D geometry has already been derived for the development of first simulation tools. When considering the B-H curve of ferromagnetic materials, the non-linear magnetostatic problem can be solved with the generalized boundary element method (BEMG), which comes to the evaluation of two equivalent scalar potentials: the surface charge density and the volume charge density. When applying the Galerkin method for the discretization of integral equations, the particularity of this numerical model lies in the implementation of high order basis functions for the interpolation of the scalar unknowns. This paper presents some first numerical results for the numerical validation of the semi-analytical model.