An algorithm is developed for the identification of Wiener systems, linear dynamic elements followed by static nonlinearities. In this case, the linear element is modeled using a recursive digital filter, while the static nonlinearity is represented by a spline of arbitrary but fixed degree. The primary contribution in this note is the use of variable knot splines, which allow for the use of splines with relatively few knot points, in the context of Wiener system identification. The model output is shown to be nonlinear in the filter parameters and in the knot points, but linear in the remaining spline parameters. Thus, a separable least squares algorithm is used to estimate the model parameters. Monte-Carlo simulations are used to compare the performance of the algorithm identifying models with linear and cubic spline nonlinearities, with a similar technique using polynomial nonlinearities.