Recently Fraser and Schoen showed that the solution of a certain extremal Steklov eigenvalue problem on a compact surface with boundary can be used to generate a free boundary minimal surface, i.e., a surface contained in the ball that has (i) zero mean curvature and (ii) meets the boundary of the ball orthogonally. In this paper, we develop numerical methods that use this connection to realize free boundary minimal surfaces. Namely, on a compact surface, Σ, with genus γ and b boundary components, we maximize σ j (Σ, g) L(∂Σ, g) over a class of smooth metrics, g, where σ j (Σ, g) is the j-th nonzero Steklov eigenvalue and L(∂Σ, g) is the length of ∂Σ.