In this note we extend the sampling method to deal with the computation of eigenvalues of Sturm-Liouville systems. To do so, we use the 2-dimensional version of the sampling theorem to find a representation for eigencurves. These are easily approximated by computing a finite matrix of size N × N and we show that its truncation error is of order.