We have performed high-resolution numerical simulations of supersonic slip surfaces to confirm and illuminate earlier analytic nonlinear stability calculations of such structures. This analytic work was in turn inspired by earlier computer simulations reported in Woodward (1985) and Woodward et al. (1987). In particular Artola & Majda (1987) examined the response of a supersonic slip surface to an incident train of small-amplitude nonlinear sound waves. They found analytic solutions which indicate that nonlinear resonance occurs at three angles of incidence which depend upon the Mach number of the relative motion. The two-dimensional simulations described here numerically solve this problem for a Mach-4 flow using the piecewise-parabolic method (Colella & Woodward 1984; Woodward & Colella 1984). The simulations show that sound waves incident at a predicted resonance angle excite nonlinear behaviour in the slip surface. At these angles the amplitude of the reflected waves is much greater than the incident wave amplitude (i.e. a shock forms). The observed resonance is fairly broad, but the resonance narrows as the strength of the incident waves is reduced. The nature of the nonlinear kink modes observed in the simulations is similar to that discussed by Artola & Majda. Most of the modes move in either direction with speeds near the predicted value. Speeds of other than this value are observed, but the disagreement is not serious in view of the strongly nonlinear behaviour seen in the simulations but not treated in the analytic work. The stationary modes seen in the analytic results are perhaps observed as transient structures. They may eventually dominate the flow at late times (Woodward et al. 1987). The role of the kink modes in the stability of slab jets is discussed, and it is argued that the stationary modes are more disruptive than the propagating modes.