The connection between extremal eigenvalues of the Steklov operator and free boundary minimal surfaces has been the subject of numerous works over the past decade. In this talk, we propose to study how this formalism can also enable the numerical approximation of these new geometric objects. We will be particularly interested in the implementation of numerical models allowing approximation in arbitrary genus.