Let K be a number field, and S a finite set of places in K containing all infinite places. We present an implementation for solving the S-unit equation x + y = 1, x, y ∈ O × K,S in the computer algebra package SageMath. This paper outlines the mathematical basis for the implementation. We discuss and reference the results of extensive computations, including exponent bounds for solutions in many fields of small degree for small sets S. As an application, we prove an asymptotic version of Fermat's Last Theorem for totally real cubic number fields with bounded discriminant where 2 is totally ramified. In addition, we improve bounds for S-unit equation solutions that would allow enumeration of a certain class of genus 2 curves, and use the implementation to find all solutions to some cubic Ramanujan-Nagell equations.1 See also the recent translation [16] by Fuchs.