A numerical method for reconstructing an impedance in a Sturm-Liouville operator from finitely many eigenvalues is investigated. The method constructs an impedance that has the given eigenvalues by finding a zero of a nonlinear finite dimensional map. A Newton scheme is investigated and numerical examples are considered.